

07398-^027 -R0-00 


RAD IO/OPTICAL /STRAPDOWN I3STEICTI-A.L 
C3-XJID0 A.3STOE STUDY FOK A^IDVA^ZSTOEID 

kick; stage applications 

FINAL REPORT 

"VolTJurin.© II— ID et ailed Study Results (Tasks 111 and IV) 

November 1968 


PART I DERIVATION OF GUIDANCE AND CONTROL 

FUNCTIONAL AND PERFORMANCE REQUIREMENTS 


Contract No. NAS 12-141 
Prepared for 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
ELECTRONICS RESEARCH CENTER 
Cambridge, Massachusetts 



ONE SPACE PARK • REDONDO BEACH, CALIFORNIA 



Reproduced by 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

US Department of Commerce 
- Sprtngfield, YA. 22151 





07398 -6027-RO-00 


PP, Al!D IO / OPTICAL / STR, AFD O W2ST INERTIAL 
OTJID^lSrOE STUDY FOR ADVANCED 
KICK STAGE APPLICATIONS 

FIISLAL REPORT 

Volume II — IDetaj-leci Study Results (Tasks III and IV) 

November 1968 


PART I DERIVATION OF GUIDANCE AND CONTROL 

FUNCTIONAL AND PERFORMANCE REQUIREMENTS 


Contract No. NAS 12-141 
Prepared for 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
ELECTRONICS RESEARCH CENTER 
Cambridge, Massachusetts 


U% Tf# 


ini c.-STa/" 


TRW 

SYSTEMS GROUP 


ONE SPAfcE PARK • REDONDO BEACH, CALIFORNIA 



FOREWORD 


This final report documents the results of the work 
accomplished under Tasks III and IV of a study of Radio/ 
Optical/Strapdown Inertial Guidance Systems for future 
unmanned space missions, conducted by TRW Systems for 
the NASA Electronics Research Center under Contract 
NAS 12-141. This effort expands and extends the work 
accomplished previously under Tasks I and II of the same 
contract. 

Volume I summarizes both the results of the study 
and recommendations reached, including those developed 
under Tasks X and II. Volume II documents the detailed 
study results for Tasks III and IV. 



CONTENTS 


VOLUME II - PART I 

1. INTRODUCTION 1 

1. 1 Study Objectives 1 

1. 2 Study Implementation 4 

i. 2. 1 Derivation of Guidance and Control 

Functional and Performance Requirements 4 

1. 2. 2 Preliminary Modular Design 8 

1. 3 Definition of Terms 11 

4. 3. 1 Missions 12 

1. 3. 2 Vehicle Terms 13 

1. 3. 3 Mission Events 14 

1.3.4 Trajectory Terms 15 

1. 3. 5 Coordinate Systems 17 

2. MISSION CHARACTERISTICS 20 

2. 1 Introduction 20 

2. 2 Earth Low-Altitude Polar -Orbit Mission 20 

2.3 Mars Orbiter Missions 23 

2.4 Solar Probe (With Jupiter Assist) Missions 25 

3. GUIDANCE AND CONTROL SYSTEM CONCEPTUAL DESIGNS 35 

3. 1 Introduction 35 

3. 2 Guidance and Control Requirements ' 36 

3. 2. 1 Guidance and Control System Operational 

Sequences 36 

3. 2. 2 Control System Functional Requirements 39 

3. 3 Guidance and Control System Conceptual Designs 42 

3.4 Near-Earth Polar-Orbit Mission 45 

3.5 Earth-Synchronous Satellite Mission 48 

3. 6 Lunar Orbiter Mission 50 

3.7 Mars Orbiter Mission 50 

3.8 Solar Probe Missions 53 

4. GUIDANCE AND NAVIGATION PERFORMANCE ANALYSES 56 

4. 1 Introduction and Summary 56 

4.2 Powered Performance Analyses of the Near-Earth 

Polar -Orbit Mis sion 58 

4. 3 Analysis of the Sun Sighting, Perigee Burn-Time Update 62 

iii 



CONTENTS (Continued) 



4. 3. 1 

Visibility Constraints 

64 


4. 3. 2 

Accuracy Considerations 

64 

4.4 

Powered Performance Analyses of the Mars Missions 

69 


4.4. 1 

Injection Errors 

69 


4. 4. 2 

Target Misses Due to Injection Errors 

72 


4.4. 3 

Corrective Velocity Requirements 

74 


4.4.4 

Corrected Target Miss 

74 

4.5 

Summary of Powered Performance Analyses of the 



Jupiter 

Missions 

77 

4. 6 

Mars Approach Analyses 

79 


4. 6. 1 

Summary of SVEAD Results 

89 


4. 6.2 

Conversion of SVEAD Results to Approach 
Condition Errors 

89 


4.6.3 

Conclusions from Approach Navigation 
Performance Analyses 

98 


4. 6.4 

Effects of Approach Orbit Determination and 
Execution Error s on In-Plane Performance 

103 

. 

4. 6. 5 

Effect of Approach Orbit Determination Errors 
on Out-of-Plane Performance 

113 

4.7 

Mars and Lunar Orbit Determination 

114 


4. 7. 1 

Mars Orbit Determination from DSIF Tracking 
Data 

114 


4. 7. 2 

Lunar Orbit Determination From MSFN Tracking 
Data 

115 

5. CONTROL SYSTEM PERFORMANCE ANALYSES 

124 

5. 1 

Introduction and Summary 

124 

5. 2 

Control System Stability and Bending Mode 
Transmissibility Coefficients 

125 

5. 3 

Stability Analysis of the Atlas/ Centaur Vehicle 

132 

5.4 

Stability Analysis of the Saturn IB/ Centaur Vehicle 

163 

5. 5 

Analysis of the Voyager Spacecraft 

172 


5. 5. 1 

Powered Flight Stability Analysis 

172 


5.5.2 

Voyager Spacecraft Thrust Vector Pointing 
Error Analysis 

174 

APPENDICES 



A 

CONTROL SYSTEM STABILITY CRITERIA 

184 


1 . 

Rigid -Body Stability Margins 

184 


2. 

Propellant-Slosh Stability Margins 

184 


3. 

Structural Bending Mode Stability 
Margins 

184 


IV 



CONTENTS (Continued) 


B VEHICLE DYNAMICS AND CONTROL SYSTEM 

HOMOGENEOUS EQUATION MATRIX 188 

C ATLAS/ CENTAUR CONTROL SYSTEM DATA 196 


1 . 

Introduction 

196 

2. 

Aerodynamic, Mass Properties, and 
Trajectory-Related Controls Data 

201 

3. 

Propellant Sloshing Data 

202 

4. 

Vehicle Bending Data 

204 

5. 

Thrust Vector Control Data 

205 

SATURN IB/ CENTAUR CONTROL SYSTEM DATA 

211 

1 . 

Introduction 

211 

2. 

Aerodynamic Mass Properties, and 
Trajectory-Related Controls Data 

213 

3. 

Propellant Sloshing Data 

213 

4. 

Vehicle Bending Data 

213 

5. 

Thrust Vector Control Data 

226 

SATURN V (VOYAGER PAYLOAD) CONTROL SYSTEM 
DATA 

228 

1 . 

Introduction 

228 

2. 

Propellant Sloshing Data 

230 

3. 

Saturn V Thrust Vector Control 
System Data 

232 

4. 

Voyager Spacecraft Data 

249 


v 



CONTENTS (Continued) 
VOLUME II - PART II 


6. SYSTEM CONFIGURATION AND INTERFACES 1 

6. 1 Modularity Concept 1 

6. 2 Vehicle Interfaces and Mechanical Mounting 

Considerations 8 

6. 3 Guidance Equipment Mechanical Interface and 

Packaging Considerations 12 

6. 3. 1 Navigation Platform Subassembly 12 

6.3.2 Sensor Alignment 13 

6. 3. 3 Electronics Packaging . 14 

6.4 Control System Interfaces and Modularity 17 

6. 5 Thermal Control Considerations 20 

6.5. 1 Summary of Conclusions and Recommendations 2.0 

6.5.2 Vehicle Guidance /Equipment Thermal Interface 23 
Considerations 

6.5.3 Sensor Package Thermal Design Considerations 27 

7. ONBOARD COMPUTATIONAL ELEMENTS 32 

7. 1 Summary 32 

7. 2 Computer Sizing 32 

7. 2. 1 Introduction 32 

7.2.2 Major Functions to be Performed by the 

Computer 36 

7.2.3 Storage and Timing Requirement Estimates 42 

7.3 Guidance Computer Interface 45 

7.3.1 General Description 45 

7.3.2 Gomputer Description 46 

7.3.3 CIU Description 49 

7.3.4 IMU Interface 51 

7.3.5 Booster Rate Gyro Interface 52 

7.3.6 Electro -Optical Sensor Interface 52 

7.3.7 Control Electronics Interface 55 

7. 3.8 Vehicle and Spacecraft Discrete Interface 55 

7.‘3.9 Command Data Interface 57 

7.3.10 Telemetry 57 

7.3.11 CIU Interface Control 58 

7.3.12 CIU Operational Description 60 

7.3.13 Telemetry Units 62 

7. 3. 14 CIU Reliability and Physical Characteristics 63 

7.4 Guidance Computer Reliability Assessment 64 


vi 



CONTENTS (Continued) 


8. CONTROL SUBSYSTEM DESIGN 69 

8. 1 Introduction and Summary 69 

8. 2 Control Subsystem Interface Considerations 72 

8.2.1 Atlas SLV3C/ Centaur/ Spacecraft 72 

8.2.2 Saturn IB/ Centaur/Spacecraft 79 

8. 2, 3 Saturn V/Spacecraft 84 

8. 3 Digital Control Equations 87 

9. ELECTRO-OPTICAL SENSORS 99 

9. 1 Introduction and Summary 99 

9.2 Selected Electro -Optical Sensors and Operational 

Sequences 100 

9.2.1 Earth Polar Orbit 100 

9. 2. 2 Earth -Synchronous Orbit (Direct Ascent) 100 

9.2.3 Earth-Synchronous Orbit (Parking 

Orbit Injection 102 

9.2.4 Lunar Orbiter 102 

9.2.5 Mars Orbiter 103 

9.2.6 Jupiter Flyby 103 

9. 3 Summary of Sensor Accuracies 104 

9.4 Sun Sensors 104 

9.4.1 Digital Solar-Aspect Sensors 104 

9.4.2 Coarse and Fine Sun Sensors 119 

9.5 Earth Sensor 134 

9.5. 1 Rationale for Selection 136 

9.5.2 Background and Present Status 136 

9.5.3 General System Description 138 

9.5.4 Reliability 154 

9.5.5 Earth-Sensor Accuracy 160 

9.5.6 Candidate Earth-Sensor Subsystem • 

Configuration 170 

9.5.7 Mechanical Configuration and Mounting 

Characteristics 173 

9.5.8 Preliminary Specification-Earth 

Horizon Sensor 173 

9.6 Canopus Trackers 184 

9. 6. 1 Rationale for Selection 184 

9. 6.2 Selected Canopus Tracker for Lunar and 

Interplanetary Missions 188 

9.6.3 Canopus Tracker for Mars Approach 

Guidance 207 


vii 



CONTENTS (Continued) 

Page 

9.7 Planet-Approach Sensor 237 

9.7.1 Rationale for Selection 237 

9.7.2 Optical Field-of-View and Pointing 

Angle Requirements 239 

9.7.3 General System Description 242 

9.7.4 Detailed Characteristics and Estimated 

Accuracy 246 

9.7.5 Mechanical Configuration 254 

9.7.6 Reliability 254 

9.7.7 Performance Specification 257 

10. INERTIAL REFERENCE UNIT 263 

10. 1 Functional Description and Performance 

Characteristics 263 

10.2 TG-166 IRU Design and Packaging Concept 277 

10.2. 1 IRU Instrument-Block Subassembly 278 

10.2.2 Electronic Packaging Concept 280 

10.2.3 Thermal Design Considerations 281 

11. SUMMARY OF PERFORMANCE CHARACTERISTICS 

OF MODULAR DESIGN 283 

11.1 Trajectory Accuracy and Fuel Required 
For Correction of Guidance, Navigation, 

and Control Errors 283 

11.2 Summary of System Reliability, Weight, and 

Power Requirements 285 

APPENDIX F THERMAL ENVIRONMENT AND THERMAL 285 

CONTROL TECHNIQUES 

1. Introduction 288 

2. Thermal Control for Ascent and 289 

Orbit Injection Phases 

3. Thermal Control Techniques for 292 

Extended Earth Orbiting Phases 

4. Thermal Control Techniques for Translunar 294 

and Interplanetary Phases 


viii 



ILLUSTRATIONS 


2-i. Ecliptic Projection of Sample 1975 Type I Mars Mission, 

Showing Relative Heliocentric Positions of Earth, 

Vehicle, and Mars . 27 

2-2. Ecliptic Projection of Sample 1975 Type II Mars Mission, 
Showing Relative Heliocentric Positions of Earth, 

Vehicle, and Mars 28 

2-3. Time Histories of Heliocentric Orientation Angles and 

Distances for 1975 Type I Mars Mission 29 

2-4. Time Histories of Heliocentric Orientation Angles and 

Distances for 1975 Type II Mars Mission 30 

2- 5. Encounter Geometry 32 

3- 1. Composite Equipment Configuration 43 

3-2. Control System Signal and Equations Flow Diagram 46 

3-3. Basic Conceptual Design Configuration for the Near- 

Earth Polar -Orbit Mission 47 

3-4. Basic Conceptual Design- Configuration for the Earth- 

Synchronous Satellite Mission 49 

3-5. Basic Conceptual Design Configuration for the Lunar 

Orbiter Mission 51 

3-6. Basic Conceptual Design Configuration for the Mars 

Orbiter Mission 52 

3- 7. Basic Conceptual Design Configuration for the Jupiter 

Missions 54 

4- 1. Corrective Velocity Requirements for Atlas/ Burner II 

Near-Earth Polar-Orbit Mission 61 

4-2. Sun Sighting Time Update Technique for Multi-Parking 

Orbit Synchronous Satellite Mission 64 

4-3. Combinations of Q and t that Satisfy the Visibility 

Constraints 65 

4-4. Normal Component Sensitivity to Initial Azimuth 

Uncertainty 73 

4-5. Uncorrected Miss Ellipse for Mars Trajectory, Type I 75 

4-6. Uncorrected Miss Ellipse for Mars Trajectory, Type II 75 


IX 



4-7. 

4-8. 

4-9. 

4-10. 

4-11. 

4-12. 

4-13. 

4-14. 

4-15. 

4-16. 

4-17. 

4-18. 

4-19. 

4-20. 

4-21. 

4-22. 

4-23. 

4-24. 

4-25. 

4-26. 


ILLUSTRATIONS (Continued) 

Miss Ellipse after First Midcourse Correction 
for Mars Trajectory, Type I 

Miss Ellipse after First Midcourse Correction 
for Mars Trajectory, Type II 

Uncorrected Miss Ellipse for Jupiter, Solar 
Probe Trajectory 

Uncorrected Miss Ellipse for Jupiter, Out of 
Ecliptic Trajectory 

Miss Ellipse after First Midcourse Correction 
for Jupiter, Out of Ecliptic Trajectory 

Miss Ellipse after First Midcourse Correction 
for Jupiter, Solar Probe Trajectory 

Measurement Angles Versus Time to Encounter, 

1975 Type I and II Trajectories 

Optical Angle Measurements 

Position Uncertainty Versus Time to Encounter 

Position Uncertainty Versus Time to Encounter 

Position Uncertainty Versus Time to Encounter 

Velocity Uncertainty Versus Time to Encounter 

Velocity Uncertainty Versus Time to Encounter 

Velocity Uncertainty Versus Time to Encounter 

Impact Parameter Standard Deviation, 1975 Type I and 
II Trajectories 

Approach Velocity Standard Deviation, 1975 Type I and 
II Trajectories 


78 

78 

80 

80 

82 

82 

83 

84 

90 

91 

92 

93 

94 

95 

99 

99 


Orbital Inclination Standard Deviation, 1975 Type, I and 
II Trajectories 

Ascending Node Longitude Standard Deviation, 1975 Type 

I and II Trajectories 

Argument of Periapsis .Standard Deviation, 1975 Type I 
and II Trajectories 

Time -to -Encounter Standard Deviation, 1975 Type I and 

II Trajectories 


100 

100 

101 

101 


x 



4-30. ±lcr Spread in Parameters of Deboosted Orbit Due to 

Uncertainty in Approach Impact Parameter, B 112 

4-31. Position Uncertainties versus Time for Spacecraft in 

Mars Orbit 116 

4-32. Velocity Uncertainties Versus Time for Spacecraft in 

Mars Orbit 117 

4-33. Uncertainties in Orbital Parameters Versus Time for 

Spacecraft in Mars Orbit ■ 118 

4-34. Uncertainties in Orbital Orientation Parameters Versus 

Time for Spacecraft in Mars Orbit 119 

4- 35. Comparison of L03 Trajectory Data with Trajectory 

Computed Using Triaxial Gravity Model 122 

5- 1. Typical Gain- Phase Diagram With Margins 126 

5-2. Effects of Bending Mode Transmissibility Upon Stability 

Margins Without Compensation 126 

5-3, First Stage Booster Bending Transmissibility Coefficients 128 

5-4. Atlas/Centaur Higher Mode Transmissibility Coefficients 129 

5-5. Saturn V/ Voyager First Bending Mode Transmissibility 

Coefficients Without Filtering 130 

5-6. Saturn V/ Voyager Second Bending Mode Transmissibility 

Coefficients Without Filtering 131 

5-7. Saturn IB/ Apollo CSM First-Stage Booster Bending 

Transmissibility Coefficients 133 

5-8. Atlas S-LV3/ Burner First-Stage Booster Bending 

Transmissibility Coefficients 133 

5-9. Gain-Phase Plot Atlas/ Centaur AC=13 T=0 Analog 

Autopilot 134 

5-10. Gain-Phase Plot Atlas/Centaur AC- 13 T=40 Analog 

Autopilot 135 

5-11. Gain- Phase Plot Atlas/Centaur AC- 13 Max Q T=80 

Analog Autopilot 136 

5-12. Gain-Phase Plot Atlas/Centaur AC- 13 T=i22 Analog 

Autopilot 137 


xi 



ILLUSTRATIONS (Continued) 


5-13. Gain-Phase Plot Atlas/ Centaur AC- 13 T=154 Analog 

Autopilot 138 

5-14. Atlas/ Centaur AC-13 Rigid Body Gain Margins, Analog 

Autopilot, and Complex Poles at 15. 5 Rad/ Sec 140' 

5-15. Gain-Phase Plot Atlas/Centaur AC- 13 T=0 Analog Autopilot 

WF = 11 Rad/ Sec 141 

5-16. Gain- Phase Plot Atlas/ Centaur AC -13 T=80 Analog Auto- 
pilot WF = 11 Rad/ Sec 142 

5-17. Gain- Phase Plot Atlas/Centaur AC- 13 T=154 Analog 

Autopilot WF = 11 Rad/ Sec 143 

5-18. Atlas/ Centaur AC- 13 Rigid Body. Gain Margins, Analog 

Autopilot, and Complex Poles changed at T = 70 144 

5-19. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T-0 Atlas Rate Gyros Run 1 146 

5-20. Gain-Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=40 Atlas Rate Gyros 147 

5-21. Gain- Phase Plot Atlas/Centaur AC-13 Digital Autopilot 

Max Q T=80 Atlas Rate Gyros ■ 148 

5-22. Gain-Phase Plot Atlas/Centaur AC- 13 Digital Autopilot 

T=112 Atlas Rate Gyros Run 149 

5-23. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=154 Atlas Rate Gyros 15 0 

5-24. Atlas/Centaur AC-13 Body Gain Margins, ROI Guidance 
System With Atlas Rate Gyros, and Complex Poles at 
V = 0.3 Rad/Sec 151 

5-25. Gain-Phase Plot Atlas/Centaur AC- 13 Digital Autopilot 

T=80 Atlas Rate Gyros WF = 15. 5 + Notch 152 

5-26. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=80 Atlas Rate Gyros Notch 8. 8 153 

5-27. Gain-Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=0 Spacecraft Sensors WF =11 154 

5-28. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=80 Spacecraft Sensors WF = 11 .155 

5-29. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=80 Spacecraft Sensors Compensation 1 157 

xii 



ILLUSTRATIONS (Continued) 


5-30. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=80 Spacecraft Sensors Compensation 2 158 

5-31. Atlas/Centaur AC-13 First Mode Antinode Location versus 

Flight Time 159 

5-32. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=0 Atlas Rate Gyros at Station 708 160 

5-33. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=80 Atlas Rate Gyros at Station 708 161 

5-34. Gain- Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 

T=B0 Atlas Rate Gyros at Station 708 162 

5-35. Gain- Phase Plot Saturn IB/Centaur at T-0 With Com- 
pensation 1 164 

5-36. Gain- Phase Plot Saturn IB/Centaur Max Q T=80 With 

Compensation 1 165 

5-37. Gain- Phase Plot Saturn IB/ Centaur at T=B0 With Com- 
pensation 1 166 

5-38. Gain- Phase Plot Saturn IB/ Centaur T=0 Position Sensor 

at Spacecraft Compensation 1 167 

5-39. Gain- Phase Plot Saturn IB/ Centaur T=80 Position Sensors 

at Spacecraft Compensation 1 168 

5-40. Gain-Phase Plot Saturn IB/ Centaur T=B0 Position Sensors 

at Spacecraft Compensation 1 169 

5-41a. Saturn IB/ Centaur Rigid- Body Gain Margins 171 

5-41b. Rigid-Body Phase Margin versus Flight Time With 
Position Gain as a Parameter for the Saturn IB/ 

Centaur Stage I Digital Control System 172 

5-42. Gain-Phase Plot, MCC, Capsule, Integral 

Attitude Feedback, Kj = 0. 125 175 

5-43. Gain-Phase Plot, Start MOI, With Capsule, Integral 

Attitude Feedback, Kj = 0. 125 176 

5-44. Gain- Phase Plot, End MOI, With Capsule, Integral 

Attitude Feedback, Kj = 0.125 177 

5-45. Gain-Phase Plot, MOT, With Capsule, Integral 

Attitude Feedback, Kj = 0.125 178 


xiii 



TABLES (Continued) 


4-IX 

Ninety-Five Percent AV Midcourse (5 Days) 
Requirements for 1975 Mars Missions 

76 

4-X 

Saturn IB/Centaur Jupiter Missions Injection Errors 
(RTN Coordinates) 

77 

4-XI 

Ninety-Five Percent AV Midcourse (5 Days) 
Requirements for the Two 1971 Jupiter Missions 
(TG-266 System) 

81 

4- XII 

Radio/ Optical/Inertial Error Model Mars Mission 

85 

4-xni 

Optical Error Model A 

86 

4- XIV 

Optical Error Model B 

87 

4-XV 

Optical Error Model C 

88 

4-XVI 

Normalized Covariance Matrices and 1-cr Component 
Values at 2 Hr Before Encounter 

102 

4- XVII 

Ninety-Five Percent Orbit Trim AV Requirements 
(M/Sec) 

ill 

4-xvm 

Maximum Errors in Deboosted Orbit Parameters 

113 

4 -XIX 

Orbital Parameters for Mars Orbit (Martian 
Equatorial Coordinate System) 

114 

4 -XX 

Lunar Orbit Navigation Uncertainties 

120 

5-1 

Saturn IB/Centaur Control System Gains and Relative 
Stability Margins — Aft S-IVB Rate Gyros 

170 

5 -II 

Saturn IB /Centaur Control System Gains and Relative 
Stability Margins — Spacecraft Sensors 

170 

5 -III 

Saturn IB/Centaur Closed Loop Slosh Poles, 
S-IB Phase 

173 


xvi 



ILLUSTRATIONS (Continued) 


5-46. e ,p Without Integrators 179 

5-47. Vrj, Without Integrators 179 

5-48. e With Attitude Feedback Integration 180 

5-49. Vrp With Attitude Feedback Integration 180 

5-50. € ^ With Attitude and Engine Angle Feedback Integration 181 

5-51. Vrj-, With Attitude and Engine Angle Feedback Integration 181 


xiv 



TABLES 


2-1 

Radio/Optical/Strapdown Inertial Task III Mission 
Summaries 

' 21 

2 -II 

Atlas SLV- 3A/Burner II Sequence of Events 

22 

2 -III 

Saturn V Launch and Injection Trajectory 
Characteristics 

24 

2 -IV 

Saturn V Sequence of Events (Type II Transfer 
T rajectory) 

25 

Z-V 

Characteristics of 1975 Earth-Mars Trajectories 

26 

2- VI 

Saturn IB/Centaur Launch and Injection Trajectory 
Characte ristics 

31 

2 -VII 

Saturn IB /Centaur Sequence of Events 

31 

2 -VIII 

Characteristics of 1972 Jupiter Probes 

33 

3-1 

Equipment Utilization 

-44 

4-1 

Error Models Used for Strapdown Inertial Guidance 
Performance Analysis 

57 

4 -II 

Atlas /Burner II Near-Earth Polar-Orbit Injection 
Errors 

58 

4-HI 

Summary of Atlas /Burner II Near-Earth Polar- 
Orbit Mission Error Analysis Results 

59 

4-IV 

Summary of Atlas /Burner II Near-Earth Polar- 
Orbit Mission Error Analysis Results 

60 

4- V 

Summary of Synchronous Mission Performance 
Analysis Results 

63 

4- VI 

Summary of Long Coast Synchronous Satellite 
Mission Error Analysis Results (TG-166 System) 

70 

4 -VII 

Summary of Long Coast Synchronous Satellite 
Mission Error Analysis Results (TG-266 System) 

71 

4 -VIII 

Saturn V Mars Mission Injection Errors (Initial 
Azimuth Alignment Error = 20 Arc Sec) 

72 


xv 



i. INTRODUCTION 


This TRW Systems final report documents the detailed results of 
the work accomplished for the NASA Electronics Research Center under 
Tasks III and IV of a study of "Radio/Optical/Strapdown Inertial Guidance 
Systems" for application to future unmanned space missions. This effort 
extends and refines the work previously carried out and documented under 
Tasks I and II (see Ref. -1-i). 

I.’ 1 STUDY OBJECTIVES 

The primary objective of the study was to evaluate the feasibility 
of the " integrated modular design" concept for the guidance and control 
of launch vehicles and spacecraft for specified NASA unmanned space 
missions by means of analysis and design of a responsive system. 
Detailed study objectives were to: 

« Establish the guidance and control requirements for a 
selected group of future NASA space missions. 

• Investigate possible guidance concepts based on the 

appropriate use of radio, inertial, and optical techniques, 
with .the further objective of establishing the -functional 
role, the capabilities, limitations, and constraints of 
■’ each of these elements in the overall guidance system 
concept. 

-Define feasible radio/ optical/strapdown inertial navigation, 
guidance and control system "conceptual designs." 

® Perform analyses to establish the feasibility (per- 
formance) of the selected conceptual designs and to 
'establish the significant performance characteris- 
- tics of each component and subsystem. 

» Perform a "Preliminary modular design" of the radio/ 

. - optical/ strapdown inertial system meeting the com- 
posite requirements of all the missions considered, 
configured so that specific components may be 
interchangeably combined into given operational 
systems. 

Perform preliminary design studies of the inertial and 
electro -optical sensor subsystems and indicate areas of 
technology where state-of-the-art advances are 
necessary. 
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• Establish the performance capabilities of the preliminary 
modular design and verify by performance analyses that 
this design meets the performance requirements for each 
mission. 

•A "conceptual design" is a functional representation of the guidance and 
control system component configuration responsive to a given mission, and 
includes 1) a functional schematic blocking out each component subsystem, 
the mechanization of the various operational computations, all data flow, 
and all moding and switching functions, 2) functional descriptions, 
performance characteristics and development status for each component 
subsystem. 

A "preliminary modular design" is a selection of specific components, 
meeting the composite requirements for all the missions considered, that 
may be interchangeably combined into given operational systems for 
specific applications. Such a design includes 1) block schematics of the 
complete complement of guidance and control components selected on the 
basis of the analysis leading to, and the evaluation of, the various con- 
ceptual designs, 2) functional descriptions, physical characteristics, 
performance specifications and interface characteristics for each of the 
modular elements, 3) specification of the mechanical and electrical inter- 
faces between the modular elements of the system and between the system 
and the launch vehicle or spacecraft. 

The initial objective of Task III was to formulate the requirements 
for an integral modular guidance, navigation, and control system capable 
of meeting the mission requirements of Earth low-altitude polar and 
synchronous equatorial orbits, lunar orbit. Mars orbit, and solar probe 
(Jupiter flyby) missions. The results of Tasks I and II (see Ref. 1-1) 
provided the basis for this Task III formulation. Conceptual designs re- 
sponsive to these requirements were then to be developed. Parametric 
variations of the performance characteristics of each of the critical com- 
ponents and subsystems of these conceptual designs were to be analyzed 
so as to permit the establishment of specific performance requirements 
relative to mission accuracy, fuel expenditure, system reliability, and 
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weight. These analyses were to be used under the Task IV effort in speci- 
fying a "Preliminary Modular Design" and in assuring a technically sound 
rationale for the equipment specifications. 

The study constraints and the scope of work applicable to Task III 
can be summarized as follows: 

a) The representative missions to be studied were 

1) Earth-Polar Orbit -Injection Mission utilizing Atlas/ 
SLV3A/Burner II. 

2) ' Synchronous Equatorial Earth-Orbit Mission utilizing 

the Atlas SLV3X/ Centaur. (Both direct ascent and 
parking orbit modes were to be considered. ) 

3) Mars Orbiter Mission (Voyager spacecraft launched 
by Saturn V). 

4) Lunar Orbiter Mission (Lunar orbiter spacecraft 
launched by Atlas SLV3C/ Centaur) . 

5) Solar-Probe Mission using Jupiter assist (advanced 
planetary probe spacecraft launched by Saturn IB/ 

Centaur). (Close-in solar probe (0. 1 AU) and out- 
of-ecliptic missions were to be considered.) 

b) The resultant guidance and control instrumentation for a 
given set of launch vehicles, upper stages, unmanned space- 
craft, and missions was to be based upon the boost phase 
(launch through injection) requirements as well as those for 
midcourse, target approach, encounter, and deboost into 
orbit phases of flight, 

c) - .The choice of inertial systems was to be limited to strapdown 

systems. 

d) Only the existing NASA and DOD radio tracking systems . 
were to be considered. (See Ref. 1-1.) 

..e) .. Specific control system. design, concepts and interfaces with 
existing boost-vehicle control system elements were.to.be 
established for each of the launch vehicles. No attempt was 
to be made to optimize the total control system, design. ' 

f) • -• Onboard -computational requirements (memory capacity,' - 
yjord length, and execution time) w.ere to be. established , 
utilizing the NASA-ERC United Aircraft computer concept 
' described in Ref. 1-2. Sizing studies were to be based 
on. guidance equations preyiously developed by TRW plus 
the control equations developed in this study. 
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g) Computer interfaces were to be defined with respect to 
the NASA-ERC UAC computer concept defined in Ref. 1-2. 
Interface hardware (input/ output) preliminary design was 
to be accomplished, but no specific design information 
for the computer was required. 

h) The planet tracker used in the approach guidance system 
for the Mars mission was to be the NASA-ERC Kollsman 
sensor currently under advanced development. This is 
the only practical approach to the problem of planet 
tracking to have been developed to date. 

1. 2 STUDY IMPLEMENTATION 

In accordance with the objectives stated above, the study effort was 
divided into two groups: 

• Derivation of Guidance and Control Functional and 
Performance Requirements (Vol. II, Part I) 

• Definition of Mission Characteristics 

• Conceptual Design 

• System Performance Analyses 

• Preliminary Modular Design (Vol. II, Part II) 

« System Configuration and Interfaces 

• Subsystem Design Studies 

• Performance Analyses of Modular Design 

Vol. II of this report contains the detailed study results obtained under the 
Task III and IV effort, and is published in two parts as indicated above. 
The following paragraphs describe briefly the implementation of each of 
these groups of tasks. 

1. 2. 1 Derivation of Guidance and Control Functional and Performance 
Requirements (Vol. II, Part I, Secs. 2 Through 5) 

1. 2. 1. 1 Mission Characteristics (Sec. 2) 

Reference trajectories for the five basic missions were developed 
by TRW through the use of its Multivehicle N-Stage (MVNS) and Space 
Navigation Simulation (SNS) precision integration programs (Refs. 1-3 
and 1-4). However, for the earth- synchronous-orbit and lunar -orbit 
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missions, utilizing the Atlas /Centaur launch vehicle, TRW used the' 
trajectories generated under the previous study effort (see Ref. 1-1). New 
trajectories generated under Task III were: 

Reference Powered Trajectories 

a) Atlas Burner II — Low altitude earth-circular polar orbit 
mission launched from WTR. 

b) Saturn V— Launch-to-inj ection trajectory, with 
earth -injection conditions determined to match the 
interplanetary trajectories defined below. 

c) Saturn IB/Centaur — Launch-to-inj ection trajectory 
with earth-injection conditions chosen to match the 
interplanetary trajectories defined below. 

Reference Interplanetary Trajectories 

a) Mars Orbiter Missions — Based on trajectory and 
mission analyses conducted for the Mars 1975 launch 
opportunity under the TRW Voyager Task D study 
(Ref. 1-5), Types I and II reference trajectories 
were' selected for the two Mars orbiter missions to 
be considered. The rationale for selection of the 
reference trajectories is presented together with a 
comparison of the heliocentric trajectory character- 
istics of both mission types in sec. 2 of this report. 

Injection state vectors for these Mars missions were 
computed analytically assuming a Saturn V launch 
vehicle and a 100-n. mi. , short-coast, circular 
parking orbit. 

b) Jupiter Flyby Missions — Trajectory data for Jupiter 
flyby missions during the 1972 launch opportunity 
were generated for the two'specified flyby missions. 
Reference trajectories were selected and an analytic 
computation of the injection state vector was per- 
formed assuming a Saturn IB /Centaur launch vehicle 
and a 100-n. mi. , short-coast, circular parking orbit. 

In addition to generation of the analytic state vector 
required at injection, the vehicle 1 s position with 
respect to the sun, earth, _tar get planet, and Canopus 
was determined for all reference trajectories ana- 
lyzed. Time histories of these quantities were de- 
veloped for both the near-earth and heliocentric 
phases of the missions. Target planet approach 
geometry was defined for all reference trajectories 
and capture conditions and orbit orientation geom- 
etry were developed for the Mars orbit missions. 
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Characteristics of the powered flight and interplanetary trajectories are 
described in sec. 2 of this volume and in subsecs. 2. 3 and 2.5 of Ref. 1-1. 
Characteristics of the Lunar and Mars orbits used for orbit determination 
performance studies are described in subsec. 4.7. 

1. 2. 1. 2 Guidance and Control System Conceptual Designs (Sec. 3) 

The implications of mission objectives on variable versus fixed 
time-of-arrival midcourse guidance schemes were examined for the 
Mars Type I mission, including tradeoffs between midcourse correction 
capabilities and requirements, and for the particular guidance schemes. 

Earth-based tracking and computation was established as the 
primary navigation mode for the lunar and interplanetary missions and 
for the determination of orbital parameters for the Mars orbiter missions. 

The booster and spacecraft attitude control system concepts were 
examined, and a digital system was selected as the most appropriate for 
the applications considered. Control system interface tradeoff studies 
were conducted to define the functional interfaces between the ROI guid- 
ance and control system and the existing or modified boost vehicle control 
electronics and thrust vector and reaction control systems. 

Special emphasis was placed on studies relating to attitude -fixed 
versus attitude-maneuvering spacecraft/payloads and gimbaled versus 
fixed optical sensors. For the translunar and interplanetary coast phases, 
body-fixed optical sensors were selected as the most appropriate space- 
craft attitude references. For the Mars approach guidance, high- 
precision optical measurements are required. Gimbaled Canopus and 
planet sensors were chosen as the most appropriate for this application. 

F-or each mis sion/booster/payload, an overall functional description 
and schematic of the radio/optical/ strapdown inertial guidance system 
were developed. These include the general signal flow, and moding and 
switching functions. Detailed mechanization equations were defined as 
required to define the data flow between subsystems and the operational 
moding and sequencing functions. 
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1. 2. 1.-3 Guidance and Navigation Performance Analyses (Sec. 4) 

. ,a) Sun- Sighting Time -Updating Technique 

The time -updating technique for the multiparking orbit 
earth- synchronous missions was analyzed in detail. 

The. accuracy of this method and the impact on overall 
system accuracy were as.sessed, 

b) Powered Flight Performance Analysis 

For those missions in which the guidance, navigation, 
and control system under study has prime control over 
the boost and injection phase, the GEAP II error 
analysis program (Ref. - 1-6) was used to evaluate injec- 
tion accuracy and to establish the requirements for 
midcourse velocity corrections. Parametric tradeoff 
studies involving strapdown inertial instrument quality, 
and prelaunch azimuth alignment errors, were 
performed. Midcourse correction, deboost maneuver, 
and orbital transfer maneuver accuracies were also 
evaluated. 

c) Interplanetary and Approach Navigation Analysis 
The SVEAD program (Ref. 1-7) for estimating navi- 
gation accuracy was modified to give it the capability 
of handling closed orbits around Mars. The analyses 
made earlier under Task II for the Mars mission 
were extended both to incorporate variations in 
optical sensor accuracies and to examine the impli- 

. cations of Type I versus Type II trajectories. 

1.2. 1.4 Control System Performance Analyses (Sec 5) . 

Bending modes were generated for the Saturn Y/Voyager vehicle 
configuration, and existing bending data, propellant sloshing data, aero- 
dynamic and mass properties data, and thrust vector control characteris- 
tics for each launch vehicle were assembled for use in subsequent control 
system analyses (see apps. C,D,E). 

Stability margins of the linear control system for the first stages of 
the selected boost-vehicle configurations were determined (see app. A). 

A comparison was made between the use of first-stage rate gyros and 
upper-stage gyros, and the digital compensation required under these 
conditions established. Stability margins for the Voyager spacecraft were 
also determined. 

Coast-flight attitude-r.efer.ence acquisition,- maneuvers, and normal 
mode operations were analyzed. 
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1-2.2 Preliminary Modular Design (Vol. II, Part II, Secs, 6 Through 11) 

Preliminary modular designs were developed for each mission based 
on the conceptual designs. Interface definitions were established for the 
onboard computer; the control system; and the telemetry, tracking, and 
command system. Detailed equipment descriptions and specifications 
were developed for the electro -optical sensors and the inertial reference 
unit. 

1.2. 2.1 System Configuration and Interfaces (Sec. 6) 

a) Modularity Concept 

An equipment modularity concept for the total radio/ 
optical/ strapdown inertial guidance system was 
established in accordance with the basic ROI Study 
objectives. For each of the missions, TRW established 
the equipment utilization concept, and defined the 
interconnections and interfaces of the various units 
comprising the system. 

b) Vehicle Interfaces and Mechanical Mounting 
Considerations 

Physical locations and interconnections of the modular 
radio/ optical/ strapdown inertial guidance system 
components were established for each of the five 
launch vehicle /mission combinations, considering the 
optical sensor line -of- sight requirements and other 
location and mounting constraints. Interfaces with 
existing vehicle control system elements were 
established in accordance with the control system 
conceptual and modular design studies. 

c) Guidance Equipment Mechanical Interface and 
Packaging Considerations 

Sensor mounting provisions (necessary for adequate 
mounting stability) were established including the 
requirements for precision navigation base assemblies. 

An electronics packaging modular design concept was 
also established. 

d) Thermal Design Considerations 

For each mission, the expected thermal environment 
conditions and constraints were established for the 
guidance and control equipment at the appropriate 
location in the launch vehicle, upper stage, or space- 
craft. A survey was conducted to establish the 
approximate operating temperature range for the 
most critical optical sensors, and thermal control 
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concepts established to the extent possible using 
available design data on various boost vehicles and 
spacecraft. (See Appendix F) 

Temperature control requirements and concepts were 
established for such units as the IRU, where the 
required performance can be achieved only through 
precise thermal control of critical elements. 

1. 2. 2. 2 Onboard Computational Elements (Sec. 7) 

Onboard computational requirements (memory size, word length, 
and speed requirements) were established for the NASA-ERC UAC 
computer concept (Ref. 1-2), These studies were based on equations 
previously developed by TRW for the LM Abort Guidance System (attitude 
reference navigation computations) (Ref. 1-8), Advanced Centaur Studies 
(steering and guidance computations) (Ref. 1-9) plus the control 
equations developed in this study. 

Major emphasis was placed on defining in detail the I/O interfaces 
between the computer and the electro-optical sensors, the inertial refer- 
ence unit, the control system components, and the telemetry, tracking 
and command subsystems. A conceptual design of a computer interface 
unit (CIU) was developed providing interface compatibility with the 
NASA-ER.C UAC Advanced Kick Stage Guidance Computer (Ref. 1-2). A 
reliability estimate for this computer was developed for use in mission 
reliability studies. 

1. 2. 2. 3 Control System Design (Sec, 8) 

Tradeoffs were made between control system digital autopilot equa- 
tion complexity and computational -time and memory-storage requirements. 
Several digital compensation filters were considered to determine the cost 
of added flexibility in the digital control system. 

The interface between the computer and the control system hardware 
was defined with considerations given to signal levels issued to the thrust- 
vector actuation system and to the receipt of signals from interfacing 
gyro packages. An evaluation was made of the signal mixing requirement, 
either within or outside the computer for differential roll control; the 
problem of interfacing with a varying number of boost-vehicle engines was 
also addressed. 
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Functional schematics of the Atlas /Centaur, Saturn V/Voyager, and 
Saturn IB/Centaur digital control systems were generated showing signal 
flow, and moding and switching functions. 

1.2. 2. 4 Electro-Optical Sensor Designs (Sec. 9) 

The optical sensors selected under Tasks I and II (Ref. 1-1) were 
reviewed both in light of recent state-of-the-art developments and of 
new requirements resulting from present mission specifications. Spe- 
cifically, the applicability of gimbaled Canopus and planet approach 
sensors and the use of a very narrow field sun sensor were considered. 

Based on this review, TRW 'chose a set of sensors appropriate to 
the study effort and established a configuration for each mission. Sensor 
specifications were prepared covering functional description, accuracy, 
physical performance, and reliability. Emphasis was placed both upon 
a preliminary description of data interface characteristics and upon the 
accuracy of the sensor configuration for the Mars orbit mission neces- 
sary for support of the guidance accuracy studies. The state of develop- 
ment of each of the sensor elements was evaluated. 

For each sensor required by the several missions, a preliminary 
design was generated using available data on existing equipment, where 
applicable, plus additional preliminary design effort as required. The 
following characteristics were established for each sensor: 

a) Sensor operating modes 

b) Sensor accuracy 

c) Final data interface characteristics 

d) Weight, dimensions, electrical power requirements 

e) Sensor reliability models and numerical parameters 

f ) Mechanical and electrical mounting characteristics 
consistent with required physical interchangeability 

g) Physical description consisting of a preliminary 
design drawing for each sensor. 
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1.2. 2. 5 Inertial Reference Unit (Sec. 10) 


A preliminary design was generated of a strapdown IRU meeting the 
performance requirements of the several missions based on the previous 
studies carried out under Tasks I and II (see Ref. 1-1). The following 
characteristics were established for this unit: 

a) Sensor and electronics accuracy including environmental 
sensitivities (linear and rotational acceleration 
and vibration environments). 

•b) Data interface characteristics 

c) Weight, dimensions, electrical power requirements 

d) - Mechanical mounting characteristics 

e) Mechanical electrical packaging and thermal con- 
trol concepts 

f) Reliability estimate. 

1 . 2. 2. 6 Performance Characteristics of Modular Design (Sec. 11) 

A performance analysis summary for the preliminary modular 
design was established based on the recommended sensor selections and 
specifications, demonstrating that the preliminary modular design satisfies 
the guidance and control requirements for the five missions studied. The 
overall system performance characteristics were related to trajectory 
accuracy and fuel required for correction of guidance, navigation and con- 
trol errors. 

• Weight, power, and total failure-rate estimates were made for each 
of the elements comprising the modular system and the results used to 
estimate the overall system reliability, weight, and power for each of 
the five missions considered. 

1. 3 DEFINITION OF TERMS 

6 

Certain of the definitions pertaining to the missions, the launch 
vehicle, mission events, and trajectories used throughout this report 
are summarized below. 
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1. 3. 1 Missions 


In general, the term "mission” is used in this report to encompass 
and describe the events which are associated with directing the launch 
vehicle or the spacecraft from the earth and which terminate with the 
accomplishment of the mission objectives. In the analysis of the various 
missions described in the ROI Study, the following terms are used: 


Synchronous Earth In the synchronous earth orbit 

Orbit Mission mission, the launch vehicle is 

used to place the satellite payload 
into an earth-synchronous (24-hr 
period) equatorial orbit at a desired 
longitude. The injected payload 
(satellite) is assumed to have orbit 
trim and stationkeeping capability. 


Orbiter Missions In an orbiter mission, approximately 

at the time when the spacecraft is 
closest to the target body (moor or 
planet), its trajectory is deliberately 
altered by a propulsive maneuver 
so that it remains in an orbit about 
the target body as a satellite. 


Solar Probe Mission 


Flyby Mission 


Solar Probe with 
Planetary Swingby 


In a solar probe mission the space- 
craft is injected into a heliocentric 
orbit that passes within a specified 
distance of the sun. This is an 
untargeted mission requiring no 
trajectory alterations subsequent 
to injection. 

In a flyby mission, the spacecraft 
passes close to the target planet. 

No propulsion forces are employed 
to alter the trajectory so as to remain 
in the vicinity of the target planet. 

The spacecraft departs from the re- 
gion of the target planet, although its 
trajectory will have been perturbed. 

In this type of mi^hion the spacecraft 
passes close to a planet with the 
purpose of significantly altering 
the spacecraft trajectory. After 
departure from the target planet, 
the spacecraft continues on a helio- 
centric trajectory to within a pre- 
scribed distance from the sun. No 
propulsive forces are employed to 


12 



Solar Probe with 
Planetary Swingby 
(Continued) 


1. 3; 2 Vehicle Terms 


Launch Vehicle 


Kick Stage 


High Energy Upper 
Stage (HEUS) 


. Spacecraft 


alter the trajectory in the vicinity 
of the target planet. For a given 
distance of closest approach to the 
sun, this technique may be used 
to significantly reduce the launch 
vehicle AV requirements, usually 
at the expense of considerably 
longer mission durations. 


The launch vehicle includes the 
multistage boost vehicle which 
injects the spacecraft into the 
desired trajectory and includes 
all hardware up to the interface 
where the spacecraft is mated and 
the payload shroud attaches which 
protects the spacecraft. Generic- 
ally, the launch vehicle system 
also includes all appropriate 
ground support and' test equipment. 

For the purposes of this study, 

"kick stage" refers to the final 
powered stage of the launch vehicle 
(the payload spacecraft is assumed 
to have only limited velocity capa- 
bility for incremental orbit correc- 
tions). The kick stage is assumed 
to provide complete three-axis 
guidance, navigation and control 
capability for all launch vehicle 
stages except- for the Saturn V 
(Mars orbiter mission). 

This is a particular kick stage con- 
cept using an advanced propulsion 
system burning high-energy propel- 
lants such as H 2 /F 2 - Typical gross 
weight is 3200 kg. The thrust to 
weight ratio is approximately 1. 

The spacecraft system encompasses 
the payload itself and all its compo- 
nent subsystems, the science pay- 
load, the adapter which is mounted 
to the kick stage, and limited propul- 
sion capability for orbital corrections. 
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Launch Operations 
System. 


Mission Operations 
Systems 


1.3.3 Mission Events 


The launch operations system does 
not include any flight hardware, but 
constitutes the operational responsi- 
bility for supporting and conducting 
the launch of the combined launch 
vehicle and spacecraft through the 
separation of the spacecraft from 
the launch vehicle. 

Operational responsibility for sup- 
porting and conducting the mission 
after the spacecraft is separated 
from the launch vehicle is borne 
by the mission operations system. 


In the analysis .of the various mission events described in the ROI 


Study, the following terms are used: 
Prelaunch 

Launch 

Liftoff and Ascent 


Collectively, all events before 
liftoff. 

Collectively, all events from 
liftoff to injection. 

Departure of the combined launch 
vehicle -spacecraft from the 
ground and ascent to a parking 
orbit of specified altitude (typi- 
cally 185 km (100 n. mi). 


Injection (synchronous Thrust termination of the kick stage, 

earth orbit mission) placing the kick stage /payload into 

a transfer trajectory to synchronous 
altitude from the parking orbit or, 
alternately, into the final syn- 
chronous earth orbit. 


Injection (lunar or 
interplanetary mission) 


Separation (shroud) 


Separation (spacecraft) 


Thrust termination of the lower 
stages of the launch vehicle, plac- 
ing the kick stage /payload into an 
interplanetary or translunar trajec- 
tory, from the parking orbit. 

Detachment of the nose fairing from 
the launch vehicle during ascent. 

Detachment of the spacecraft from 
the spacecraft kick stage adapter 
after injection. 
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Orientation Maneuver 


A programmed alteration of the 
injection stage or spacecraft attitude 
to cause it to return to a desired 
orientation such as the cruise 
orientation. 


Midcourse Trajectory 
Correction Maneuver 


A propulsive maneuver performed to 
compensate for inaccuracies or 
perturbations so as to redirect the 
spacecraft toward the intended aim- 
ing point. Generally, it requires 
orientation to a specific attitude, 
operation of the rocket engine, and 
reorientation to the cruise attitude. 
The time of this maneuver is dur- 
ing the interplanetary or translunar 
flight, but not necessarily at the 
midpoint. 


Encounter 


Orbit Insertion 


Generally, encounter encompasses 
events occurring when the spacecraft 
is near the target planet. Specifically, 
it refers to the time when the space- 
craft is at its point of closest 
approach (periapsis). 

The propulsive braking maneuver by 
which the (orbiter) spacecraft tra- 
jectory at the target planet is changed 
from approach (hyperbolic) to orbital 
(elliptical). 


i.3.4 Trajectory Terms 

In discussing the trajectories possible for the various missions of 
the ROI Study, the following terms are used: 


Direct Trajectory An interplanetary trajectory from 

the earth to a target planet, in which 
no intermediate planets (or satellites) 
are approached closely enough to 
significantly influence the trajectory. 

Swingby Trajectory An interplanetary trajectory from the 

earth to a target planet, in which an 
intermediate planet is passed suffi- 
ciently closely to exploit the effect 
of its gravitational attraction. This 
exploitation may provide reduced 
mission duration, reduced launch 
energy, or an opportunity for scien- 
tific observations of the intermediate 
planet. 
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Launch Opportunity 


Launch Period 


Launch Window 


Geocentric 

(heliocentric, 

planetocentric) 


C3, Launch Energy, 
Injection Energy 


Asymptote 


DLA 


ZAL 


ZAP 


The time during which trajectories 
to a target planet may be initiated 
from the earth, with reasonable 
launch energies. A launch oppor- 
tunity is usually identified by the 
year in which it occurs, and the 
target planet. 

The space in arrival date-launch 
date coordinates in which earth- 
planet trajectories are possible in 
a given launch opportunity; speci- 
fically, the number of days from 
the earliest possible launch date 
to the latest. 

The time in hours during which a 
launch is possible on a particular 
day. 

Described or measured with respect 
to inertial coordinates centered with 
the earth (sun, planet). Pertaining 
to the portion of the flight in which 
the trajectory is dominated by the 
gravitation of the earth (sun, planet). 


Twice the geocentriq energy-per- 
unit mass, of the injected space- 
craft. This is equivalent to the 
square of the geocentric asymp- 
totic departure velocity. 

The line that is the limiting position 
which the tangent to a hyperbolic 
(escape) trajectory approaches at 
large distances from the attracting 
center. 

Declination of the outgoing geocentric 
launch asymptote. 


Angle between the outgoing geocen- 
tric asymptote and the sun-earth 
vector. 

Angle between the incoming planeto- 
centric asymptote (at the target 
planet) and the planet-sun vector. 
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ZAE 


or V 


HP 


Parking Orbit 


Type I, Type II 

Interplanetary 

Trajectories 


Angle between the incoming planeto- 
centric asymptote (at the target 
planet) and the planet-earth vector. 

Planeto centric asymptotic approach 
velocity. 

An unpowered, geocentric, approxi- 
mately circular orbit, separating 
the powered portions of the launch 
and injection sequence. 

Type I transfers are defined as 
those in which the vehicle traces a 
central angle of less than 180 ° 
about the Sun between departure 
from the Earth and arrival at the 
planet. In Type II transfers, the 
angle is greater than 180 °. 


± 4 3^ 5 Coordinate Systems 

The various coordinate systems used in specifying performance re- 
quirements and powered flight performance analysis results obtained 
during the HOI Study are defined as follows: 

EC I (Earth-Centered- This is a right-handed coordinate 

Inertial) system, in which Z lies along the 

earth’s polar axis and X and Y lie in 
the earth's equatorial plane. The X- 
axis passes through the Greenwich 
meridian or in the direction of the 
Vernal Equinox at the time of launch, 
(specified in text). 


RTN (Radial-Tangential- 
Normal) 


A right-handed orthogonal coordinate 
system in which R lies in the direc- 
tion of the nominal position vector 
from the center of the earth, and 
N lies in the direction of the orbital 
angular momentum. T forms a 
right-handed orthogonal set with R 
and N. 


(X, Y, Z) Selenographic Moon-Centered Inertial Coordinates. 

This is a right-handed orthogonal 
coordinate system in which Z lies 
along lunar polar axis, and X, Y lie 
in the lunar equatorial plane with X 
passing through zero lunar longitude 
(Sinus Medii). 
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B • R, B • T 


For a given interplanetary trajec- 
tory, the impact parameter vector 
B specifies in which direction from 
the planet and what distances the 
approach asymptote lies. B is 
commonly expressed in components 
B ■ R and B * T, where R, S, T are 
a right-hand set of mutually orthogo- 
nal unit vectors aligned as follows: 

S is parallel to the planet centered 
approach asymptote, T is parallel to 
the plane of the ecliptic and positive 
eastward, and R completes the set 
and has a positive southerly 
component. 
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2. MISSION CHARACTERISTICS 


2. 1 INTRODUCTION 

The five basic missions to be investigated under Task III of this 
study determined the specific launch vehicle/payload combinations 
described in this section. Table 2-1 summarizes the mission- related 
data pertaining to the launch vehicles and the location of the Radio/ 
Optical/Strapdown Inertial Guidance System (ROI) for each of these mis- 
sions. The table also delineates the specific TRW assumptions made 
relative to the guidance regime. 

Those specific mission characteristics added to Task III or different 
from those used in the Tasks I and II portions of the overall study are 
examined in detail in the following subsections. However, the earth- 
synchronous satellite and lunar orbiter missions are substantially the 
same as those used in Tasks I and II; descriptions of their' characteristics 
are contained in subsecs 2. 3 and 2. 5 of Ref. 2-1. 

2. 2 EARTH LOW-ALTITUDE POLAR-ORBIT MISSION 

The earth low-altitude polar -orbit mission typifies one that might be 
used for earth resources studies and was added to the repertory of the 
basic Tasks I and II missions for more complete coverage of the spectrum 
of possible unmanned space missions. This study assumed that the orbiting 
satellite payload is capable of correcting for orbit-insertion errors. 
Typically, this AV capability can be on the order of 10 to 20 m/sec. The 
modular guidance system must then provide the guidance function from 
launch through orbit insertion with accuracy sufficient to ensure that the 
payload AY capability is not exceeded. 

For this mission, the Atlas/Burner II launch vehicle is assumed to 
be launched from the Western Test Range (WTR). The Atlas stages inject 
the Burner II/payload combination into a coast up to the apogee altitude 
of 927 km. At that altitude the Burner II provides the velocity increment 
for circularizing the payload orbit. The actual sequence of events is 
summarized in Table 2 -II. Basic data used to define the launch and injec- 
tion trajectory and this sequence of events were obtained from Refs. 2-2 
through 2-5. 
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TABLE 2-1 


RADIO/ OPTICAL/ ST RAPD OWN INERTIAL TASK III MISSION SUMMARIES 





Guidance 

Approximate 


Mission 



Package 

Payload 


Trajectory Characteristics 

Booster 

Location 

Weight 

Guidance Regime 

Earth Low -Altitude 

WTR Launch; ~927 km near polar 

- Atlas SLV-3A/ 

Burner II 

2, 500 

Launch through insertion of payload into 

Polar Orbit 

circular orbit 

Burner II 



desired earth orbit (payload assumes 
orbit trim and stationkeeping functions 

Earth-Synchronous 

Same as used in Tasks X and II 

Atlas SLV-3C/ 

Centaur 

400 

after separation from last booster stage) 

Orbit 

(See Ref. 2-1) 

C entaur 




a) Direct Accent 

b) Parking Orbit Accent 






Lunar Orbitcr 

Same as used in Tasks I and II 

Atlas SLV-3X/ 

Payload 

2 , 000 

Launch through insertion into final 



Centaur 


desired sclenocentric or areocentric 
orbit, including all midcourse cor- 


Mars Orbitoi, 1975 

Saturn V injects spacecraft with 

Saturn V 

Payload 

40 . 000 

rections and orbit change maneuvers 

a) Typo 1 Traioctory 

desired C 3 into Typo I or Type II 
interplanetary trajectory; space- 




b) Type 11 Trajectory 

craft performs M/C and deboost 
into 1100 x 10 , 000 -km orbit and 
subsequent injection into 500 -km 
orbit 





Jupiter Flyby 

S-IB/Contaur injects spacecraft 

Saturn IB/ 

Payload 

800 

Launch through injection into inter- 

a) 0.1— AU Probe 

onto a high-energy interplanetary 
trajectory (C 3 = iEl km^/sec^, 

C entaur 


planetary orbit and pro-oncountcr 
midcourse correction s) ; post- 


b) Cross Ecliptic Probe 

T i 464 days) post Jupiter tra- 
jectory determined by targeted 
B • T. B • E 




encounter attitude control only 



TABLE 2 -II 


ATLAS SLY -3 A/ BURNER II SEQUENCE OF EVENTS 


Event 

Description 

Time 

(sec from Liftoff) 

T 

LO 

Liftoff 

0 . 0 

BECO 

Booster engine cutoff 
(sustainer operation) 

148. 4 

JBP 

Jettison booster package and 
s hr oud 

151.4 

SECO 

Sustainer engine cutoff 

361.4 

VECO 

Vernier engine cutoff (begin coast 
to apogee) 

381. 1 

B21G 

Burner II ignition 

1107. 1 

B2BO 

Burner II burnout 
(circular orbit injection) 

1153. 1 


The characteristics of the actual orbit obtained from the TRW/ 
N-Stage program include the following: 


Injected weight 
Inclinati on 
Apogee /perigee 
Eccentricity- 
Orbital period 


2513 lb 
99° 

954/ 900 km 
0. 0037 
103. 54 min 


This is not a perfectly circular orbit. Since the above orbit was adequate' 
for error analysis purposes, further iterations of the N-Stage program to 
achieve a more circular orbit were not attempted. 
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2. 3 MARS ORBITER MISSIONS 


Two 1975 Mars orbifcer missions, corresponding to Type I and Type 
II transfers, have been selected from the optimum 20-day launch periods 
identified during the Voyager Task D study (Ref. 2-6). These two types of 
trajectories were chosen to examine the sensitivity of the trajectory de- 
termination errors (and hence the fuel required for corrective maneuvers) 
to guidance and control errors. For either' type of trajectory, the basic 
mission phases listed below are identical: 

a) Launch, parking orbit, and injection into interplanetary 
traj ectory 

b) Separation from booster and first-cruise phase 

c) Midcourse execution 

d) Subsequent cruise and midcourse corrections 

e) Approach 

f) Deboost velocity application (into 1, 100 x 10, 000-km orbit) 

g) Doppler tracking in elliptic orbit 

h) Transfer into 500-km altitude circular orbit 

Within each launch period, the critical mission was identified as that 
Earth-Mars trajectory which requires the maximum short coast Earth 
parking orbit. Table 2-III summarizes the Saturn V launch vehicle char- 
acteristics used to compute these coast times. The basic booster data 
was obtained from Refs. 2-3 and 2-7. Specific launch sequence event 
times for the Type II transfer are summarized in Table 2-IV. 

Table 2-V lists the pertinent trajectory characteristics of each 
critical mission; Figures 2-1 and 2-2 illustrate the heliocentric transfer 


Type I transfers are defined as those in which the vehicle traces a 
central angle of less than 180° about the sun between departure from the 
Earth and arrival at Mars. In T : ype II trajectories, the angle is greater 
than 180°. The two types are effectively noncontiguous: when the helio- 
centric central angle is very near 180°, the position of Mars out of the 
ecliptic causes the interplanetary trajectory to be highly inclined to the 
ecliptic, leading to excessive launch energy requirements. 
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geometry of each mission. Time histories of the following trajectory char- 
acteristics are displayed in Figures 2-3 and 2-4 for the transit phase of 
each mission: 


a) Sun- spacecraft distance . 

b) Sun -Mars distance 

c) Spacecraft-Earth distance 

d) Earth-Mars distance 

e) Spacecraft-Mars distance 

f) Sun-spacecraft-Earth angle 

g) Sun-Mars -Earth angle 

In addition, the Sun-Mars distance, Earth-Mars distance, and Sun-Mars - 
Earth angle plots have been extended to include the first 200 days of the- 
orbiting phase of each mission. 

TABLE 2-IH 


SATURN V LAUNCH AND INJECTION 
TRAJECTORY CHARACTERISTICS 


1975 MARS TRANSFERS 


Phas e 

Dui'ation 

(min) 

Angle 

(deg) 

Altitude 
(n. mi.) 

Type I Mars Transfer 




Total powered flight 

17.28 

49. 80* 


Circular parking orbit 

58.02 

236.84"' 

100.0 

.Inj ection 


8.0 

180.0 

Type II Mars Transfer 




Total powered flight 

17.20 

49.20* 


Circular parking orbit 

24.83 

101. 34"' 

100.0 

Inj ection 


8.0 ** 

180.0 


Angle traversed, measured in earth-centered inertial coordinates. 
Flight path angle at injection, measured (+) above the local horizontal. 
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TABLE 2-IV 


SATURN V SEQUENCE OF EVENTS 
(TYPE II TRANSFER TRAJECTORY) 


Event 

Description 

Time 

(sec from liftoff) 

T 

LO 

Liftoff 

0.0 

IECO 

S-IC inboard engine cutoff 

154.6 

OECO 

S-IC outboard engine cutoff 

158.6 

T 

IG2 

S-II stage ignition 

164. 1 

T 

JFI 

Jettison S-IC /S-II forward interstage 

194. 1 

T 

x JHS 

Jettison heat shroud 

214. 1 

S2CO 

S-II stage cutoff 

538. 1 

T 

IG3 

S-IVB stage first ignition 

543.6 

T 

x B03 

S-IVB stage first cutoff 
(parking orbit injection) 

686. 2 

1 2IG3 

S-IVB stage second ignition 

2184.9 

S4CO 

S-IVB stage final cutoff 
(transfer orbit injection) 

2491.3 


The distances plotted in Figures 2-3 and 2-4 affect communications 
characteristics (spacecraft- earth distance) and relate to solar radiation 
and wind intensities ( spacecraft- sun distance). The sun-spacecraft-earth 
angle is significant because of its effect on the transfer of attitude refer- 
ence from earth to sun for the performance o£ midcourse maneuvers. 

2.4 SOLAR PROBE (WITH JUPITER ASSIST) MISSIONS 

It has been shown (Ref. 2-8) that the gravitational field of Jupiter 
may be employed to obtain solar probe and out-of- ecliptic postencounter 
trajectories following a close flyby past that planet. A 1972 solar impact 
mission and a 1972 90° out- of- ecliptic mission have been analyzed assum- 
ing the Saturn IB /Centaur launch vehicle characteristics given in Table 
2 -VI (Refs. 2-3, 2-9, and 2-10). 
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TABLE 2-V 


CHARACTERISTICS OF 1975 EARTH-MARS TRAJECTORIES 



Type I 
Transfer 

Type II 
Transfer 

Departure date 

1975 September 19 

1975 September 22 

Arrival date 

1976 May 1 

1976 September 5 

Time of flight, days 

224. 75 

348.32 

Departure asymptote 
{from earth) 



km/sec 

4.45 

3.85 

C^, km^/sec^ 

19.76 

. 14.83 

Angle to equatorial 
plane, deg . 

50'. 12 

5. 13 

Angle to sun- earth 
line, deg 

248.94 

255. 14 

Heliocentric Orbit 

- 

% 

True anomaly at 
departure, deg 

1. 565 

: 0.899 

True anomaly at 
arrival, deg 

7. 204 

-8.558 

Heliocentric transfer 
angle, deg 

150.68 

203.32 

Inclination to ecliptic, deg 

3.751 

- 2.083 

Perihelion distance 
from sun, AU 

1. 003 

1.003 

Aphelion distance 
from sun, AU . 

1. 705 . 

1.675 

Eccentricity 

0. 2594 

0.2510 

Approach Asymptote (to Mars) 



km/sec 

3. 09 

2.80 

Angle to plane of Mars' 
orbit, deg 

-20.22 

26.83 

Angle to Mars -Sun line, deg 

138. 76 

54.71 
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Figure 2-4. Time Histories of Heliocentric Orientation Angles and Distances for 1975 Type XI Mars Mission 


ANGLES (DEG) 






TABLE 2-VI 


SATURN IB /CENTAUR LAUNCH AND INJECTION 
TRAJECTORY CHARACTERISTICS 


Phase 

Duration 

(min) 

Angle 

(deg) 

Altitude 
(n. mi.) 

Total powered 
flight 

18.12 

. 53.70 


Circular parking 
, orbit 

1.08 

4.42"' 

100.0 

Injection 


siesj: 

12.7 

313.0 


Angle traversed, measured in earth- centered inertial 
coordinates. 


'r # 

Flight path angle at injection, measured (+) above the 
local horizontal. 

Specific launch sequence event times for the Saturn IB/ Centaur are 
summarized in Table 2 -VII. 


TABLE 2 -VII 

SATURN IB/ CENTAUR SEQUENCE OF EVENTS 


Event 

Description 

Time 

(sec from liftoff) 

T 

LO 

Liftoff 

0.0 

IECO 

S-IB inboard engine cutoff 

139.5 

OECO 

S-IB outboard engine cutoff 

146. 0 

T 

IG2 

S-IVB ignition 

151.5 

JHS 

Jettison heat shr.oud 

181. 5 

4BCO 

S-IVB cutoff 

643. 0 

Cl 

Centaur first ignition 

654. 5 

CECO ’ 

Centaur first shutdown 
(parking orbit injection) 

660 . 4 

CI 2 

Centaur second ignition 

725.4 

ceco 2 

Centaur second cutoff 

1152.4 
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Following the Centaur second cutoff, the payload coasts in the helio- 
centric transfer ellipse to Jupiter encounter. The earth-centered and 
heliocentric transfer trajectory characteristics of both missions are 
essentially the same. The altitude of closest approach at Jupiter and the 
components B • T and B ■ R of the impact parameter B determine the 
postencounter trajectories. The impact parameter, B, is defined as a 
vector originating at the center of the target and is perpendicular to the 
incoming asymptote, V {see Figure 2-5). A unit vector T is defined as 
lying in a plane parallel to the ecliptic according to 


T = 


Y x k 
oo 

I^CO * 


where k is a unit vector normal to the ecliptic plane and pointing towards 
the north. The R axis is defined by 


R = 


V 


co 



x T 
x T I 


The impact parameter, B, lies in the R-T plane and has components 
B * T and B • R. 
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Table 2 -VIII summarizes the pertinent trajectory characteristics of 
each mission. The velocity of the solar probe as it becomes enveloped in 
the sun's photosphere is 617.45 km/sec; the total flight time beginning 
from injection is 2. 762 yr. 

TABLE 2 -VIII 

CHARACTERISTICS OF 1972 JUPITER PROBES 
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3. GUIDANCE AND CONTROL SYSTEM CONCEPTUAL DESIGNS 


3. 1 INTRODUCTION 

This section presents the recommended guidance and control system 
conceptual designs for each of the five missions considered in this study. 

A conceptual design is defined as a functional representation of the com- 
ponent configuration responsive to a specific mission, and consists of the 
following: 

1) A functional schematic of the complete guidance, 
navigation, and control system indicating all in- 
formational loops. 

2) Performance descriptions of each component and 
component subsystem in terms of its functional 
description, accuracy, physical parameters, and 
reliability. 

3) Statement of development status of each component. 

Functional schematics for each of the missions are presented in subsec. 

3. 3 following. Component descriptions, performance characteristics, and 
development status supporting the conceptual designs are presented in 
secs. 7, 9, and 10. 

The guidance and control conceptual designs presented in this sec- 
tion are based on the operational sequences and the guidance performance 
requirements developed under Tasks I and II (Ref. 3-1). These require- 
ments have been refined and extended to reflect the revised mission de- 
finitions and the five specific launch vehicle/ payload combinations defined 
in subsec. 1. 1. Subsec. 3.2 summarizes the guidance and control require- 
ments on which the conceptual designs are based. Subsec. 3.3 presents 
the conceptual designs for each mission. 

The guidance system core concept adopted during the Tasks I and 
II studies was retained in this study. However, some of the basic func- 
tional concepts have been modified. In particular - , the utilization of the 


' The characteristics of the missions that differ from those that were used 
in the Task I and II studies (Ref. 3-1) are described in sec. 2 of this 
volume. A summary of the characteristics of all the missions is given in 
sec. 3 of vol. I. 
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inertial measurement unit and. digital computer of the core configuration 
was extended to cover the launch and boost phases of all the missions. 

This modification of concept was made to examine the feasibility of per- 
forming launch/boost/ injection guidance and control with a strapdown 
inertial package. In most instances, it is difficult to justify (on the basis 
of cost, performance requirements, and payload weight and size) the need, 
or use, for a complete three-axis inertial measurement unit to be used 
solely for attitude control and midcourse velocity corrections in inter- 
planetary missions. The addition of the launch and boost-phase guidance 
and control functions to the total set of functions to be performed by the 
system thus provides a tenable basis for including the three-axis inertial 
measurement unit for these missions. 

The recommended conceptual guidance system and control configu- 
rations developed in this study for the boost vehicles considered herein 
ignores the basic fact that all these boosters already have highly developed 
or proven guidance packages of their own. However, it was not intended 
to propose replacement of the existing systems with the strapdown system 
of this study. Rather, the boosters used in this study served primarily as 
vehicles or bases from which the analytical and preliminary design studies 
could proceed. 

3. 2 GUIDANCE AND CONTROL REQUIREMENTS 

This subsection summarizes the guidance and control requirements 
which form the basis for the conceptual designs presented in subsec. 3. 3. 
Guidance system performance (accuracy) requirements are summarized in 
vol. I, sec. 3. 

The. control subsystem conceptual designs presented in this section 
are based on the functional requirements defined in this subsection, the 
performance (stability) requirements defined in app. A, and the control 
system interface considerations and tradeoffs discussed in sec. 8. 

3. 2. 1 Guidance and Control System Operational Sequences 

The functional requirements for the ROI Guidance and Control Sys- 
tem are developed from the operational sequences necessary to accomplish 
the objectives of the various missions with the postulated launch vehicles. 
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The operational sequences are presented in sec. 3 of vol. I for each of 
the five missions considered. A brief summary of these sequences is 
given here. 

For the lunar and interplanetary orbiting missions, the general 
operational sequence may be summarized as follows: 

a) Launch and Boost into Parking Orbit. The ROI strap- 
down inertial guidance and control subsystem will 
provide the guidance and control functions for this 
phase for all the missions. 

b) Coast in Parking Orbit . The spacecraft, together 
with injection stage, will coast in the parking orbit 
until translunar or interplanetary orbit injection. 

The inertial guidance subsystem serves as an 
attitude reference during this phase. 

c) Translunar or Interplanetary Orbit Injection. The 
injection stage is ignited to inject the spacecraft 
into the translunar or interplanetary trajectory. 

Attitude and burn control will be provided by the 
strapdown inertial guidance subsystem. 

d) Coast Until First Midcourse Correction. Following 
injection burn, the spacecraft is separated from the 
injection, stage, a celestial reference acquisition sequence 
is initiated, and the spacecraft becomes attitude 

fixed to the sun and star Canopus via body -fixe d sun 
and star sensors. The strapdown accelerometers 
can be turned off (except for heaters), and the flight, 
computer algorithm for updating the direction cosines 
can be placed in a standby mode. 

Deep-Space Network (DSN) tracking is used during 
this coast phase to determine the orbit and compute 
the midcourse velocity correction required to reduce 
the effects of injection errors. The midcourse 
thrust vector pointing and magnitude commands and 
time of execution command are transmitted to the 
ROI guidance system. 

e) First Midcourse Correction . At a predetermined time 
from injection, the first midcourse correction is 
executed. Ten to thirty minutes prior to the time of 
execution, the strapdown inertial reference unit is 
turned on, the direction cosine solution algorithm is 
initialized, and vehicle rotations are commanded to 
orient the thrust vector in the required inertial direc- 
tion. The midcourse burn is initiated at the correct 
time when the proper attitude is achieved, and is de- 
signed to null selected miss components at lunar or 
planetary intercept. 
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f) Second Coast Phase and Second Midcourse Correction, 

After completion of the first midcourse correction, the 
spacecraft is "unwound" to the original Sun / Canopus 
reference attitude and continues in a cruise phase iden- 
tical to the first. A second midcourse burn may be 
performed prior to lunar or planetary injection if re- 
quired. 

g) Coast Until Deboost Maneuver Into Intermediate Lunar 
or Planetary Orbit. This phase is identical to the other 
coast phases. 

h) Deboost Into Intermediate Orbit. Based on the DSN tracking 
data obtained during the previous phase, the deboost maneu- 
ver is calculated and the spacecraft is deboosted into an 
interinediate orbit, under control of the inertial guidance 
system. 

1) Coast in Intermediate Orbit. The spacecraft is tracked by 
DSN stations to determine orbital parameters and tne retro - 
maneuver required to place the payload into the final orbit. 

j) Retro Into Final Orbit. Based on the orbital estimates 
obtained from DSN tracking data and controlled by the 
strapdown inertial guidance system, the spacecraft is 
injected into the final orbit. 

For the Solar Probe (with Jupiter assist) missions, the operational 
sequence is similar except that no further powered maneuvers are required 
subsequent to the last midcourse corrections. 

The operational sequence for the earth orbiting missions consists 
basically of: 

a) Launch and boost into parking orbit 

b) Coast in parking orbit 

c) Injection into transfer orbit 

d) Transfer orbit coast 

e) Injection into final orbit 

Functionally, the operations in each of these phases is the same as 
for the corresponding phases in the lunar and interplanetary missions except 
that earth-based tracking is not utilized. 
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For the low-altitude polar-orbit mission utilizing the Atlas/Burner 
II, a direct ascent mode is utilized (no parking orbit is required). 

3. 2. 2 Control System Functional Requirements 

The conceptual guidance and control system designs defined in subsec. 
3. 3 are based on the requirement for stabilizing and controlling each stage 
of the boost vehicle and also any spacecraft stage that employs the system 
during both powered and coast phases. Booster operation is assumed to 
include all mission phases between liftoff and separation of the spacecraft. 

In the event that the last stage of the boost vehicle is used as a spacecraft, 
the spacecraft functional requirements would be applicable to that stage. 

In addition to the functional requirements defined below, the control 
system designs must meet the specific performance (stability) requirements 
defined in app. A. No attempt has been made to optimize interfaces be- 
tween the ROI guidance and control system and the existing boost vehicle 
control systems components. Existing boost vehicle control electronics 
and -actuators are used where possible, however, some modifications are 
su gg es ted to simplify the overall control system and to achieve interface 
compatibility. Interface tradeoffs are discussed in sec. 8 of this volume. 

3. 2. 2. 1- Boost Vehicle Functional Requirements 

Multistage launch vehicles often employ coast phases between the 
upper -stage powered flight phases, particularly if third or fourth stages 
are used. Requirements for control during these coast phases are included 
within the booster stage control requirements. 

Powered Flight. The ROI Guidance and Control System must maintain 
vehicle control throughout powered booster flight from liftoff to 
separation of the spacecraft. Adequate stability margins, as described 
in app. A, must be maintained throughout this phase. Shortly after lift- 
off, the system must initiate a pitch program to achieve the desired 
vehicle trajectory. Closed-loop guidance steering during first- stage 
operations is optional and dependent upon mission requirements. 

Vehicle control must be maintained in the presence of disturbances 
such as those due to winds, nose -fairing jettison, and stage separation. 

A control system autopilot design intended to relieve structural 

loads resulting from high-wind shears is required only for those vehicles 
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presently incorporating a load relief design. Of the candidate vehicle 
configurations, only the Saturn IB booster employs such a system. 

The control system autopilot must have a minimum capability of two 
changes in the gain and filter coefficients per stage, one change during 
each stage operation and one change at the start of each stage operation, 
excluding a change at liftoff. Increases in these requirements are to be 
made to satisfy the stability requirements given in app. A. 

In addition to control loop gains and stability compensation filters, 
the autopilot must contain attitude -command error limits and thrust-vector 
deflection-command limits which are capable of being changed at the start 
of each stage operation. The autopilot control laws must be capable of 
accomplishing vehicle control through linear thrust-vector deflections as 
well as pulsed reaction jet control. Linear thrust-vector control techniques 
include thrust deflection by means of liquid injection, gimbaled engine, 
and gimbaled nozzle control. However, for the boost vehicles under con- 
sideration here, only gimbaled engine control is required. The pulsed 
jet control laws must be applicable for coast phase as well as powered 
phase operation through changes made in autopilot coefficients. 

Coast Flight Operations. During the coast phases of boost flight, 
the control system must maintain three-axis attitude control of the booster, 
through the use of reaction control jets. It must also be capable of accom- 
plishing attitude change maneuvers to attain the desired orientations for 
subsequent powered flight operations. The coast phase limit cycle per- 
formance prior to spacecraft separation must be within acceptable staging 
requirements for the particular launch vehicle design. Such requirements 
will be strongly dependent upon the means adapted for spacecraft separation. 

In addition to attitude control, the coast phase control system must 
be capable of accomplishing forward translational acceleration for propel- 
lant settling operations to facilitate main engine operations and for adding 
low forward velocity increments. 

3. 2. 2. 2 Spacecraft Functional Requirements 

Spacecraft functional requirements are given below with both digital 
and analog control modes included. The analog control modes serve as 
either backup modes selectable upon ground command or as primary modes 
during mission phases when the digital computer is deactivated. 
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Launch Vehicle Separation. The tip off rates imparted to the space- 
craft during separation from the launch vehicle are removed by the ROI 
Guidance and Control System through the use of the spacecraft reaction 
control devices. 

Digital Control Inertial Hold . Subsequent to separation from the 
launch vehicle, the digital control system must provide three-axi3 attitude 
control of the spacecraft employing signals from the inertial sensors in 
an attitude hold mode. 

During all phases of the mission in which the digital computer is 
active, the digital control inertial hold mode must be capable of selection 
through ground commands or onboard sequencing commands. 

Acquisition of Celestial References. The ROI Guidance and Control 
System must reorient the spacecraft to the desired orientation to enable 
acquisition of the celestial references (i. e. , earth, sun, stars, planets) 
through onboard preprogrammed command angles or through ground com- 
manded angles. In the event that one of these references is not directly 
obtainable, a search mode must be provided in which the vehicle is either 
librated or slowly spun about the search axis. 

Acquisition of the celestial references through both digital and analog 
control must be provided by the ROI Guidance and Control System. The 
analog control would be a backup system selectable upon ground command. 
Acquisition of celestial sensors would be accomplished or completed 
through use of the electro -optical sensors. 

Analog Control Inertial Hold . During all phases of the mission in 
which the digital computer is inactive, an analog control inertial hold 
mode must be capable of selection through ground commands or onboard 
sequencing commands. This inertial hold mode must provide three -axis 
attitude control of the spacecraft through use of signals from inertial 
instruments. 

Digital and Analog Celestial Hold . After acquisition of the celestial 
references and through use of the electro -optical sensors, the ROI Guid- 
ance and Control System must maintain a celestial hold mode through both 
digital and analog control selectable upon ground command. During earth 
and planetary orbits, if one of these bodies were employed as a reference. 
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the celestial hold mode would entail one spacecraft rotation per orbit- to 
maintain a specific control axis aligned to the planet geocenter. 

Reorientation Maneuvers. In addition to maneuvers for celestial 
reference acquisitions, the ROI Guidance and Control System must accom- 
plish additional spacecraft reorientation maneuvers through digital or back- 
up analog control employing inertial instruments. The additional maneu- 
vers include reorientation for: powered flight firings; acquisition of solar 
power; improvement of communication strength; experiment pointing or 
tr acking. 

Thermal Control Mode . As a means of reducing the thermal gradi- 
ent between the sunlit and shadowed portions of the spacecraft, the ROI 
Guidance System must be capable of maintaining a low- spin rate on the 
spacecraft. This mode is to be employed over long-mission durations in 
which the digital computer is deactivated, hence, analog control must be 
employed with electro -optical sensing for spin control. 

Spacecraft Powered Flight Operation . During the phases when 
trajectory and orbit changes are accomplished, the ROI Guidance and Con- 
trol System must maintain three-axis attitude control through the use of 
spacecraft reaction control devices or through deflection of the_main engine 
thrust. The use of the inertial sensors would be made for these firings. 
These firing phases include orbit transfer, orbit circularization, orbit- 
trim, translunar and interplanetary orbit injection, midcourse correction, 
and lunar and planetary deboost firings. 

The powered flight operation may include an attitude hold, or con- 
stant angular rate modes, as well as guidance commanded operation em- 
ploying explicit guidance laws. 

3. 3 GUIDANCE AND CONTROL SYSTEM CONCEPTUAL DESIGNS 

The composite conceptual equipment configuration depicted in Figure 
3-1 was developed on the basis of the above premises and the performance 
analyses in secs. 4 and 5 and Ref. 3-1. The utilization of equipment for 
each mission is summarized in Table 3-1. 

A digital control system concept was chosen for these studies be- 
cause it provides a flexible means of implementing the control functions 
on a.per flight basis without requiring hardware modifications. In this 
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Figure 3-1. Composite Equipment Configuration 
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study, full advantage was taken of the digital computer existing within 
the ROI guidance system to accomplish the stabilization and control of the 
boost vehicle as well as of the spacecraft. With the use of the computer 
a single autopilot can be employed to control all the booster and spacecraft 
stages, thereby eliminating the need for multiple autopilots, which are 
often used in multistage space boosters. Other benefits of this approach 
include: 1) the elimination and/ or simplification of certain items of booster 
control system hardware and 2) the ease with which the ROI guidance and 
control system can be adapted to the various booster/ spacecraft TVC and 
RGS systems with a minimum of special-purpose interface hardware. 

The recommended control system conceptual design configuration is 
shown in Figure 3-2. Control system interface considerations and trade- 
offs are discussed in sec. 8. Specific adaptations of this concept for each 
of the missions is given below. 

3.4 NEAR-EARTH POLAR- ORBIT MISSION 

The powered and coast phases of the near-earth polar-orbit mission 
up to injection of the payload into the design orbit is of short duration 
(19.2 min) with no inordinate demands exceeding state-of-the-art guidance 
capabilities. Electro-optical sensors are not required for any mission 
phase; therefore, guidance system for this mission is comprised of the 
core package of Table 3-1. 

The integrated guidance and control configuration is indicated in 
Figure 3-3. The basic guidance package is installed in the Burner II and 
provides the guidance function for the Atlas stages as well. A control 
electronics package is required on the Burner II to interface between 
1) the primary ROI computer and the Burner II attitude control system 
and 2) the ROI computer and the Atlas components, indicated in Figure 
3-3, which are part of the existing Atlas system. 

Autopilot stability studies indicate that the Atlas rate gyros should 
be retained, with considerations of possible relocation (see sec. 5). 
However, the Atlas position gyro functions can be taken over by the ROI 
core package. These comments pertaining to the Atlas hold for the two 
missions discussed in subsecs. 3.5 and 3. 6. 
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•WITHOUT POSITION GYROS. 


Figure 3-3 


Basic Conceptual Design Configuration for the 
Near-Earth Polar-Orbit Mission 
















The guidance performance analysis of this conceptual design can be 
found in subsec. 4. 2. 

3. 5 EARTH-SYNCHRONOUS SATELLITE MISSION 

The integrated guidance and control conceptual configuration for the 
earth- synchro nous satellite mission is indicated in Figure 3-4. An earth 
horizon scanner and a solar aspect sensor (see subsecs. 9.4 and 9. 5 
respectively) have been added to the core package. The core package, the 
electro -optical sensors, and an interface electronics package are installed 
on the Centaur. No changes are made to the basic Centaur control actua- 
tion system, and the Atlas control system configuration is the same as in 

• , 'f 

the previous mission. 

The functioning of the various sensors can best be described with 
reference to the basic mission profile. During the Atlas and first Centaur 
burns to parking orbit injection, guidance and steering are controlled 
inertially. For the direct-ascent mission, the second Centaur burn is 
initiated at first equatorial crossing, approximately half an hour after 
launch. During the intermediate coasting period, constant attitude is 
maintained and the second Centaur burn for Hohmann transfer from parking 
orbit altitude to synchronous altitude is again controlled inertially. For 
this direct-ascent mission, no external attitude or timing update informa- 
tion is required (Ref. 3-1). 

However, for the long parking orbit coast case, both an attitude and 
timing update are highly beneficial prior to the second Centaur, or perigee, 
burn (see Table 4-V) . Both these updates can be obtained with the com- 
bination of the earth sensor and solar aspect sensor shown in Figure 3-4. 

During the long (approximately 5. 25-hr) Hohmann transfer coast to 
apogee at synchronous altitude, attitude is maintained inertially. How- 
ever, prior to the third Centaur, or apogee, burn, an attitude update is 
accomplished again with the aid of the earth and sun sensors. 

The performance achievable with this system configuration is dis- 
cussed in subsec. 4.3. 


The control system design is discussed further in subsec. 


8 . 2 . 
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so ' 



Figure 3-4 









3 . 6 LUNAR ORBITER MISSION 

For the Lunar orbiter mission, the ROI guidance package is installed 
.in the orbiter spacecraft, and data and signal transfer to the Atlas control 
system configuration is effected through a Centaur electronics interface 
package (Figure 3-5). The canopus tracker and sun sensor (see subsecs. 
9.4 and 9. 6 respectively) replace the earth sensor and solar aspect sen- 
sor of the previous mission. These sensors are used to establish the 
celestial attitude reference only during the translunar coasting' phases. 

Performance analyses for this mission were conducted during the 
Tasks I and II phases of this overall study. A summary of the translunar 
orbit injection analysis for the two strapdown inertial systems studied and 
for the Centaur gimbaled inertial guidance system is given in vol. 1, 
subsec. 5. 4. 

3. 7 MARS ORBITER MISSION 

The major difference in the system elements for the Mars orbiter 
mission as compared to those of the lunar orbiter mission is the possible 
addition of the planetary approach sensor. As shown in subsec. 4. 6, data 
from this sensor in conjunction with data from the sun and Canopus sensors 
can be utilized by ground-based stations to improve the quality of the de- 
termination of the spacecraft approach orbit to Mars. However, for mis- 
sion requirements comparable to those in use up to now, it is not clear 
that this improvement in approach orbit determination is absolutely essen- 
tial. Thus, the planetary approach sensor shown in Figure 3-6 is included 
conditionally so that the implications on preliminary modular design can 
be investigated for applications to possible future missions with high- 
• accuracy requirements. 

Except for the planetary approach sensor, the functions and utiliza- 
tion of the total guidance and control system substantially parallel the 
functional operations of the lunar orbiter mission. As shown in sec. 5, 
the Saturn V rate gyros are retained to simplify the autopilot design prob- 
lem. 
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Figure 3-5. Basic Conceptual Design Configuration for the 
Lunar Orbiter Mission 
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Basic Conceptual Design Configuration for the 
Mars Orbiter Mission 


















3. 8 SOLAR PROBE MISSIONS (WITH JUPITER ASSIST) 

Up to Jupiter encounter, the solar probe missions closely resemble 
the lunar mission. Therefore, the conceptual configuration. Figure 3-7, 
is very similar to that of the lunar mission, insofar as the guidance sys- 
tem is concerned. 
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4. GUIDANCE AND NAVIGATION PERFORMANCE ANALYSES 


4. i. INTRODUCTION AND SUMMARY 

In secs. 7 and 8 of the Tasks I and II report {Ref. 4-1) , performance 
analyses for the following miss ions/ miss ion segments were obtained: 

a) Injection of a payload into the earth- synchronous 
satellite orbit using Atlas/ Centaur 

b) Translunar orbit injection with the Atlas/ Centaur 

c) Midcourse and planetary orbit insertion maneuver s 

Under the present study effort, the following performance analyses were 
conducted: 

a) Injection accuracy and 95%-corrective AV require- 
ments for the Atlas/ Burner II near-earth polar 
mission (see subsec. 4.2) 

b) Analysis of the sunsighting, perigee burn-time update 
technique for the synchronous satellite mission (see 
subsec. 4. 3) 

c) Injection accuracy, 95% first midcourse AV, and 
.target miss analysis for the Saturn V/Mars orbiter 
mission (see subsec. 4.4) 

d) injection accuracy, 95% first midcourse AV, and 
target miss analysis for the Saturn- IB/ Centaur/ 

Jupiter missions (see subsec. 4. 5) 

e) Mars approach analysis (see subsec. 4. 6) 

The inertial instrument error models for these analyses are the 
same as those derived in sec. 4 of Ref. 4-1, and are summarized in 
Table 4-1. Both the TG-166 and TG-266 inertial system error models 
were used in the analyses. Either system : is capable of providing adequate 
performance; the TG-166 is the less accurate of the two systems, but is 
available’ at lower cost, weight, size and power requirements. For the 
Earth polar and- Mars missions, the error models were modified to the 
extent that initial azimuth alignment, error source No. 7, was varied 
parametrically to establish requirements 'on azimuth alignment. • The 
error sources attributable to the electro -optical sensors were modified to 
reflect those developed in the present study. 
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TABLE 4-1 


ERROR MODELS USED FOR STRAPDOWN INERTIAL GUIDANCE 

PERFORMANCE ANALYSIS 


Number 

TG- 

166 

TG- 

266 

Units 

Type 

Description 

1 

3 

3 

m 

Initial 

Vertical position 

2, 3 

15 

15 

m 

Initial 

East, north position 

7, 8, 9" 

20 

20 

arc sec 

Initial 

Orientation 

10 

600 

600 

arc sec 

Opt ical 

Roll axis at apogee 

20,30 

180 

180 

arc sec 

Optical 

Yaw, pitch axes at apogee 

11,31 

720 

720 

arc sec 

Optical 

Roll, pitch axes at per igee 

21 

180 

180 

arc sec 

Optical 

Yaw axis at perigee 

39*" 

4. 0^ 

4. of 

sec 

Optical (Sun 
Sensor) 

Update time 

40, 51,62 

21 

14 

gg 

Accelerometer 

Bias 

73,77,81 

‘ 75 

24 

^g/g 

Accelerometer 

Scale factor 

74 

12 

10 

arc sec 

* Accelerometer 

X accelerometer input axis 
rotation toward y-axis 

75 ■ 

■ 12 

10 

arc sec 

Accelerometer 

X accelerometer input axis 
rotation toward z-axis 

78 

12 

10 

arc sec 

Accelerometer 

Y accelerometer input axis 
rotation toward z-axis 

82,83,84 

-15 

10 

^g/g • 

Accelerometer 

Pendulous axis g sensitivity 

85,86,87 ‘ 

1 

1 

4g/g 

Accelerometer 

Output axis g sensitivity 

91, 97, 103 

50 

30 

8-g/g 2 

Accelerometer 

Input -pendulous g-product 
sensitivity 

92,98, 104 

■ 0.5 

0.5 

t*g/ g 2 ’ 

Accelerometer 

Input- output g-product 
sensitivity 

230, 241, 252 

0. 187 

0. 09 

deg / hr 

Gyro 

Bias drift 

263, 266,269 

0. 627 

0. 16 

deg/ hr/ g ’ 

Gyro 

Input axis g-sensitive drift 

264,267,270 

0. 627 

. 0. 16 

deg/hr/ g 

Gyro 

Spin axis g-sensitive drift 

265,268,271 

0. 02 

.0 . 

deg/ hr/ g 

Gyro 

. Output axis g-sensitive. drift 

275, 281,287 

0. 04 

0. 04 

deg/ hr/ g 2 

Gyro 

Anisoelastic drift 

290,294,298 

57 

26 

ppm 

Gyro 

Scale factor. 

291,292,293 

10 

.10 

arc sec 

Gyro 

Gyro input axis rotations 
toward each of other two axes 

295,296,297 

10 - 

10 

arc sec 

Gyro 

Gyro input axis rotations 
toward each of other two axes 


Initial azimuth alignment error (source No. 7) was varied parametrically between 10 
and 180'arc sec. 

Applies only to the multi-parking orbit-synchronous satellite missions. 

+0. 238 sec used in previous error analyses (Ref. 4-1). 
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The error analyses were conducted with the aid of the GEAP II and 
SVEAD error analysis programs, as were Tasks I and II. The basic 
powered trajectories are those described in sec. 2. 

4. 2 POWERED PERFORMANCE ANALYSES OF THE NEAR-EARTH 
POLAR- ORBIT MISSION 

For the near-earth polar-orbit mission, the basic guidance system 
is located in the Burner II. The guidance system configuration contains 
only the core package, comprised of the strapdown IMU package and the 
computer. No optical sensors are required; thus, error source Nos. 10, 
11, 20, 21, 30, 31 and 39 of Table 4-1 do not apply. 

For an initial azimuth alignment error of 20 arc sec, the resulting 
1-crrss position and velocity component errors in radial, tangential, 
normal (RTN) coordinates for the two inertial systems are as indicated in 
Table 4 -II. 


TABLE 4 -II 

ATLAS/BURNER II NEAR-EARTH POLAR-ORBIT INJECTION ERRORS 



Position 

(km) 

Y eloc ity 
(m/ sec) 

System 

R 

T 

N 

R 

T 

N 

TG-166 

1.28 

1. 59 

4.01 

3. 28 

1.50 

6.05 

TG-266 

0. 54 

0. 96 

1.32 

1. 56 

0. 86 

1. 79 


The major contributors to these totals are summarized in Tables 4-III and 
4-IV for the TG-166 and TG-266 systems, respectively. 

To relate these injection errors to mission performance, the 95% 
corrective AV required to correct the payload orbit was computed for both 
systems and for a range of initial azimuth accuracies {error source No. 7). 
The orbit correction scheme was assumed to be the same as that used in 
Ref. .4- 1, i. e. , 

1) Perform a Hohmann transfer from the perigee of the 
imperfect orbit to the desired 926. 5-km (500 n. mi.) 
circular orbit altitude. 
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TABLE 4- III 

SUMMARY OF ATLAS/BURNER II NEAR-EARTH POLAR-ORBIT MISSION ERROR ANALYSIS RESULTS' 

{TG-166 System with Initial Azimuth Alignment Error = 20 arc sec) 



Position 

(km) 


O' 

•0. 70 
0. 16 
0 

-0.09 

0 

0 

1.31 

0 

-0.07 

0.31 

-0.44 


-0. 81 
-0.07 
-0.31 
0 
0 

0 - 
- 0 . 57 
0 

1.02 
-3. 76 
0 
0 


Velocity 

(m/sec) 

Orientation 
(arc sec) 

R 

T 

N 

Yaw 

Roll 

Pitch 

0 

0 

1 

O 

-J 

19. 8 

-22. 6 

0 ' 

0. 73 

-0. 58 

0 

-2. 5 

-2. 6 

19. 7 

O 

O* 

i 

0. 14 

0 

CO 

l 

t 

-12. 9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1. 17 

0 

0 

0 

0 

0 

0 

0 

-1.46 

192. 8 

38. 0 

-8. 7 

-2. 82 

1.28 

0 

-9. 1 

-3. 3 

-215.0 

0 

0 

1.02 

-38. 2 

193. 0 

0 

o 

I 

0 

-5. 73 

518. 0 

180. 5 

24. 6 

-0. 52 

0. 30 

0 

0 

0 

-20. 6 

0. 55 

-0.-40 

0 

0 

0 

' 6. 3 












cr* * 
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TABLE 4- IV 

SUMMARY OF ATLAS/BURNER II NEAR-EARTH POLAR-ORBIT MISSION ERROR ANALYSIS RESULTS 

(TG-266 System with Initial Azimuth Alignment Error = 20 arc sec) 



Position 

(km) 


Velocity 
(m/ sec ) 


Or ientation 
(arc sec) 


R 

' T 

, N 

R 

T 

N 

Yaw 

Roll . 

Pitch 

7 

• -0.05 

0 

-0. 81 

0 

0 

1 

o 

-a 

19. 8 

-22. 6 

0 

8 

0. 29 

O 

R- 

• 

O 

! 

r- 

o 

o 

1 

0. 73 

-0.58 

0 

0 

0 

19. 7 

9 

-0.09 

0. 16 

-0. 31 

-0. 20 

0. 14 

0 

-14. 8 

-12. 9 

0 

40 

0. 06 

0 

0 

0 

0 

0 

0 

0 . 

0 

62 

0. 03 

-0.06 

0 

0 

0 

0 

0 

0 

0 

73 

0. 23 

0 

0 

0.37 

0 

0 

0 

0 

0 • 

230 

-0. 02 

0 

-0. 27 

0 

0 

0. 70 

93. 0 

18. 0 

0 

.241 

-0.44 

0. 63 

0 

-1. 37 

0. 61 

0 

0 

0 

-104.0 

252 

.0.03 

0 

0.49 

0 

0 

0.49 

-18.0 

93. 0 

0 

263 

-0. 05 

-0. 02 

-0. 96 

0 

0 

1.46 

132. 2 

46.0 

-6. 3 

294 

0. 12 

-0. 20 

0 

0. 24 

0 

0 

0 

0 

0 



2) Circularize when the desired 926. 5-km altitude is 
reached. 

3) Perforin the inclination correct at equatorial crossing. 

The results from a 1000-run Monte Carlo analysis of the above correction 
sequence are indicated in Figure 4-1. Both the average AV and 95% AY 
requirements are shown. As expected, the TG-266 system shows a per- 
formance effectiveness two to three times better than the TG-166. The 
"knee" of the 95% AV curves occurs near 20 arc sec, however, no really 
significant degradation occurs before i arc min. Operationally, a value 
of 20 arc sec can be achieved by optical means and is a recommended 
prelaunch value. 



Figure 4-1. Corrective Velocity Requirements for Atlas /Burner II 
Near-Earth Polar-Orbit Mission 
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4. 3 ANALYSIS OF THE SUN SIGHTING, PERIGEE BURN-TIME UPDATE 

In the Task II report {Ref. 4-1) , it was recommended that a time 
update be provided to ensure that Centaur ignition for the transfer from 
the 185 -km parking orbit to synchronous altitude be executed at equatorial 
crossing. The benefits of this concept are evident, comparing the orbit 
trim 95% AV requirements shown in Table 4-V, which were extracted 
from the referenced report. 

It was proposed that the time update be achieved by means of a 
sun/ earth sighting technique. The pertinent geometry for this scheme is 
illustrated in Figure 4-2. At equatorial crossing, the angle E measured 
in the orbital plane from the ascending node is 0 or it, depending on 
whether ignition is desired at the ascending or descending node, respec- 
tively. Then the zenith angle at equatorial crossing is determined from 

cos A = ± (cos Q cos t - sin 12 sin t cos i ) ■ (4-1) 

O 0 

wher e 

Cl = longitude of the ascending node 

T = sun angle measured in the ecliptic plane 
from vernal equinox 

i = obliquity of the ecliptic = 23. 4° 

+ is used for ignition at the ascending node 
is used for ignition at the descending node 

Thus, with known nominal values for Q and t, A q is pre computable. 
By continually comparing this value against the zenith angle. A, to the sun 
measured in orbit via the earth scanner /digital aspect sensor combination, 
ignition can be commanded when A q - A = 0. The primary constraint on 
this proposed scheme is that for a period of 5 to 15 min prior to, and up 
to, equatorial crossing, 45 < A < 90°. The timeupdate accuracy of this 
technique depends on sensor accuracies, the accuracy of the preknowledge 
of launch time (i. e. , error in t), and the accuracy of injection into, the 
185 -km parking orbit (i.e., errors in fi and i , the orbital inclination, 

q s 

nominally 28 ). 
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TABLE 

SUMMARY OF SYNCHRONOUS MISSION 


u> 


Run 

No. 

Coast 

Orbits 

System 

No. 

Pre launch 
Calibration 

Time 

Update 

1 

0 

166 

No 

No 

2 

0 

166 

No 

No 

3 

0 

166 

Yes 

No 

4 

0 

166 

Yes 

No 

5 

0 

266 

No 

No 

6 

0 

266 

No 

No 

7 

8 

166 

No 

No 

8 

8 

166 

No 

Yes 

9 

8 

166 

Yes 

No 

10 

8 

166 

Yes 

Yes 

1 1 

8 

266 

No 

No 







12 


8 


266 


No 


Yes 


Attitude Update 

95% AV (m/sec) 
Required by Payload 
for Orbit Trim 

Perigee 

Apogee 

No 

No 

73 

No 

Yes 

13 

No 

No 

75 

No 

Yes 

9 

No 

No 

35 

No 

Yes 

8 

Yes 

Yes 

163 

Yes 

Yes 

22 

Yes 

Yes 

109 

Yes 

Yes 

19 

Yes 

Yes 

83 

Yes ' 

Yes 

19 
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PLANE OF THE 
ECLIPTIC 

EQUATORIAL 
PLANE 

ASCENDING NODE 

CONSTRAINTS: 

45 £ A s 90° AT 
PERIGEE BURN AND FOR 
AT LEAST 5 MIN PRIOR 

cos A = (cos n cos r - sin a sin r cos I ) cos E 

e 

- (sin fl cos t + cos Q sin t cos i ) cos i + sin t srn t sin i sin E 
l e s e s j 

Figure 4-2. Sun Sighting Time Update Technique for Multi-Parking 
Orbit Synchronous Satellite Mission 

4. 3. 1 Visibility Constraints 

The 90° limit on A is imposed to preclude the possibility of having 
to look through the earth's atmosphere and the attendant refraction effects. 
The lower limit is imposed to ensure a reasonable angular spread between 
the lines -of- sight to the sun and earth for accurate vehicle attitude updat- 
ing (with A = 0, the earth and sun are 180° apart and no vehicle yaw 
attitude information can be obtained). 

Equation (4-1) can be manipulated to determine those combinations 
of Q and t that satisfy the 45° S A < 90° constraint. Some typical results 
are shown in Figure 4-3. For illustrative purposes, consider, the curve 
labeled (50°, A)., Any Q , t combination falling on this curve gives .an 
angle A of 50° at the ascending node and further,, 45 s A < 90° for at 
least 5 min prior to reaching the ascending node. 

4. 3. 2 Accuracy Considerations 

In general, satisfaction of A - A = 0 results in an error in E' from 

o 

the desired 0 or tt value because of the errors mentioned previously. 

The possible errors in t should be considered first. 
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Figure 4-3. Combinations of Q and t that Satisfy 

the Visibility Constraints 

If there were no error in launch time or orbital period, and if T 
were known precisely at launch, it would be easy to compute the value of 
t that would exist at any equatorial crossing. Errors in launch time or 
in orbital period would cause a deviation from this nominal value. But t 
varies nearly linearly with time, so that if t is known at any time, an 
accurate value of t can always be obtained from the computer by peri- 
odically updating T under control of the internal clock. If this scheme 
is used, the only errors in t arise because of clock errors, and launch 
time and orbital period errors have no effect on r . 

The sun angle varies by only 360°/yr, or approximately l°/day. 

If an update were performed every minute, the maximum error would be 
about 0. 04 arc min. In practice, updates could be performed much more 
frequently if desired. The contribution of the clock error should be com- 
pletely negligible. 
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If the method for updating t described above is used, errors in t 
become uncorrelated from the injection errors. Also, errors in the sen- 
sor measurement are uncorrelated with all other errors. Errors in Q 
are a function of injection errors, and it will be shown that errors in the 
orbital inclination, i , have no effect. Thus, the total error in angle E 
can be written simply as 


6E 



(4-2) 


The various sensitivities may be determined by differentiating the 

complete expression for cos A shown in Figure 4-2 with respect to all 

the variables, and evaluating these derivatives at the nominal conditions 

(i.e., at E = 0° ) . When this is done SE/di is seen to be zero. The 

s 

other sensitivities are 


3E 

SQ 


sin Q cos t + cos Q sin t sin i 

e 


D 


K, 


(4-3) 


9E 

9t 


cos Q sin T + sin Q cos t cos i 
e 

D ' ' 



(4-4) 


3E _ sin A _ ^ 
9A D 3 

where 


(4-5) 


D = (sin Q cos t + cos Q sin t cos i ) cos i + sin t sin i sin i 

e s e s 

(4-6) 

Then, the variance of E is given by 


2 

°E 





(4-7) 


Since errors in measured A. and t are uncorrelated with each other and 
with anything else, 
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and 


< 7 = 8 T 

T 


ct a= 5A 

From the covariance matrix of RTN position and velocity error components 
at injection into the 185 -km parking orbit, <r^ may be computed. 

The angular momentum vector of a satellite is perpendicular to the 
orbit plane and is given by H = R x V, where R and V are the position 
and velocity vectors of the satellite. The cross product of a unit vector 
along the Z-axis (jx^.) with the unit vector along the H vector (jx^) results 
in a vector in a direction perpendicular to both of them (jx^. To be 
perpendicular to H, it must lie in the orbit plane, and therefore passes 
through the ascending and descending nodes. Because of the order in 
which the cross products were taken, its positive direction is toward the 
ascending node. But the angle Q is defined as the angle between the 
X-axis and the ascending node, so that 

cos Q = jx^ • jx^ (4-8) 


But 


1 

'2 


8 Q = - 


1 - { % * ii Q ) 


d(R x * R Q ) 


d( %' d( ^Q ) 


since jx x is a constant, therefore 


(4-9) 

(4-10) 


d( ttn> = 


d 


(jx z X jx H ) 
|i±Z X %l 


(4-11) 


After the differentiation has been performed, the partial derivatives 
may be written as 
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an 

ax 

an 

ay 

an 

az 

an 

ax 

an 

ax 

an 

az 


-z T1 

B 

Z T2 
B ' 

X T1 - Y T2 
B 

Z T1 
B 

-Z T2 
B 

-X T1 + Y T2 
B 


(4-12) 


where 

T1 = YZ - ZY 

T2 = XZ - ZX 
2 2 

B = T1 - T2 

• • • 

and where X, Y, Z, X, Y, Z are the components of position and velocity. 
Letting F denote the matrix of partial derivatives shown above, then 



. -(4-13) 


If C denotes the covariance of the state vector, then 


n 


S F C F‘ 


(4.-14) 


The resulting value for Oq is 5'. 2 arc min. 

It was assumed that r was updated as infrequently as every 1'5 'sec. 
Then the maximum error 8 t is approximately ‘0. 0! arc* min, ' 'Tile efror 
in measured zenith angle attributable to the solar aspect and earth’ sensors 


- - 1 
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was taken to be conservatively 20 arc min (the solar aspect sensor accuracy 
is approximately 3 arc min and the earth horizon scanner accuracy is 
approximately 12 arc min). 

Evaluation of the sensitivity coefficients, Equations (4-3) through 
(4-5), shows these to be all of approximate unity magnitude, and thus, cr_ 
is approximately 20.6 arc min. The timing error is related by 

ct T = °e /e 

or, for the 185 -km parking orbit altitude, 

QVj, = 4 sec 

In the Task II study, a value of = 0.238 sec was used; thus, the new 
value is a factor of 17 larger. 

However, even this magnitude increase does not degrade ultimate 
performance as shown in Tables 4-VI and 4-VII. These tables summarize 
the major error contributors to Runs 8 and 12, respectively, of Table 4-V. 
The underlined error source, No. 39, is the time update error and the 
values correspond to those for et„ = 0. 238 sec. Increasing these by a 
factor of 17 results in the rss errors labeled "new" in Tables 4-VI and 
4-VII. The resulting 95% A V requirements have increased by only 
1 m/sec in each case. Thus, even with the relatively coarse time update 
of 4 sec, the total performance is still significantly improved over the 
cases without time update (see Runs 7, 9, and 11 of Table 4-V). 

4. 4 POWERED PERFORMANCE ANALYSES OF THE MARS MISSIONS 

4. 4. 1 Injection Errors 

Using the basic error source values of Table 4-1 and the powered 
trajectory to injection into the heliocentric transfer orbit to Mars, the 
resulting 1-a rss position and velocity component errors in the RTN 
coordinates for the two inertial systems are as indicated in Table 4-VIII. 

The major contributors to these totals are the same as for all previous 
powered error analysis runs, viz, initial alignment, x- and y-accelerom- 
eter biases, x-accelerometer scale factor, gyro bias drift, and y-gyro 
output axis g- sensitive drift and scale factor. 
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TABLE 4 -VI 


SUMMARY OF LONG COAST SYNCHRONOUS SATELLITE MISSION 
ERROR ANALYSIS RESULTS (TG-166 SYSTEM) 



Initial orientation 
Roll attitude at apogee burn 
Roll attitude at perigree burn 
Yaw attitude at apogee burn 
Pitch attitude at apogee burn 
Pitch attitude at perigee burn 
Time Update 

X- | i 

> Accelerometer bias < 

Y- ) ( 

X- Accelerometer scale factor 

X- j f 

Y - > Bias drift j 

Z- ) ' 

Y- Gyro OA g-sensitive drift 
Y- Gyro scale factor 


Error 

Source 


RSS (old) 
RSS (new) 


Position 

(km) 


59- 1 1 28.0 

59. 5 148.2 


Velocity 

(m/sec) 

Orientation 
(arc sec) 

R 

T 

N 

Yaw 

Roll 

Pitch 

-2. 6 

0. 6 

0. 3 

0 

0 

0 

3. 8 

-0. 1 

0 . 0 

14 

475 

257 

-0. 5 

-0. 6 

1.4 

0 

0 

0 

0. 0 

-0. 5 

-1.0 

120 

- 3 

- 2 

-0. 6 

0. 0 

0. 0, 

0 

57 

-106 

-5. 9 

0. 3 

0. 2 

0 

0 

0 

-0. 2 

0.0 

0.0 

0 

0 

0 

0. 8 

-0. 1 

-0. 1 

0 

0 

0 

-0. 2 

0. 0 

0. 0 

0 

0 

0 

4. 2 

-0. 6 

-0. 4 

0 

0 

0 

0. 8 

-0. 1 

-0. 1 

2 

113 

47 

4. 3 

-1.4 

-0. 8 

0 

47 

-113 

0. 1 

0. 5 

0. 8 

-1 28 

2 

1 

1. 0 

-0. 3 

-0.2 

0 

-2 

- 3 

-1.7 

0.4 

0. 2 

0 

0 

0 

9.9 

2.0 

2. 1 




10. 4 

2. 0 

2. 1 

















TABLE 4 -VII 


SUMMARY OF LONG COAST SYNCHRONOUS SATELLITE MISSION 
ERROR ANALYSIS RESULTS (TG-266 SYSTEM) 




Initial orientation 

9 

Roll attitude at apogee burn 

10 

Roll attitude at perigree burn 

11 

Yaw attitude at apogee burn 

20 

Pitch attitude at apogee burn 

30 

Pitch attitude at perigee burn 

31 

Time Update 

39 

X- Accelerometer bias 

40 

X- Accelerometer scale factor 
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230 


Bias drift 


241 

Z- 



252 

Y- Gyro scale factor 

294 



RSS (old) 



RSS (new) 


Position 

(km) 

Velocity 

(m/scc) 

Orientation 
(arc sec) 

95% AV 

R 

T 

N 

R 

T 

N 

Yaw 

Roll 

Pitch 

(m/sec) 

-19. 2 

1 8. 3 

9.9 

-2. 6 

0.6 

0. 3 

0 

0 

0 


0. 1 

0. 0 

0. 0 

3. 8 

-0. 1 

0. 0 

14 

475 

257 


- 3.7 

2. 1 

1.9 

-0. 5 

-0. 6 

1.4 

0 

0 

0 


0. 0 

0. 0 

0. 0 

0. 0 

-0.5 

-1.0 

120 

-3 

- 2 


0. 0 

0.0 

0. 0 

-0.6 

0. 0 

0. 0 

0 

57 

-106 


-12. 1 

117.6 

63.7 

-5.9 

0. 3 

0. 2 

0 

0 

0 


0. 3 

4. 4 

2.4 

-0. 2 

0. 0 

0. 0 

0 

0 

0 


3. 6 

-4.7 

-2. 5 

0.5 

-0. 1 

-0. 1 

0 

0 

0 


8.7 

-11. 2 

-6. 1 

1.3 

-0. 2 

-0. 2 

0 

0 

0 


- 0.2 

- 1.0 

2.0 

0.4 

0. 0 

0. 1 

1 

54 

23 


21. 3 

-10. 6 

-5.7 

2. 1 

-0. 7 

-0.4 

0 

23 

- 54 


0.2 

1. 9 

-4. 0 

0.0 

0. 3 

0. 4 

-62 

1 

0 


- 6.1 

4.7 

2.5 

-0. 8 

0. 2 

0. 1 

0 

0 

0 


33. 2 

121.0 

66.0 

8.0 

1. 3 

1.8 




19 

33. 6 

142. 2 

77, 6 

8.7 

1. 3 

1. 8 




20 












TABLE 4-VHI 


SATURN V MARS MISSION INJECTION ERRORS (INITIAL 
AZIMUTH ALIGNMENT ERROR = 20 ARC SEC) 


System 

Position 

(km) 

R T N 

Velocity 
(m / sec) 

R T N 

TG-166 

TG-266 

7.1 2.01 3.37 

3.37 1.02 1.07 

7.96 2. 89 4.84 

3. 52 1. 30 1.61 


As the initial azimuth alignment value (error source No. 7) is varied 
parametrically, only the normal components of position and velocity in 
Table 4-VIII vary. The variation is as indicated in Figure 4-4. The 
20-arc sec value again appears to be a satisfactory compromise between 
operational feasibility and system performance. 

4.4.2 Target Misses Due to Injection Errors 

The sensitivity matrix, C, relating injection errors to target misses 
in the commonly accepted B * T and B * R coordinates is 

9 (B * T ) 5 (B ' T) a (B ’ T r ) 

dX dY SZ 

a (B • r ) 

ai J 

The numerical values of the sensitivity coefficients are obtained from the 
"Interplanetary Search Program" (ISP), and for the Type I and Type II 
transfer trajectories they are approximately 


C = 


a (B * T ) 5 (B T) 5 (B T) 

ax 5Y a Z 

a (B • R ) 
ax 


0. 1188E5 

-0. 336E5 

-0. 121E5 

0.455E8 

-0. 260E7 

0. 241E8 

-0. 155E5 

_ 

0. 315E5 

0. 156E5 

-0. 429E8 

-0. 153E7 

-0. 284E8 

-0. 962E5 

0. 164E6 

0. 943E5 

-0.224E9 

-0. 156E9 

0.685E8 

0.267E5 

L 

-0.424E5 

-0. 274E5 

0.610E8 

0. 387E8 

0. 183E8 


where 

the E denotes powers of 10, and the units are km/km 
and km/km/sec. 
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NORMAL POSITION (KM) AND VELOCITY ERRORS (M, 




Evaluation of Eq. (4-15) yields, for the Type I trajectory, 


S 


T 


7. 8954E10 
-7. 2562E10 


-7. 2562E10 
6. 6956E10 


km 2 


(4-16) 


and for the Type II trajectory 


r 


T 


1. 3754E12 
-3.7192E11 


-3.7192E11 
1. 0058E11 


km 


(4-17) 


The resulting error ellipses corresponding to the covariance matrices of 
Eqs. (4-16) and (4-17) are depicted in Figures 4-5 and 4-6, respectively. 

4. 4. 3 Corrective Velocity Requirements 

To determine the midcourse velocity required to correct for the 
injection errors, one of the options of GEAP II is used. This option 
requires the inputting of the 3x6 guidance sensitivity matrix G, which is 
composed of the partials of -the three midcourse velocity components 
(AV^, AV^, anc ^ 2 ) with respect to the six injection error components. 
The matrix G is obtained from ISP for specific times of application of the 
midcourse correction. Two types of G can he obtained: one for miss plus 
time-of-arrival correction and one for miss only. The program outputs 
are the midcourse correction velocity covariance matrix, and the 

expected 95% AV value. The results obtained from this GEAP II option 
are summarized in Table 4-IX. 

4.4,4 Corrected Target Miss 

Inasmuch as the midcourse velocity corrections will not be applied 
perfectly, the target miss ellipses of par. 4.4, 2 will not be eliminated 
completely. From subsec. 8.4 of Ref. 4-1, the velocity application error 
along the thrust and cross axis directions is 
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B • R 



I I I I i I I i I ! I 

10 864202468 10 (KMX10 5 ) 


Figure 4-5. Mars Trajectory, Type I (Uncorrected Miss Ellipse) 


B • R 



Figure 4-6. Mars Trajectory, Type II (Uncorrected Miss Ellipse) 
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TABLE 4- IX. 


NINETY-FIVE PERCENT AV MIDCOURSE (5 DAYS) 
REQUIREMENTS FOR 1975 MARS MISSIONS 


Trajectory 

Type 

Correction 

Type 

System 

95% AV 
(m/ sec) 

I 

M + T"* 

TG- 166 

77.8 

I 

*'i'*T* 

M 

TG- 166 

70. 5 

II 

M + T 

TG- 166 

77. 1 

II 

M 

TG-166 

57. 6 

I 

M + T 

TG-266 

35.4 

ft 

I 

M 

TG-266 

32. 0 

II 

M + T 

TG-266 

35. 7 

II 

M 

TG-266 

26. 6 


Miss plus time-of-flight correction. 
Miss correction only. 


e (AV) 
AV 


0.43 x 10“ 3 



e (AV . ) 

cross axis 

AV 



and the execution error covariance matrix is 


S 

e 




) 2 


AV 




I 


where 



expected value of V squared 


I = unit matrix 







The target miss covariance matrix following the imperfect midcourse 
correction is then 


S TM 


M E M T 
e 


where M is the sensitivity matrix relating midcourse velocity errors to 
target misses in the B ’ T , B * R coordinates. At 5 days after injection, 
the numerical results for the Type I and Type II trajectories are, 
respectively. 


and 


E TM 


-1.616E5 

-1.484E5 


-1.484E5 
1. 370E5 


km 2 


S TM 


2.814E6 
-7. 610E5 


-7. 610E5 
2. 058E5 


km 


2 


The resulting error ellipses are depicted in Figures 4-7 and 4-8. 

4. 5 SUMMARY OF POWERED PERFORMANCE ANALYSES OF THE 
JUPITER MISSIONS 

The Jupiter missions were analyzed for only one value of initial 
azimuth misalignment, viz, 20 arc sec. For this value, the Saturn 
IB /Centaur injection errors are those indicated in Table 4-X for the 
TG-266 system. 


TABLE 4-X. 

SATURN IB /CENTAUR JUPITER MISSIONS 
INJECTION ERRORS (RTN COORDINATES) 


Position 

(km) 

Velocity 
(m/s ec) 

RTN 

RTN 

0.97 1.40 1.89 

4.12 2. 02 5.46 
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B • R 



1600 1200 800 400 


0 ' 400 800 1200 1600 (KM) 


Figure 4-7. Mars Trajectory, Type I (Miss Ellipse 
After First Midcourse) 


B ■ ft 



Figure 4-8. Mars Trajectory, Type II (Miss Ellipse 
After First Midcourse) 
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The resulting uncorrected target miss covariance matrices for the 
solar probe and the cross ecliptic probe are, respectively, 


and 



2. 551E10 

7. 144E9 

T t = 

7. 144E9 

2. 121E9 


" 



2. 328E10 

6. 65E9 

S T ” 

6.65E9 

2. 023E9 


km 


km 


The corresponding error ellipses are shown in Figures 4-9 and 4-10, 
The midcourse 95% velocity requirements at 5 days after injection to 
correct for the injection errors are tabulated in Table 4-XI. 

The corrected target miss ellipses are, respectively, 


v = 

5. 22E4 

1.462E4 

km 2 

TM 

1.462E4 

4. 339E3 

y = 

4.763E4 

1.361E4 

km 2 

TM 

1.361E4 

4. 14E3 


and the corresponding error ellipses are depicted in Figures 4-11 and 4-12. 

4. 6 MARS APPROACH ANALYSES 

The navigation performance analyses for interplanetary and planet 
approach phases conducted under Tasks I and II of this study have been 
modified and extended to the extent summarized below: 

a) New approach trajectories to Mars have been 
incorporated (nominal perifocus altitude of 
1100 km for both Type I and II trajectories, 
see sec. 2). 
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TABLE 4 -XI 


NINETY-FIVE PERCENT AV MIDCOURSE {5 DAYS) 
REQUIREMENTS FOR THE TWO 1971 JUPITER MISSIONS 

(TG-266 SYSTEM) 


Mission 

Correction 

Type 

95% AV 
(m/ sec) 


M + T* 

10. 5 

Solar Probe 




M 

9. 1 

Cross Ecliptic 

M + T 

10. 5 

Probe 

M 

9.1 


Miss plus time -of -flight correction. 
Miss correction, only. 


b) Approach navigation performance was investi- 
gated starting from 50 days before encounter, 
rather than throughout the interplanetary 
trajectory. (Previous results indicated that 
significant orbit determination improvements 

did not occur until about 10 days before encounter. ) 

c) For the doppler plus onboard optical measure- 
ment cases, optical measurements were utilized 
only in the region where the Mars subtense angle 
was between 0. 2 and 3°. (This corresponds to 
roughly the time period from 350 hr down to 11 hr 
prior to encounter for both trajectories, Figure 4-13). 
See par. 4. 6. 3 for further discussion. 

d) The position and velocity errors as given in the 
Tasks I and II report were converted to errors in 
basic approach parameters. 

e) The effects of errors in approach trajectory deter- 
mination and deboost execution errors on total fuel 
requirements were investigated. 
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B • R 


AV EXECUTION ERRORS a - 522.1 KM 
b = 278.6 KM 



600 500 400 300 200 100 0 100 200 300 400 500 600 (KM) 

Figure 4-11. Jupiter, Out of Ecliptic Trajectory (Miss Ellipse 
After First Midcourse) 


b • R 



i i i i i i i i i i i 1 1 

600 500 400 300 200 100 0 - 100 200 300 400 500 600 (KM) 

Figure 4-12. Jupiter, Solar Probe Trajectory (Miss Ellipse 
After First Midcourse) 
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CLOCK AND CONE ANGLES (DEG) 


CLOCK ANGLE TYPE I 



Figure 4-13. Measurement Angles Versus Time to Encounter, 1975 Type I and II Trajectorie 



For this study, five different configurations were considered: 

a) Case 1 — doppler only 

b) Case 2 — doppler plus optical Model A' 

c) Case 3 — ~ doppler plus optical Model B 

d) Case 4 — doppler plus optical Model C 

e) Case 5 — same as Case 4 less Mars subtense 

angle measurement 

The clock and cone angle measurements used with the optical instruments 
are illustrated in Figure 4-14. The models are summarized in Tables 
4-XII through 4-XV. The nonoptical error model given in Table 4-XII 
was modified from that used in the previous study to include a more 
detailed model of the gravitational constant and second zonal harmonic of 
Mars. These parameters were included in the total error model and 
switched in when the vehicle approached the sphere of influence of Mars. 



Figure 4-14, Optical Angle Measurements 

’’’Same as that used in Ref. 4-1. 

^ Model discussed in sec. 5. 
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TABLE 4 -XII 


RADIO/ OPTICAL/ INERTIAL ERROR MODEL MARS MISSION 


Error 

(Variance) ^ 2 

Time Constant 

Random acceleration acting on spacecraft' 
(models uncertainty in the dynamic model 
of the solar systemj i. e. , errors in solar 
pressure forces, gravitational constants, 
etc. ) 

0. 531 x 10 ^m/secp 
(0. 174 x 10 ft/sec ) 
(causes a 200 -km position 
error in 176 days) 

1 week 

Tracking system errors 



• Range rate bias 

10 2 m/sec (0. 0328 ft/sec) 

1/3 day 

• Uncorrelated noise on doppler rate 

0. 732 x 10" 2 m/sec 



(0. 024 ft/ sec) 

(equivalent to 0. 12 ft/ sec 
per 1-sec sample, 25 
meas. averaged) 


Vehicle errors at injection (3 hr) 



• Position 

2 km (6560 ft) 


• Velocity 

2 m/sec (6. 56 ft/sec) 


Mars ephemeris error (relative to Earth) 



• Position 

220 km (7. 22 x 10 2 ft) 


• Velocity 

0. 05m/sec (0. 164 ft /sec) 


Radius of Mars 

20 km (6. 56 x 10 4 ft) 

1 day 

Uncertainty in gravitational constant of Mars 

8.59556 km 2 / sec 2 


Uncertainty in second zonal harmonic of Mars 

0.48 x 10~ 3 



Equivalent error averaged over 25 measurements. 




TABLE 4 -XIII 


OPTICAL ERROR MODEL A 


Error 

(Var iance) 

1/2 

Time Constant 

Sun sensor bias 

1. 746 x 10“ 3 rad 

(6 arc min) 

1/ 2 week 

Sun sensor uncorrelated noise 

0. 349 x 10 4 rad' 
(0. 1746 x 10-3 rad ) 

(0. 12 arc min) 


Mars sensor bias 

1. 92 x 10~ 3 rad 

(6. 6 arc min) 

1/2 week 

Canopus sensor bias 

0. 873 x 10~ 3 rad 

(3 arc min) 

1/2 week • 

Mars sensor uncorrelated noise 

0. 349 x 10“ 4 rad" 
(0. 1746 x 10-3 rad)"' 



Canopus sensor uncorrelated noise 

0. 1746 x 10 “ 4 rad* 
(0. 873 x 10-4 rad ) 



Mars subtense measurement 




Lower limit on 
(variance) ^ ' 2 of bias 

0. 873 x 10 3 rad 

(3 arc min) 

1/2 week 

Error proportional to 
subtense angle' 

1% 


1/2 week 

Uncorrelated noise 

0. 1746 x 10~ 4 rad*..,. 
(0. 873 x 10-4 rad)"'”" 

0. 06 arc min.,..^ 
(0. 3 arc min) 

- 


Equivalent error of 25 measurements averaged. This value was used in the error analysis. 
Single measurement error. 

' 'Percent of subtense angle (a) contributing to the standard deviation of the bias error adding to 
subtense angle, i. e. , <j^ = (0. 873 x 10-3) 2, + (0.01a) 2 . 




TABLE 4-XIV 


OPTICAL ERROR MODEL B 


Error 

(Variance) 

1/2 

Sun sensor bias 

0.407 x 10~ 3 rad 

(1.4 arc min) 

Sun sensor uncorrelated noise 

0. 349 x 10 4 rad ‘ „ 
(0. 1746 x 10-3 rad)"'" 

(0. 12 arc min) 

Mars sensor bias 

0.391 x 10~ 3 rad 

(1.35 arc min) 

Canopus sensor bias 

0. 391 x 10 3 rad 

(1,35 arc min) 

Mars sensor uncorrelated noise 

0.349 x 10 4 rad' 

(0. 1746 x 10" 3 rad)"'" 

- 

Canopus sensor uncorrelated noise 

0. 1746 x 10 4 rad?,..,. 
(0. 873 x i0“ 4 rad)'"'" 


Mars subtense measurement 



Lower limit on 
(variance)*' 2 of bias 

0.485 x 10“ 4 rad 

(0. 17 arc min) 

Error proportional to 
subtense angle' 1 

0% 


Uncorrelated noise 

0. 1745 x 10 “ 4 rad* ; , 
(0. 873 x 10“ 4 rad) 

0.06 arc min 
(0. 3 arc min) 


Time Constant 
1/E week 

1/2 week 
1/2 week 


1/2 week 


Equivalent error of 25 measurements averaged. This value was used in the error analysis. 
Single measurement error. 

Percent of subtense angle {a) contributing to the standard deviation of the bias error adding to 
subtense angle, i. e. , = {0. 873 x 10“3) 2 + (0.01a) 2 . 




TABLE 4 -XV 


OPTICAL ERROR MODEL C 


Error 

(Var iance) 

1/2 

Time Constant 

Sun sensor bias 

0 . 407 x 10 1 rad 

(1.4 arc min) 

1/ 2 week 

Sun sensor uncorrelated noise 

0.349 x 10 4 rad‘ 
(0. 1746 x 10-3 r ad) 

(0. 12 arc min) 


Mars sensor bias 

6. 153 x 10 4 rad 

33 arc sec 

1/2 week 

Canopus sensor bias 

Mars sensor uncorrelated noise 

Canopus sensor uncorrelated noise 

0. 727 x 10 -4 rad 

0. 349 x 10 4 rad‘ . 
(0. 1746 x 10 “3 r ad) ' r 

0. 1746 x 10 4 rad,,..,. 
(0. 873 x 10-4 r ad) "' r 

15 arc sec 

1/ 2 week 

Mars subtense measurement 

Lower limit on 
(variance) 1 / 2, of bias 

0.485 x 10~ 4 rad 

(0. 17 arc min) 

1/2 week 

Error proportional to 
subtense angle 

0% 



Uncorrelated noise 

0. 1745 x 10“ 4 rad’L 
(0. 873 x lO- 4 rad) 

0. 06 arc min ;t , ;1 , 
(0. 3 arc min) 



Equivalent error of 25 measurements averaged. This value was used in the error analysis 
Single measurement error. 







4. 6. 1 Summary of SVEAD Results 


The position and velocity error covariance matrices generated by 
the SVEAD program are expressed in ECI coordinates, with Z along the 
earth's spin axis and X along the vernal equinox. The 1 ~<r values of the 
position and velocity components for the last 10 days of the mission are 
shown in Figures 4-15 through 4-20 for both the Type I and II trajectories. 

One immediate conclusion from these plots is that the results for 
Cases 4 and 5 are almost identical for all components, i. e. , angle sub- 
tense measurement, at least with the accuracies quoted, does not improve 
overall navigation accuracy. Also, except for the X component, there is 
no apparent significant difference in results for Type I versus Type II 
trajectories. 

4. 6. 2 Conversion of SVEAD Results to Approach Condition Errors 

For purposes of further error analysis, it is convenient to express 
the SVEAD error quantities in terms of the orbital approach parameters. 
The pertinent equations are summarized below. 

The transformation from ECI to MCI coordinates is (see Ref. 4-2): 


-MCI fT3 ~ECI 

where 


the MCI coordinate system is chosen with Z 
along the Martian spin axis and X along the 
Martian vernal equinox, x denotes any vector, 
and [T] is the coordinate transformation matrix, 
which for 1975 is numerically equal to 


[T] 


0. 8941E-3 
-0. 9037406 
0. 4279037 


0. 9028834 
-0. 1847545 
-0. 3879994 


0.4297079 
0. 3860002 
0. 8161376 


(Ref. 4-2) 
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TYPE I TRAJECTORY 
TYPE II TRAJECTORY 



TIME TO ENCOUNTER (DAYS) 


Figure 4-15, Position. Uncertainty Versus Time to 
Encounter (Type I and II Trajectories) 


EE 





Figure 4-16. Position Uncertainty Versus Time to 
Encounter (Type I and II Trajectories) 
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Figure 4-19. 


Velocity Uncertainty Versus Time to 
Encounter (Type I and II Trajectories) 






Present radial distance: 

2 2 2 ? 
R - X + Y + Z 

where 


X, Y, Z are obtained from the SVEAD trajectory 
printout. 


Velocity magnitude at R: 


where 


V^ = X + Y + Z 


X, Y, Z are obtained from the SVEAD trajec- 
tory printout. ' 

Velocity direction relative to R: 

XX + YY + ZZ 

COS Y = 


Hyperbolic approach velocity: 

V 2 = V 2 - 2^ 

CO R 


where 


p = Mars' gravitational constant. 
Impact parameter: 


B = R V sin y/V ro 


Semimajor axis, approach angle, and eccentricity of the approach hyper- 
bola: 


a = 



tan v = a/B 
e = esc v 

Present hyperbolic and true anomaly: 

cosh H = (1 - R/a)/e 

tan f = B sinh H/a (e - cosh H) 
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tan X = tan A sin cp 
tan u = tan cp sec A 

Time to periares: 

AT = a {e sinh H - H)/V 

The first-order differentials of the above equations yield the error 
expressions for the usual orbital elements {a, e, i, 0, «), AT), or 
equivalently, (B, V^, i, fi, co, AT): 



The approach orbit determination error covariance matrix is then 

£ od = [ A ] [T] S Pj y [T] T [A] T (4-18) 

where 

,, is the error covariance matrix from SVEAD 
Jr, V 

expressed in ECI coordinates. 
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4.6.3 Conclusions from Approach Navigation Performance 
Analyses 

The results of the conversion indicated in Eq, (4-18) are plotted in 
Figures 4-21 through 4-26 for Cases 1 and 4 (Type I) and Case 4 (Type II). 
Between 50 and 2 hr before encounter, the 1-cr values of hyperbolic approach 
velocity estimation error are essentially identical for the two Type I cases, 
i. e. , the optical angle measurements, although improving the quality of 
component velocity estimation (see Figures 4-18, 4-19 and 4-20), do not 
improve the overall approach velocity estimation. This behavior is attri- 
butable to the different manner in which the errors are correlated in the 
doppler-only and doppler -plus -optical cases. In the former, all errors 
are highly correlated as indicated in the normalized Spy covariance 
matrixes 2 hr before encounter (Table 4-XVI). 

In general, there are no significant differences in the result for 
Type I and Type II trajectories. The choice of the type of trajectory is 
then dictated by mission considerations other than accuracy. 

Another important conclusion to be made from this analysis is that 
in the time span from about 5 hr down to 2 hr prior to encounter, the 
position uncertainty as given by Og (Figure 4-21) is not significantly 
improved by the addition of -optical measurements to doppler measurements. 
That is, the added complexity of a planetary approach sensor to the ROI 
guidance system may not be warranted if only the size of the orbits to be 
established about Mars is of concern (see also par. 4. 6. 3). However, if 
it is desired to accurately establish cross plane parameters, viz, orbital 
inclination and longitude of the ascending node for the approach phase, the 
planetary approach sensor can be useful (see Figures 4-23 and 4-24 and 
par. 4.6.4). 

Additional onboard optical measurements during the last 11 hr would 
result in only minor improvement in knowledge of the approach trajectory 
characteristics as indicated by Figures 4-21 through 4-26. Tracking was 
thus terminated at this point to simplify the design of the planet tracker 
(see subsec. 9. 7) and to provide the necessary time for ground compu- 
tation and transmission of corrections to the orbit insertion maneuver - 
based on the doppler tracking data and telemetered onboard optical 
measurements . 
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UNCEKIA1NIY IN APPROACH VELOCITY IM/SEC) UNCERTAINTY IN IMPACT PARAMETER, 8 (KM) 




(DEG) ' <7|(DEG) 



HOURS BEFORE ENCOUNTER 

Figure 4-23. Orbital Inclination Standard Deviation, 
1975 Type I and II Trajectories 



HOURS BEFORE ENCOUNTER 

Figure 4-24. Ascending Node Longitude Standard Deviation, 
1975 Type I and II Trajectories 
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TABLE 4-XVI. 


NORMALIZED COVARIANCE MATRICES AND i-ff COM- 
PONENT VALUES AT 2 HR BEFORE ENCOUNTER 


Doppler only (Case 1): 

°X = 56.04, cr^ = 72.78, cr z 
cr^ = 0.801, oy = 1.603, cr^ 

= 455.4 
= 3.962 

km 

msec 

i -0.605 

-0.687 0.327 

-0. 749 

0.687 

1 

0.994 -0.949 

0. 980 

-0. 994 


1 -0.911 

0. 996 

-1 

Symmetric 

1 

-0. 870 

0. 910 



1 

-0.996 




1 

Doppler + best optic 

s (Case 4) 



cr x = 34. 59, 

0y “ 7*12; (T^ 

= 31.39 

km 

cr^. = 0.264, 

o-y = 0.187, 

= 0.275 

msec 

1 0.647 

-0.356 -0.957 

-0. 846 

0.366 

1 

0.436 -0.805 

-0. 165 

-0.428 


1 0. 107 

0. 768 

-0. 981 

Symmetr ic 

1 

0. 688 

-0. 12 



1 

-0. 775 




1 


102 








4.6.4 Effects of Approach Orbit Determination and Execution 
Errors on In-Plane Performance 


On the basis of the nominal values of the impact parameter B and 
hyperbolic approach velocity V o determined during the approach phase, 
the velocity increment ^^jeboost’ ^OT^ed to establish the 1100 x 
10,000 km (periares and apoares altitudes) orbit, can be determined from 
the following’'' (refer to Figure 4-27) : 



NOTES: 1. v AND r ARE PREDICTED APPROACH PERIARES, 

’ DISTANCE AND VELOCITY . 

2. V IS VELOCITY REQUIRED AT r TO REACH 
P P 

DESIRED APOARES R q 

3* AV dboost = v p - V p APPLIED ANTiPARALLEL 
TOv 

P 

Figure 4-27. Approach Parameters 


*It is assumed that the desired orbits about Mars are in the same plane 
as the approach orbit plane. 
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Semimajor axis of approach hyperbola: 


a = 



Approach angle : 

tan v = a/B 

Eccentricity of approach hyperbola: 


Periares distance: 


Velocity at periares: 


e = esc v 


r = a(e - 1) 
P 


V = BV /r 
P 00 P 


Velocity at periares required to attain desired apoares R^; 

9 R /r 

V 2 = 2 p 5 - - P 
p ^ R_ + r 


a p 


Required deboost velocity: 


AV , , = v - V 

deboost p p 

The time from some reference point until AV , , , is to be applied can 

c deboost rr 

be obtained from 


where 


T = dr- (e sinh H - H ) 

V co o O 


cosh H “ (1 + r /a)/e 
o o 

r = reference distance, 
o 


It is apparent from the above that errors in the determination of 
B and V^ will propagate to an error in the magnitude of an d 

nontangential application of AV^^^^ as illustrated in Figure 4-28. 


• 104 




NOTES: 1 . R, V R ARE ACTUAL DISTANCE AND VELOCITY 

PREVAILING AT PREDICTED TIME OF V DEBOOST 
APPLICATION 

2. V IS RESULTING VELOCITY AT R 

Figure 4-28. AV^ ^ Application Error Geometry 


In Figure 4-28, the angles of interest are 


Av — rotation of the line of apsides of the approach 
hype rbola 


f — true anomaly angle at which AV , , 

■T • , • • i » -» Cl 6 D O O S t 

application is applied 




flight path angles before and after AV 
application at (R, f) 


deboost 
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The rotation of the line of apsides of the approach hyperbola due to errors 


in estimation of B and V_ is 

CO 


AB AV ° 

Av = -sm v cos v — g- -2 sm v cos v -y— 


For the nominal approach conditions of 


B = 7829 km, 

- 3. 0922 km/sec, and 
v = 29° 50' , 

| Avj is less than 1° for [ AB | < 200 km and | AV^ j < 2 msec. 

The true anomaly angle measured from the true line of apsides is given 
by 


f _ ^/(e 4- Ae)^ - 1 sinh H 
an (e + Ae) - cosh H 

where 

(V a +AV w ) 

(e + Ae) sinh H - H = 8(AT ) — ^ — 


H = hyperbolic anomaly on the true trajectory at the 

time of predicted perifocus passage on the assumed 
trajectory 

6 (AT) = difference in perifocus passage time for the true and 
assumed approach trajectories 


Ae, Aa = the error in predicted approach hyperbola eccentricity 
and semimajor axis due to AB and AV m . 

The flight path angle prior to ^Vd e k oos t- a Ppli ca -ti on is 


sin y 


(B + AB) (V M + AV m ) 
RV 
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where 


V 2 = (V + AV ) 2 + ^ 

R = p/ [l + (e + Ae) cos f] 

p = (e + Ae + 1) (r + Ar ) 
V P P 


Following A¥ deloost 
V sin r = V 

v cos r = v 


application, 

, sin y - (AV 
cos y - (AV 


deboost 

deboost 


+ 8 V) cos (Av+ f + 7] ) 
+ 6 V) sin (Av + f + Ip 


where 


6 V = magnitude of execution error 
T] = execution direction error 

The various parameters of the resulting elliptic orbit include the following 
Semimajor axis: 

a Q = R/(2 - V R 2 R/jx) 

Eccentricity: 

2 , „ 2, r 2 . 2 „ . 
e = 1 - R V_ sin T a u. 
o R o' 

Periares and apoares: 

r = a (1 - e ) 
p o o 

r = a (1 + e ) 
a o o 

Rotation of line of apsides: 

cos(Ai|j) = ((r r a / a 0 ) " R )/ e 0 R 

Thus the orbit established by the application of -^ as 

following errors from the desired: 

Periares error: 

6 r = r' - R = r' - (1100 + 3410) km 
P P P P 
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Apoares error: 


5r = r - R = r - (10, 000 + 3410) km 
a a a a 


Rotation of the line of apsides: 

Aip = cos 1 {({r r /a ) - R)/e R) 

P «. U U 


For this study, it was assumed desirable to correct to the nominally 
desired 1, 100 x 10, 000-km orbit with the desired direction of the line of 
apsides. From par- 2. 4. 4 of the Tasks I and II report, the following 
velocity increments are required to make these corrections: 


Apoares correction: 


A V 

a 



R /R 
a p 

R + R 
a p 



Feriares correction: 


AV b 



R /R 
a' p 

R + R 
a p 



6 r 


a 


_RJ 

+ V 


Rotation correction: 


AV c ‘\pprt - V R a> sin ^ 

’ a p ^ 

Substituting the desired numerical values, 

AV = 0. 14 [fir m/sec 

a I a i 

AV, = 0. 105 ISr | m/sec 

b j p| 

where 

6r and 6r are expressed in km, and 
Si p 

AV = 2490 sin m/s ec 

c L 
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A Monte Carlo simulation was run to determine the effects of AB, 
AV , .<5v, and 7| on Sr , <5r , AiJj, and d(AV), where S(AV) is the 

03 3- p 

difference in total velocity increment requirements from the ideal nominal 
case. For this case, the total AV is that required to deboost into the 
1100 x 10, 000-km orbit and subsequently transfer down to and circularize 



where 


R = desired circular orbit radial distance, 
c 

For the parameters of interest, 


AVppQjp^.^ = 1, 577 + 0. 07 + 0. 80 = 2. 45 km/sec 

The Monte Carlo runs were made not only for the nominal approach impact 
parameter value of B = 7 829 km, but also for high- and low-approach cases 
of B = 8029 and 7629 km, respectively. For these cases, 6 (AV) is 
affected further by the fact that are 599 and 1. 554 km/ sec, 

respectively, as compared to 1. 577 km/ sec. 

The results from the Monte Carlo runs are shown in Figures 4-29 
and 4-30. The 1-cr values used for the runs were 

10 < Ug < 100 km 

CTy = 0. 1 m/sec 

CD 

<Tgy = °s = 9^3% of applied AV 

cr„ = cr =0. 001 rad 
11 P 
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CUMULATIVE DISTRIBUTION OF DISPERSION 
IN ORBITAL TRIM AV 



CHANGE IN AV FROM NOMINAL (MSEC) 


& Figure 4-29. Cumulative Distribution of Dispersion 

in Orbital Trim AV 


The effects of these errors on the change in total AV requirements, 
i. e., the AV required for trimming the deboosted elliptic orbit to the 
desired 1, 100 x 10, 000-km orbit, are illustrated in Figure 4-29 in the form 
of cumulative distributions. The 95% AV requirements are summarized 
in Table 4-XVH. 
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TABLE 4-XVII. 

NINETY-FIVE PERCENT ORBIT TRIM AV REQUIREMENTS (M/SEC) 


Approach 

<7-0 = 10 km 

-D 

50 km 

100 km 

High 

85 

125 

185 

Nominal 

45 

85 

145 

Low 

45 

75 

135 


If the final approach orbit estimation is delayed until approximately 
2 hr before encounter, then Og - 55 km for the doppler -only case and 
35 km for the doppler -plus -be st- optics case (Type I trajectory). Inter- 
polation of data in Table 4-XVII indicates that the reduction in orbital 
trim requirements afforded by the use of the optical measurements would 
be slight and would not warrant the added complexity of the onboard 
planetary approach sensor. 

Furthermore, examination of Figure 4-30 indicates that the 
parameters of the deboosted orbit do not differ significantly from those 
nominally desired. For (Tg = 50 km, the maximum differences are 
shown in Table 4-XVHI. For these magnitude errors, orbital trim would 
most likely not be applied. The only change in AV requirements would 
be caused by transfer from the non-nominal elliptic orbit to the 500-km 
circular orbit. In this event, the maximum change in AV would be less 
than ±50 msec (out of 2.45 km/ sec) for any case. 

Thus, if only in -plane orbital parameters are of concern, the use 
of a planetary approach sensor would be difficult to justify. Consideration 
of cross-plane characteristics shows that the approach sensor might 
possibly be a valuable adjunct to doppler tracking orbit determination in 
certain missions. 
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TABLE 4-XVIII 

MAXIMUM ERRORS IN DEBOOSTED ORBIT PARAMETERS 


Approach 

Apoares 

(km from 10, 000) 

Periares 
{km from 1100) 

Rotation 

(deg) 

High 

-100 to +200 

+150 to +230 

3. 8 

Nominal 

-100 to +200 

-35 to +34 

3.8 

Low 

-100 to +200 

-130 to -210 

3.8 


4. 6. 5 Effect of Approach Orbit Determination Errors 
on Out-of-Plane Performance 

Referring to Figure 4-23, the uncertainty in the inclination of the 
approach orbit at 2 hr before encounter is seen to be about 2. 5° and 0. 25° 
for the doppler-only and doppler -plus -approach- sensor cases, respec- 
tively. If an inclination correction is attempted simultaneously with the 
deboost AV application, subsequent tracking in orbit would reveal the 
inclination error of the approach orbit determination. The corrective 
AV required to adjust the inclination error is 


AV. 

x 



R /R 

a p 

R 4- R 
a p 


/ 


1 - cos Ai 


{Ref. 

par. 


4-1, 

2. 4.4) 


= Ai 


R /R 
a p 

R + R 
a p 


where 


: 2665 Ai m/sec 


Ai is the negative of the inclination uncertainty 
of the approach determination function. 

For the two cases considered, the AV. penalties are 

116 m/sec; doppler only 

AV. = 

11.6 m/sec; doppler + approach sensor 
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Thus the approach sensor can provide significant benefits if either 


a) precise inclination information is desired 
very early, or 

b) approach inclination corrections are to be 
applied simultaneously with the deboost AV. 

However, if neither of the above is a requirement, then it would be better 
to doppler track while the payload is in orbit about Mars and forego the 
approach sensor. 

4. 7 MARS AND LUNAR ORBIT DETERMINATION 

4. 7. i Mars Orbit Determination from DSIF Tracking Data 


The accuracy of orbit determination while the spacecraft is in an 
areocentric orbit was obtained using the SVEAD computer program (see 
Ref. 4-1, app. D) for the nominal 1, 100 x 10, 000-km orbit obtained from 
the Type I heliocentric transfer orbit. The orbital characteristics are 
shown in Table 4-XIX. 


TABLE 4-XIX 

ORBITAL PARAMETERS FOR MARS ORBIT 
(MARTIAN EQUATORIAL COORDINATE SYSTEM) 


Parameter 

Value 

Semimajor axis (a) 

8960 km 

Eccentricity (e) 

‘ 0. 496 

Inclination (i) 

3 6. 6° 

Longitude of 
ascending node (£2) 

143. 1° 

Argument of 
perigee (co) 

-12.3° 

Period (T) 

7. 15 hr 


For the 1, 100 x 10, 000-km orbit obtained, the spacecraft goes behind 
Mars 9 min after periapsis and is visible again 33 min later. 

The initial state vector errors used in the analysis were those ob- 
tained at the end of the approach orbit determination phase (doppler 
tracking only). All other error models were the same as those used in 
the approach orbit determination phase (refer to Table 4-XII). 
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The resulting behavior of the unc ertaintie s in the spacecraft position 
and velocity in RTN coordinates are illustrated in Figures 4-31 and 4-32 
for slightly more than one complete orbit. The corresponding orbital 
elements are illustrated in Figures 4-33 and 4-34. 

3 ? 

This analysis assumes an uncertainty of 8. 6 km /sec^ in the gravi- 

_3 

tational constant of Mars and an uncertainty of 0.48 x 10 in the second 
zonal harmonic, as indicated in Table 4-XII. The results shown in Figures 
4-31 through 4-34 indicate approximately an order of magnitude reduction 
in the initial errors over a period of one orbit. These results are valid 
only if no significant local gravity anomalies or other unknown disturbing 
accelerations are present. The only method of validating this assumption 
is by analysis of actual tracking data obtainedfor a spacecraft in Mars orbit. 

The maximum allowable 3c r orbit determination errors, specified at 
orbit insertion plusfour returnsto periapsis are (from Table 3- VII, voll): 

Semimajor axis 3cr =10 km 

C 1 

-4 

Eccentricity 3tr =10 

Uncertainty in time of 3cr = 5 sec 

periapsis in passage ^ 

Comparison of the orbital determination errors given in Figure 4-33 
with these values shows that the errors have been reduced to approxi- 
mately one-third of the desired values after one orbit. Although the 
specific analysis was not carried out, it is expected that the desired ac- 
curacies would be achieved after four orbits. 

4.7.2 Lunar Orbit Determination From MSFN Tracking Data 

It has recently become clear that the orbit-determination errors 
for low trajectories around the moon deviate significantly from 
previous predictions (Ref. 4-4). A recent report published by JPL 
(Ref. 4-5) concluded that the moon is gravitationally rougher than anti- 
cipated in the sense that comparatively high-degree terms in the spherical 
harmonic expansion would be required for effective representation of the 
lunar gravity field. This conclusion is based upon an extensive investi- 
gation of the radio tracking data for the Lunar Orbiter Missions received 
by the NASA Deep Space Network. 
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(KM) 



Figure 4-31. Position Uncertainties Versus Time for 

Spacecraft in Mars Orbit (DSIF Tracking) 
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Figure. 4- 33. Uncertainties in Orbital Parameters Versus’ Time . 

for Spacecraft in Mars Orbit (DSIF Tracking) 
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(DEG) 



HOURS FROM FIRST PERIAPSIS 

Figure 4-34. Uncertainties in Orbital Orientation Parameters Versus Time 
for Spacecraft in Mars Orbit (DSIF Tracking) 




The effect of the observed gravitational anomalies on orbit- 
determination accuracies is under intensive investigation at the present 
time in support of the Apollo program. The results of studies completed 
to date, while preliminary, give a reasonably good estimate of the orbit- 
determination capabilities under various circumstances. 

Table 4-XX shows the estimated navigation uncertainties based on 
the Lunar Orbiter ground tracking data postflight analysis for one pass 
of tracking data {from Ref: 4-6]. 

Table 4-XX indicates that the uncertainties in the position compo- 
nents of the spacecraft state vector grow rapidly with time after tracking 
is terminated. 


TABLE 4-XX 

LUNAR ORBIT NAVIGATION UNCERTAINTIES 



At Time of Tracking 

Rate of Growth 

Propagated for 
Two Orbits 
(no tracking) 

^R 

1000 ft 

1200 ft/rev 

2600 

^T 

3000 

4000 ft/rev 

8544 

°N 

300 

500 ft/rev 

1044 

% 

0. 1° 

0 . 005°/r ev 

0.11 ° 


0.02° 

0 

0.02° 

°R 

7. 3 ft/sec 

- 

7. 3 ft/ sec 

T 

2. 2 ft/sec 

- 

2. 2 ft/sec 

°”n 

9 . 2 ft/sec 

- 

9 . 3 ft /sec 


Note : 

o' ^ = error in flight path angle 

o'q = error in direction of angular momentum vector. {See par. 1.3. 4 
for definition of RTN coordinate system) 
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Figure 4-35 (from. Ref. 4-7) compares a trajectory calculation 
using a triaxial gravity model with a Lunar Orbiter 3 (L03) trajectory. 

The L03 trajectory was generated by processing the nondestruct doppler 
measurement data by the modified Apollo Powered Flight Data Processing 
Program (APOP) (Ref. 4-8). The data were obtained at Goldstone, 

Woomera and Guam. Since the modified APOP' can model a large random 
uncertainty in the gravity model by treating it as a system noise, the 
trajectory generated by the modified APOP is free to deviate from that 
predicted by the triaxial gravity model and thus closely approximates 
the actual trajectory. 

Comparison of the results shown in Figure 4-35 with the uncertainties 
given in Table 4- XX shows reasonable agreement for the inplane compo- 
nents. Figure 4-35 also shows somewhat larger out-of-plane position errors. 
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(b) Downrange Position Error Versus Time 



(c) Crossrange Position Error Versus Time 


Figure 4-35. Comparison of L03 Trajectory Data with Trajectory 
Computed Using T'riaxial Gravity Model 
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5. CONTROL SYSTEM PERFORMANCE ANALYSES 


5.1 INTRODUCTION AND SUMMARY 

Control system analyses made in support of the conceptual and 
modular design of the Radio/ Optical/Strapdown Inertial Guidance and 
Control system are presented in this section. Several vehicle configura- 
tions were evaluated in terms of their bending transmis sibility in addition 
to performance of detailed stability analyses. On the basis of these studies, 
suitable locations of the rate gyro packages in the launch vehicles have 
been determined. Filtering requirements for stability compensation also 
have been determined and used to formulate the general digital control 
equations presented in sec. 8 for the modular guidance and control system 
design. Control system stability criteria followed in this study are given 
in Appendix A. 

Gain-phase stability plots, control system gains, and relative stabil- 
ity margins are given for alternate rate gyro locations and compensation 
filters for the Atlas/ Centaur, Saturn IB/ Centaur, launch vehicles and for 
the Voyager Spacecraft during midcourse and Mars orbit insertion burns. 

In addition, a spacecraft thrust vector pointing error analysis was per- 
formed for the Voyager spacecraft to estimate the TVC pointing errors 
and crossrange velocity errors incurred during the powered burns. 
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5. 2 CONTROL SYSTEM STABILITY AND BENDING MODE 
TRANSMISSIBILITY COEFFICIENTS 

In the design of the boost vehicle first stage, minimum gains for 
control system stability are dictated primarily by the aerodynamic 
moment acting on the vehicle and the control moment available to counter- 
act it. Given a vehicle configuration, its trajectory, and the desired 
aerodynamic stability margin, the control gains are constrained by a 
determinable lower limit. 

Similarly, an upper constraint on the control gains exists dictated 
primarily by the effects of vehicle structural resonance modes. When the 
structural bending sensed by the control system sensors becomes exces- 
sive, the control system gains may be over -constrained such that vehicle 
stability with conventional control systems or, for that matter, any con- 
trol system becomes impossible. 

An illustration of the stability problem is given by Figures 5-1 
and 5-2. In the first figure, a typical gain-phase plot is shown in which 
the control system frequency response and stability margins are indi- 
cated. The desired operating point is defined by the intersection of the 
dashed nominal gain line and the gain axis. The response of the control 
system to increasing frequencies, indicated by the arrows, is shown to 
encircle the operating point, producing an island of stability. The 
associated stability margins are termed: aerodynamic or low-frequency 
rigid body gain margin (A), rigid body phase margin (C), high-frequency 
rigid body gain margin (B), and the first bending mode phase margin (D). 

As the first bending mode transmissibility increases, an erosion of 
the high-frequency rigid body gain margin ensues as shown in Figure 
5-2, The increase in bending transmissibility is shown to dissolve the 
island of stability to a point where even the addition of stability compensa- 
tion networks would be of no avail. It is estimated that a desirable trans- 
missibility would be less than 10 and certainly no greater than 30, this 
latter case requiring compensation networks which include notch filters. 

A notch filter could be designed at the high-frequency rigid body fre- 
quency thereby producing attenuation and increased gain margin; however. 
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GAIN (DB) 



A LOW FREQUENCY AERODYNAMIC GAIN MARG IN 
B HIGH FREQUENCY RIGID BODY GAIN MARGIN 
C RIGID BODY PHASE MARGIN 

D 1 ST BENDING MODE PHASE MARG IN (LOW FREQUENCY) 

E 1 ST BENDING MODE PHASE MARG IN (HIGH FREQUENCY) 

F 2ND BEND ING MODE GAIN MARG IN (IRRESPECTIVE OF PHASE) 

G PHASE TO THE PEAK OF 2ND BENDING MODE RESPONSE 
H 3RD BENDING MODE GAIN MARGIN (IRRESPECTIVE OF PHASE) 

I PHASE TO THE PEAK OF 3RD BENDING MODE RESPONSE 

Figure 5-i. Typical Gain- Phase Diagram With Margins 



Figure 5-2. Effects of Bending Mode Transmissibility Upon 
Stability Margins Without Compensation . 
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its use is constrained by its associated phase shifting effects in which 
increased phase lag and phase lead are accrued about the design frequency. 
This produces an erosion of both ridid body and first mode bending phase 
margins, an erosion which may not be tolerable. 

An effective and normally used method to maintain low transmis- 
sibility coefficients, thereby easing the control compensation problem, is 
the location of control system sensors on the vehicle such that structural 
bending effects are minimized. Location of the rate gyro package on or 
near the first bending mode antinode to accomplish this minimization is 
a common practice. Location of the angular position sensors is less 
critical since a smaller contribution to bending transmissibility is usually 
associated with this feedback loop. In concurrence with the purposes of 
this design study, the angular position sensors were assumed located in 
the spacecraft atop the boost vehicle. 

Several vehicle configurations were evaluated in terms of their 
bending transmissibility in addition to the performance of detailed 
stability analyses. Figure 5-3 shows bending transmissibility time 
histories of three vehicle configurations. The gains assumed were those 
employed on present vehicle programs. Gain changes were not included 
although typically for these liquid propellant boosters, a decrease in gain 
by a factor of two or three is usually executed at approximately 110 sec 
after liftoff. In the figure, the effects of locating the rate gyro package 
at three separate vehicle stations are shown for the Saturn V and Atlas/ 
Centaur vehicle and at two stations for the Saturn IB/Centaur vehicle. 

The unusual drop in transmissibility at mid first stage flight for the 
Atlas/Centaur vehicle results from the abrupt change in the structural 
damping values which were obtained from past design study reports. In 
actuality, a representative smooth curve is to be expected. For the 
Atlas/Centaur vehicle, a rate gyro located at station 708 appears desir- 
able, a location at station 583 (possibly the present location) appears 
tolerable, and location within the spacecraft appears intolerable. 

For the Saturn v/Apollo vehicle, location of the rate gyro package 
in the aft compartment of the S-IVB stage appears desirable, location in 
the instrument unit appears tolerable, and location in the spacecraft 
appears intolerable. 
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Figure 5-3. First Stage Booster Bending Transmissibility 
Coefficients (Position Sensors in Spacecraft) 

For the Saturn IB/ Centaur vehicle, location of the rate gyro package 
in the spacecraft, IU, or aftS-IVB compartinent appears satisfactory with 
preference given to the aft S-IVB location. 

Approximate filtering regions associated with transmissibility 
coefficient values are indicated on the left side of Figure 5-3. The 
lower part of each region may require use of only one filter whereas, for 
the upper part of each region, a filter change during the flight may be 
required. In the notch filter design, additional complex pole filters are 
assumed included for attenuation of higher bending modes. The nonlinear 
filtering region, shown -unbounded, also has limits above which stability 
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is impossible. Indications are that this limit may occur at a transmis- 
sibility less than sixty. 

The significance of the structural damping value and its influence 
on control system compensation is evident, and precise values from 
ground tests are preferable to the use of assumed values. A value of 1% 
was employed (with the exceptions noted) following the procedure that has 
been used in the past. 

The transmissibility coefficients for the Atlas/Centaur higher modes 
are shown in Figure 5-4. Stabilization of the higher modes is more 
readily achieved as indicated by the lower coefficients. Satisfactory sec- 
ond mode stability margins for the design using spacecraft rate gyros 
will be more difficult to achieve as is shown in subsequent discussions. 

No attempt was made to relate higher order transmissibility values to 
filtering requirements although probably such a relationship can be 
formulated. 

The first and second mode transmissibility coefficients for the 
Saturn v/Voyager vehicle are given inFigures 5-5 and 5-6, respectively. 



'Figure 5-4. Atlas/Centaur Higher Mode Transmissibility 
Coefficients 
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Figure 5-5. Saturn V/Voyager First Bending Mode Transmissibility • 
Coefficients Without Filtering 
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Figure 5-6. Saturn V/Voyager Second Bending Mode Transmissibility 
Coefficients Without Filtering 




Location of the rate gyros in the aft S-II compartment is near ideal; in 
the aft S-IVB compartment acceptable; in the ITJ marginal; and in the 
spacecraft intolerable. Location in the IU is considered marginal rather 
than unacceptable in view of the gain change that would occur near 110 sec 
of S-IC flight time. This design condition will, however, pose a consider- 
able filtering challenge. The stabilization of the second mode will also 
prove to be difficult if the rate gyros are located within the IU. 

Two other vehicle configuration transmissibility coefficients are 
shown in Figure 5-7 for the Saturn IB/ Apollo CSM and in Figure 5-8 for 
the Atlas S LV3/Burner 2. The sharp contrast in transmissibility coef- 
ficient between the rigid 9 -tank-clustered Saturn IB vehicle and the highly 
elastic balloon-designed Atlas vehicle is clearly evident. 

5. 3 STABILITY ANALYSIS OF THE ATLAS/ CENTAUR VEHICLE 

The equations employed in the booster stability analyses are given 
in app. B. The analysis data employed are given in app. C for the 
Atlas/ Centaur AC- 13 vehicle configuration. 

The first step of the AC-13 vehicle design study was to obtain a 
baseline for design comparison purposes by investigating the stability 
margins of the existing design. The available data placed the rate gyro 
location at station 583, the position gyros within the booster stage com- 
mand, and the control pod at station 99 L Since an autopilot design for 
this new vehicle configuration was not available, the AC -5 vehicle auto- 
pilot parameters, consisting of complex poles at 15. 5 rad/ sec with pole 
damping of 0. 5 and a position gain of 1. 8, were employed. The analog 
control system gain-phase plots at five booster flight times (0, 40, 80, 

112, and 154 sec) are given in Figures 5-9 through 5-13. In the first 
figure, no stability margins for the first bending mode exist and 
stable slosh modes are shown. In the 40 -sec case, a small first 
mode bending margin exists and an Atlas LOX-slosh instability occurs, 
which is typical for the Atlas/Centaur vehicle. (This instability exists 
over a brief duration and is not usually of importance. ) In the 80-sec 
case, the bending margins are considerably improved and the Atlas LOX- 
slosh mode is stable. In the 112-sec case, the bending margins continue 
to improve, and problems with the Centaur LOX-slosh mode loom. In 
the booster engine cutoff (BECO) case at 154 sec, the bending stability 
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Figure 5-7. Saturn IB/ Apollo CSM First-Stage Booster 
Bending Transmissibility Coefficients 
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Figure 5-9. Gain -Phase Plot Atlas/ Centaur A C= 13 
T=0 Analog Autopilot 
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Figure 5- 10. Gain-Phase Plot.Atlas'/ Centaur A C- 13 
T=40 Analog Autopilot 
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Figure 5-11. Gain-Phase Plot Atlas/ Centaur AC- 13 
Max Q T=80 Analog Autopilot 
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Figure 5-12.. Gain-Phase Plot Atlas/ Centaur AG- 13 
’T= 122 Analog Autopilot 
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Figure 5-i3. Gain-Phase Plot Atlas/ Centaur AC- 13 
T=154 Analog Autopilot 
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margins are ample and the interaction between the Atlas and Centaur 
BOX- slosh modes are heightened; therefore resulting in slosh margin 
reductions. A reduction in control gain to produce some slosh margin 
appears desirable. Figure 5-14 shows the rigid body gain axis crossover 
points plotted as a function of time. The desirable gain schedule is shown 
to involve linear gain variations with time, which will require a rather 
complex implementation ‘compared with the simple gain changes normally 
used. 

Instead of employing linear gain variations, a filter change used 
during the first stage of flight would be a simpler design implementation. 
The effects of changing the complex pole filter frequency to 1 1 rad/sec, 
with the filter pole damping maintained at 0. 5, are shown in Figures 5-15 
through 5-17 for the liftoff, maximum q, and BECO cases, respectively. 
Bending stability at liftoff is obtained at the expense of slosh instabilities 
later in flight. Combining these results with the original filter margins 
produces the composite gain axis crossover curves shown in Figure 5-18. 
A filter change at 70 sec of flight appears desirable, and a single gain 
change at 110 sec appears to produce compromising gain margins. Mar- 
gins in excess, of 6 dB are desired and the 4 to 5 dB level of high- 
frequency crossover margins is considered rather skimpy, but tolerable, 
for flight-proven vehicles. Certainly, improvements in the design mar- 
gins are needed. - 

The location of the rate gyros for the AC-13 vehicle is not optimum 
for control stability, but the possible movement of this package is dis- 
cussed subsequently. After evaluation of the analog AC-13 control sys- 
tem, which is employed as a baseline design, the effects of substituting 
a digital control system were investigated. The ROI guidance system 
was located in the spacecraft on top of the Centaur and three distinct rate- 
sensing configurations were studied: 1) use of the present Atlas rate 
gyros at station 583, 2) use of the ROI rate gyros in the spacecraft, and 
3) use of the Atlas 'rate gyros at station 708 {the optimum first mode 
stability location). The effects of a 25-samples/ sec computer sampling 
rate and a 10 -msec computational delay were included in all of the analyses 
performed. The roots of the equation matrix given in app. A were 
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FLIGHT TIME (SEC) 


Figure 5-14. Atlas/Centaur AC- 13 Rigid Body Gain Margins, 
Analog Autopilot, and Complex Poles at 15. 5 
Rad/Sec 
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ATLAS/CENTAUR AC-13 T = 0 ANALOG AUTOPILOT WF = 1 1 RAD/SEC ’ ■ ' 
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Figure 5-15. Gain-Phase Plot Atlas/ Centaur AC- 13 T=0 
Analog Autopilot WF =11 Rad/ Sec 


141 
















Gain-Phase Plot Atlas/ Centaur AC-13 
T=80 Analog Autopilot WF =11 Rad/ Sec 


142 




Figure 5-17. Gain-Phase Plot Atlas/ Centaur AC- 13 

T= 154 Analog Autopilot ¥F =11 Had/ Sec 
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-FLIGHT TIME (SEC) • • 

Figure' 5- 18, Atlas/ Centaur AC- 13 Rigid Body Gain Margins, Analog 
Autopilot, ’arid Complex Poles changed at T '= 70 



transformed first into the Z -plane and then into the W-plane in which 
control compensation was added. Gain-phase plots of the open-loop 
frequency response were obtained in the W-plane and appear similar to 
those for the analog system. 

With the rate gyros located at station 583, the output signals routed 
upstage to the ROI guidance system, and the use of a digital complex pole 
filter equivalent to the analog filter, the corresponding gain -phase .plots 
for five times of flight are given in Figures 5-19 through 5-23. Compared 
to the analog counterpart of the first bending mode shown in Figure 5-9, 
more phase stabilization is obtained because of the relocation of the posi- 
tion sensors from aft to forward of the antinode of the first bending mode. 
This trend persists through the other four plots, though the slosh mode 
responses remain unaffected by the design change. 

The rigid body gain margin plot of the desirable position gain 
schedule is shown in Figure 5-24. To maintain gain margins above 4 dB, 
a two-step gain schedule is required. The use of a digital notch filter in 
addition to the complex pole filter was next investigated. Two notch 
filter designs at V = 0. 216 rad/sec and V - 0. 174 rad/ sec corresponding 
to S-plane frequencies of 1 1 rad/ sec and 8. 8 rad/ sec were included in the 
respective gain-phase plots of Figures 5-25 and 5-26 for the maximum q 
condition. Damping of the zeros was selected to be 0, 2, and 0. 3 for the 
poles. Additional lead also was included to counteract the phase lag 
effects of the notch, resulting in a net 3 dB of attenuation at the notch 
frequency. The first of the two designs is more desirable since the notch 
lies nearer the gain axis, and with the use of the added notch filter, a 
single gain change may suffice. The other flight times were not investi- 
gated since good stability margins are to be expected from employing a 
filter change which can be readily implemented within the digital computer. 

The second ROI guidance design investigated discarded the Atlas 
rate gyros and used the strapdown rate gyro information for the control 
process. The placement of the control rate sensors at the top of the 
vehicle is probably the worst location because the bending transmis sibility 
is increased by two to three fold over the already marginal baseline design. 

The liftoff and maximum q (T = 80) cases are shown in Figures 5-27 
and 5-28, respectively. A complex pole filter at V - 0, 216 rad/sec 
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Figure 5-19. Gain-Phase Plot Atlas/ Centaur AC - 13. Digital Autopilot 
T=0 Atlas Rate Gyros Run 1 
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Figure 5-20. Gain-Phase Plot Atlas/ Centaur AG- 13 Digital Autopilot 
T=40 Atlas Rate Gyros 
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Figure 5-21. Gain-Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot- 
Max Q T=80 Atlas Rate Gyros 
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Figure 5-22. Gain-Phase Plot Atlas/ Centaur AC-13 Digital Autopilot 
T=112 Atlas Rate Gyros Run 
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Figure 5-2 3 .. „.Gain- Ehas e „Elo t, A tla s /..C entau r .. A.G-- 1,3.. Di gltal' A.uto.pilo.t. 

T=154 Atlks Rate Gyros P- 


150 

t . 






GAIN (DB) 





: F-igure-5-25. Gain-Phase Plot Atlas /Centaur AG- i 3 DigitalAutopilot' 
*■»-' a- - T=80 Atlas Rate Gyros WF = 15.-5‘ + Notch 
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Figure 5-26. Gain-Phase Plot Atlas./ Centaur AC-13 Digital Autopilot 
T=80 Atlas Rate Gyros Notch 8.8 
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Figure„5-27._. Gain-Phase .Plot A tla s / Centaur- -A-G - 1-3 Di-gita-l -Autopilot 
T='0 Spacecraft'. Sensors WP •= ' 1 1' 
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Figure 5-28. Gain-Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 
T=80 Spacecraft Sensors WF =11 
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(equivalent to S- plane 11 rad/ sec) with 0. 5 damping is included in these 
responses. The high-frequency rigid body gain margin is shown to be 
negative or nonexistent for the maximum q condition, and, as shown in 
Figures '5-29 and 5-30, the added notch filters do not improve this severe 
condition. . The notch pole damping was maintained at 0. 5, while the damp- 
ing of the zeros was 0. 1 and 0. 2, respectively. The frequency of both 
notches was kept at Y = 0. 174 rad/sec (S-plane frequency of 8. 8 rad/sec). 

The third ROI guidance system design utilized Atlas rate gyros that 
were relocated to station 708. This location was selected to give more 
weight to the- bending stability problem over the first two -thirds of 
the flight. A plot of the first mode antinode movement with flight time is 

I- 

shown in Figure 5-31, where, at approximately t = 110 sec, the antinode 
moves...pas.t station 708. Assuming that the relocation is acceptable from 
installation considerations, the corresponding ga-in-phase plots for the 
liftoff, maximum q, and BECO conditions are given in Figures 5-32 
through 5-34, and a complex pole filter at Y = 0. 3 rad/sec (15. 5 rad/sec 
in the S-plane) is also included. Considerable improvement in the stability 
portrait is evident for the liftoff and maxim-urn q cases. For the BECO 
condition, phase stability of the first mode is removed when the antinode 
passes the rate gyro station. To regain phase stability as well as to 
obtain gain stability at the time of gain change (T = 110 sec), a second 
complex pole filter can be added to produce a good control system design. 

The use of filter changes has been avoided in past Atlas designs 
because of the electronic complications. This may explain the selected 
location of the AC- 13 rate gyros. However, use of the digital computer 
as part of the control system permits filtering changes to be made with 
ease.. This advantage may be used to improve the control system perfor- 
mance margins over those existing in the current Atlas design. Reloca- 
tion of the rate gyros to station 70 8 may not be practical, but placement 
of this -gyr.o package at station 675, as in previous designs, will undoubt- 
edly be acceptable. The rather ample stability margins earlier in 
flight will be reduced slightly while easing the phase shifting effects 
later inflight. 
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Figure 5-29. Gain-Phase Plot Atlas /Centaur AC- 13 Digital Autopilot 
T=80 Spacecraft Sensors Compensation 1 
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Figure 5-30. Gain-Phase Plot Atlas/ Centaur AC-13 Digital Autopilot 
T=80 Spacecraft Sensors Compensation 2 
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Figure 5-3 i. Atlas/Centaur AC-13 First Mode Antinode Location 
versus Flight Time 
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Figure 5-33. Gain-Phase Plot Atlas /Centaur AC-13 Digital Autopilot 
T=80 Atlas Rate Gyros at Station 708 
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Figure 5-34. Gain-Phase Plot Atlas/ Centaur AC- 13 Digital Autopilot 
T=B0 Atlas Rate Gyros at Station 708 
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5. 4 STABILITY ANALYSIS OF THE SATURN IB/ CENTAUR VEHICLE 

The first stage stability of a Saturn IB/ Centaur vehicle carrying a 
4100-lb Jupiter probe was examined. Since the S~IB first stage consists 
of eight 70 -in. diameter tanks clustered about a 105 -in. diameter center 
tank, considerable structural rigidity of this stage results. Bending 
stability is easily achieved and location of the rate sensors either in the 
aft S-TVB compartment, instrument unit, or spacecraft is acceptable 
although the aft S-IVB location produces the largest stability margins. 

The data employed in the analysis is given in app. D. The forward- 
loop integrators previously included in the Atlas / Centaur control system 
are not used on the Saturn vehicles since they employ dynamic pressure 
feedback servovalves that have zero static leakage, and, hence do not 
induce attitude errors required to compensate for the hydraulic leakage 
flow. For the design with rate gyros located in the aft S-IVB compartment, 
the gain-phase stability for liftoff, maximum q, and burnout conditions of 
the S-IB stage are given in Figures 5-35 through 5-3 7, respectively. The 
compensation filter used consisted of double real poles at a W-plane fre- 
quency equivalent to 12. 5 rad/ sec in the S -plane. Ample stability margins 
are shown in the figures and listed in Table 5-1. Similarly, the stability 
margins, with the spacecraft rate sensors used corresponding to Figures 
5-38 through 5-40 are listed in Table 5-H. The rigid-body gain margins 
for both designs are shown in Figure 5-4ia, with the desired position 
gain schedule indicated. 

The corresponding rigid-body phase margins are shown in Figure 5-41b 
as a function of flight time. The rigid-body phase margin is showndropping 
briefly to 0. 22° as a result of the gain change; however, this is considered 
acceptable based on previous vehicle design experience. The closed-loop 
propellant slosh roots listed in Table 5-III indicate an unstable mode due 
to the slosh masses in the S-IB 70-in. LOX tanks. The slosh instability 
is small as indicated by the real part of the root, the inverse of which pro- 
duces a 700 -sec divergence time constant. Because of the large magnitude 
of the time constant and since a similar condition undoubtedly exists and is 
acceptable on Saturn IB/Apollo flights, it must be assumed to be acceptable 
for the Saturn IB/ Centaur vehicle. However, further evaluation of the pro- 
pellant slosh oscillations resulting from the instability would be warranted. 
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Figure 5-35. Gain-Phase Plot Saturn IB/ Centaur at T=0 
With Compensation i 
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Figure 5-36. Gain-Phase Plot Saturn IB/ Centaur Max Q T=80 
With Compensation 1 
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Figure 5-37. 


Gain-Phase Plot Saturn IB/ Centaur at T=B0 
With Compensation i 
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Figure 5-38. Gain-Phase Plot Saturn IB/Centaur T=0 Position 
Sensor at Spacecraft Compensation 1 
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Figure 5-39. 


Gain-Phase Plot Saturn IB/ Centaur T=80 Position 
Sensors at Spacecraft Compensation 1 
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SATURN. IB/CENT T = BO POS SENSORS AT S/C COMPENSATION 1 
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Figure 5-40. 


Gain-Phase Plot Saturn IB/Centaur T=B0 Position 
Sensors at Spacecraft Compensation 1 
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TABLE 5-1 


SATURN IB/ CENTAUR CONTROL SYSTEM GAINS AND RELATIVE 
STABILITY MARGINS - AFT S-IVB RATE GYROS 


Flight 

Time 

K R^D 

(sec) 

k d 

Rigid 

Body 

Phase 

Margin 

(deg) 

1st 

Bend 

Phase 

Margin 

(deg) 

1st 

Bend 

Ampli- 

tude 

above 

Nominal 

Gain 

(dB) 

Digital Filter D(Z) ' 
and Comment 

0 

1.2 

1.413 

40 

- 60 

8. 5 

1 1 







* (1-1.2Z + 0. 36Z ) 

80 

1.2 

1.413 

38 

-98 

5 








S -plane equivalent: 

110 

Gain change 

from Kp- 

= 1.413 to 


7 


k d = 

0.562 



(12. 5 ) 

153 

1.2 

0. 562 

38 


-10 

(s+12. 5) 2 


TABLE 5 -II 


SATURN IB /CENTAUR CONTROL SYSTEM GAINS AND RELATIVE 
STABILITY MARGINS - SPACECRAFT SENSORS 


■ 

k d 

Rigid 

Body 

Phase 

Margin 

(deg) 

1st 

Bend 

Phase 

Margin 

(deg) 

1st Bend 
Amplitude 
above 
N ominal 
Gain 
(dB) 

Digital Filter D(Z) 
and Comment 

0 

1. 189 

38 

-53 

15 

6i ( l - fr 2 " 1 ) 

(i -4 z_1 +|- z ) 

80 

1.. 189 

36 

-82 

14 

Equivalent K^/Kj^ S 

110 

Gain ch= 

k d = o. 

mge from ¥ 
562 

: D = i-i8 

9 to 

term is incorporated 
in D(Z), 

S-plane equivalent: 

153 

0.562 

36 

-130 

2 

100 (1.2 s+1) 
(s+10) 2 
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Figure 5-41a. Saturn IB/Centaur Rigid-Body Gain Margins 






(a) Booster Rate Gyros at Station 950 
(AFT S-IVB) 



Figure 5-41b. Rigid-Body Phase Margin versus Flight Time 
With Position Gain as a Parameter for the 
Saturn IB /Centaur Stage I Digital Control System 


5 . 5 ANALYSIS OF THE VOYAGER SPACECRAFT 
5 . 5 ., 1 Powered Flight Stability Analysis 

The Voyager spacecraft design consists of two vehicles; 1) the bus 
which performs midcourse corrections and Martian orbit deboost firings, 
and 2 ) the capsule which subsequently separates from the bus and lands 
on the Martian surface. 
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TABLE 5 -III 


SATURN IB /CENTAUR CLOSED LOOP SLOSH POLES, S-IB PHASE 


t 

0 

80 

153 

S-IB Fuel 

0 ±j4.838 

0 ± j6.032 

-0.003 ± j3 .288 

S-IB 70-in. LOX 

+0.0015 ±j4.863 

+0.0004 ± j6.1 12 

- 

S-IB 105-in. 
LOX 

-0.0005 ± j 3 .97 3 

-0.0003 ± j4.931 

- 

S-IVB LOX 

-0.069 ± j3.174 

-0.075 ± j3.948 

-0.088 ± j4.3l6 

S-IVB LH2 

-0.0669 ± j2 .544 

-0.075 ± j3.184 

-0.0713 ± j3 .496 

Centaur LOX 

-0.158 ± j3.631 

-0.128 ± j4.157 

-0.105 ± j4.839 

Centaur LH2 

-0.099 ± j3.860 

-0.107 ± j4.821 

-0.106 ± j5.246 


The LM Descent Stage Engine is employed at two thrust levels, 
that is, a high level for the Mars Orbit Insertion (MOI) firings; and a 
throttled level for Midcourse Correction (MCC) and Mars Orbit Trim 
(MOT) firings. The engine is gimbaled and controlled by electrome- 
chanical actuators for pitch and yaw axis stability while pulsing of jets 
are employed for roll control. 

The powered flight stability analysis performed on the Voyager 
spacecraft used the data given in app. E and included the propellant 
slosh effects for both the bus and capsule modes. The primary body 
bending modes were those caused by coupling of the bus with various 
appendage resonances, the bus essentially acting as a rigid beam with 
little deformation. The major appendage modes were those caused by 
an aft equipment module, a high-gain antenna, the cantilevered capsule, 
and a platform (PSP) employed for scanning the planet surface. The true 
bus bending modes were of much higher frequency and safely ignored. 

The analysis was previously performed for an analog control sys- 
tem design; however, very little difference with a digital control system 
is obtained if the high sampling' rate of 25 samples/ sec employed for 
booster control is also maintained during the spacecraft powered phases. 
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The use of forward loop integrators to reduce trajectory errors 
because of thrust vector misalignments was found desirable and is 
discussed in par. 5.4.2. - 

Gain-phase stability portraits for the four mission phases studied 
are given in Figures 5-42 through 5-45. The low-frequency crossover 
due to inclusion of an attitude integrator (gain of 0. 125) is indicated. A 
large but stable bus slosh loop and a smaller capsule slosh loop are 
visible followed by the aft equipment module mode, high gain antenna 
mode, capsule, and PSP modes. The desired position gain for the low- 
thrust midcourse correction and Mars orbit trim firings is shown to be 
Kp = 4, while the desired gain during the Mars orbit-insertion phase is 
kept constant at = 0.5. No compensation filters were included in the 
design; however, use of some filtering would be desirable to further 
increase stability margins and furnish added noise attenuation. 

5. 5. 2 Voyager Spacecraft Thrust Vector Pointing Error Analysis 

To estimate the thrust vector control (TVC) pointing errors and 
crossrange velocity errors incurred during midcourse correction (MCC), 
Mars orbit insertion (MOI), and Mars orbit trim (MOT), the pointing 
errors introduced by the gyros, the attitude control system, vehicle limit 
cycles, and attitude reference misalignments were combined with thrust 
vector pointing errors due to eg offsets, and thrust vector misalignments 
to obtain the total powered flight thrust vector pointing error. The 
effects of three control system designs on the thrust-vector -induced 
pointing error were evaluated. Following are the three control system 
designs studied: 

a) Without integrators 

b) With attitude feedback integration 

c) With attitude and engine angle feedback integration 

The larger eg uncertainty and the small distance between the engine and 
eg locations at MOT produce a large offset angle. Attitude and engine 
angle feedback integration reduce the steady-state pointing error for 
these effects to zero. The addition of an outer guidance loop (closed-loop 
guidance steering) will significantly reduce these errors. 
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Figure 5-42. Gain-Phase Plot, MCC, with Capsule, Integral 
Attitude Feedback, Kj = 0. 125 
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j- VOYAGER START MOI W/CAP S. A. STRUT THETA INT A = 0 Kl = 0. 125 



Figure 5-43. Gain-Phase Plot, Start MOI, With Capsule, 
Integral Attitude Feedback, Kj. = 0. 125 
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Figure 5-44. Gain-Phase Plot, End MOI, With Capsule, 
Integral Attitude Feedback, = 0. 125 
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Figure 5-45. Gain-Phase Plot, MOT, With Capsule, 


Integral Attitude Feedback, Kj. = 0, 125 
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Time histories of the ungraded control system instantaneous point- 
ing error) and velocity error, V,p, for a 1° thrust vector misalign- 

ment or a 1° thrust vector offset angle are shown in Figures 5-46 through 
5-51 for the three integrator cases. For constant thrust levels, the 
velocity error curve can be obtained from the pointing error curve by 
multiplying the thrust by the integral of . 



Figures 5-46. Without Integrators 



10 20 30 40 50 

TIME (SEC) 


Figures 5-47. Without Integrators 
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Figures 5-48. e,^, With Attitude Feedback Integration 



TIME (SEC) 

Figures 5-49. With Attitude Feedback Integration 
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Figures 5-50. e^, With Attitude and Engine Angle Feedback Integration 



TIME (SEC) 


Figures 5-51. V-j> With Attitude and Engine Angle Feedback Integration 


81 



The use of attitude feedback integration with biases for mean eg 
offsets is considered more desirable than the use of attitude plus engine 
angle feedback integration. Although the latter design removes the 
effects of eg offsets and thrust vector offsets, it does so rather slowly, 
thereby compromising its usefulness for this purpose. Also, difficulties 
in its implementation result since the engine feedback loop is unstable 
without main engine thrusting. The attitude feedback integration design 
does not exhibit this instability; hence, the feedback loop can be closed 
prior to engine startup. 

A 380-sec main engine firing (AV = m/sec) was-assumed for MCC - 
and a 330-sec firing (AY = 2 km/sec) for MOI. The total pointing error 
angle and the cross velocity error at the end of engine firing are provided 
for each of the control systems. The results are based on a constant 
coefficient simulation. Since the mass and inertia properties change ■ 
radically from the start to the end of MOI, these results are considered 
preliminary. 

It was determined that the use of integrators; considerably improved 
the steady-state total pointing angle and cross-velocity accuracy. For 
each of the unguided control systems, the total pointing angle errors, 
after reaching steady state, were as follows: 



Unguided Control Systems Total 
Steady-State Pointing 
Angle Error (°) 


MCC 

MOI 

MOT ■ 


- Mean 

3c r 


3o" 

Mean 

3c r 

Without integrators 

0. 28 

0.87 

1.5 

1.52 

0. 5 

2.04 

With attitude feedback inte- 
gration, no bias angles 

0. 22 

0.73 

0.5 

0. 67 

1.0 

1.98 

With attitude feedback inte- 
gration and with bias 
angles 

0 

0.73 

0 

0. 67 

0 

1.98 

With attitude and engine 
angle feedback integration 

0 

0.76 

0 

0.76 

0 

0.76 
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APPENDIX A 


CONTROL SYSTEM STABILITY CRITERIA 

i 

The control system stability criteria followed in this study are given 
in this appendix. In general, the stability margins employing the ROI 
guidance and control system were required to be equal to-.or greater than-' 
the margins existing with present vehicle designs. In addition, the 
specific stability requirements given below were employed. 

1. RIGID-BODY STABILITY MARGINS 

The minimum rigid -body stability margins generally used are 6 dBi. 
of gain and 30° of phase. Decreases in margins below these minimqms 
often occur for a short duration because of required gain changes. ■ Some 
relaxation in the margin requirements based upon the length of this 'dura- 
tion may be warranted since deficiencies for a few seconds would be 
tolerable at rigid-body frequencies. The curves shown in Figure A-l and 
A-2 illustrate the rigid-body stability requirements that would be repre- 
sentative of this more, flexible approach. 

2. PROPELLANT-SLOSH STABILITY MARGINS 

Propellant-slosh stability is strongly influenced^ by the propellant 
damping available in the propellant tank, particularly through the use of 
tank baffles. The damping values employed in linear stability analyses are 
generally based upon acceptable slosh-limit cycles amplitudes. The'limit- 
cycle amplitude reaches a steady-state condition when slosh stability mar- 
gins become zero, hence, zero margins could be considered as being the- 
minimum requirements for linear slosh stability. If these requirements 
are not met, divergence to a larger slosh limit cycle would be expected. J 
If the divergence is slow and occurs over a short duration, it may still be; 
acceptable; however, a simulation study to determine this would be re- *" 
quired. For the purposes of this study, zero slosh margins are considered 
as minimum requirements, 

3. STRUCTURAL BENDING MODE- STABILITY MARGINS 

Phase stabilization of the first and second structural bending modes 
is generally accepted, particularly if the mode frequencies are well below 
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(less than one-third) the actuation system frequency. Since nonlinearities 
within the actuation system produce significant phase uncertainty at high 
frequencies, phase stabilization at these frequencies is generally avoided. 

A criteria to be employed for phase stabilization should therefore be based 
upon the frequency of bending stabilization. An example of such criteria 
that can be used is given in Figure A-3 in which the allowable-phase margin 
requirements are related to the ratio of bending mode and rigid-body 
frequencies. For gain stabilized bending modes a gain margin criteria as 
exemplified in Figure A-4 could be employed. 



RATIO OF BENDING MODE FREQUENCY 
TO RIGID-BODY FREQUENCY 


Figure A-3. An Example of Acceptable Phase Margins for Phase 
Stabilized Bending Modes as a Function of Modal 
Frequency Ratio 
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ACCEPTABLE PHASE MARGIN (dB) 



RATIO OF BENDING MODE FREQUENCY 
TO RIGID-BODY FREQUENCY 

Figure A-4. Acceptable Gain Margins for Gain Stabilized 
Bending Modes as a Function of Modal 
Frequency Ratio 

For the launch vehicles considered in this study, phase stabilization 
of all first bending modes and the Saturn V vehicle second bending mode 
has been employed as criteria in the past. Gain stabilization of all other 
bending modes was a further requirement in these programs. These 
criteria are also considered applicable to the vehicle configurations used 
in this study. 
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APPENDIX B 


VEHICLE DYNAMICS AND CONTROL SYSTEM 
HOMOGENEOUS EQUATION MATRIX 

The generalized homogeneous equation matrix representing the 
vehicle and control system dynamics used by TRW for linearized control 
system analyses (see Section 6) is given in this appendix. In addition to 
the rigid body equations of motion, including aerodynamics, this repre- 
sentation includes fuel sloshing and bending (flexible body)’ dynamics. 

Table B-i and Figure B-l present the nomenclature used and, the defini- 
tion of vehicle and control system angles used in the analysis. Figure B-2 
is a generalized block diagram of the control system. Table B-2 contains 
the homogeneous equation matrix. 

In all propellant slosh data given, a spring-mass model was em- 
ployed for the propellant slosh modes as indicated by -Figure B-l. 
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Table B-l 


Definition of Variables and Coefficients 

Variables of the Matrix 

5 Laplace Operator 

0 Vehicle attitude angle from reference line 

Aqj Displacement of first- stage oxidizer slosh mass 
A^ Displacement of first-stage fuel slosh mass 
Aq 2 Displacement of second- stage oxidizer slosh mass 
Displacement of second-stage fuel slosh mass 
q^ Displacement of normalized first-bending mode 

q^ Displacement of normalized second-bending mode 

q^ Displacement of first- stage third-bending mode 

q^_ Displacement of normalized fourth-bending mode 

<5 Engine displacement angle 

a Vehicle total angle of attack 

y Vehicle velocity angle from reference line 

ijj Vehicle velocity angle of attack without winds 

a' Angle of attack due to winds 

V^ Wind velocity 

V^ Vehicle velocity with respect to air 

<5 . Engine command input to the actuation system (control system 
loop broken at this point to enable open-loop stability analysis 
and insertion of additional compensation) 

Q Total sensed gyro feedback from both position and rate gyros 

B 

5 co Engine command outputs from the autopilot 

6 Total autopilot error signal prior to autopilot compensation 
6 

f' c Command angle (shown with pseudoguidance but need not be 

used) 
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Table B-l 


Definition of Variables and Coefficients (cont'd) 


Coefficients of the Matrix'-' 


I 

M 

M. 

M 

M 


01 


fl 


02 


f2 


K 


01 


K 


fl 


K 


o2 


K f2 
i 01 
l £l 
*02 
*f2 
^01 
^fl 
^02 
^f 2 


to 

co 

co 

CO 

,T. 

T 


01 


fl 


02 


f 2' 


Vehicle moment of inertia 

Slosh mass of first- stage oxidizer 

Slosh mass of first-stage fuel 

Slosh mass of second-stage oxidizer 

Slosh mass of second-stage fuel 

Slosh spring constant of first- stage oxidizer 

Slosh spring constant of first- stage fuel 

Slosh spring constant of second-stage oxidizer 

Slosh spring constant of second-stage fuel 

Slosh mass moment arm of first- stage oxidizer 

Slosh mass moment arm of first- stage fuel 

Slosh mass moment arm of second-stage oxidizer 

Slosh mass moment arm of second- stage fuel 

Slosh mode damping of first-stage oxidizer 

Slosh mode damping of first- stage fuel 

Slosh mode, damping of second oxidizer 

Slosh mode damping of second-stage fuel 

Slosh mode frequency of first- stage oxidizer 

Slosh mode frequency of first-stage fuel 

Slosh mode frequency of second- stage oxidizer 

Slosh mode frequency of second- stage fuel 

Total vehicle thrust 

Vehicle control thrust 


positive 
values if 
. slosh 
masses are 
aft of c. g. 
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Table B-i 


V = 
M = 


m 


e 


1 

e 


1 

x 


*h 





0 ). 


CO. 


CO, 



<P 2T 

<P 3T 


Definition of Variables and Coefficients (cont'd) 

Aerodynamic normal force 
Aerodynamic normal force moment arm 
Vehicle inertial velocity 
Vehicle mass excluding slosh masses 
Total control engine mass 

Effective distance of engine c.g. from engine gimbal 
Distance between engine gimbal and vehicle c. g. 

Total control engine inertia about gimbal 
Normalized mass of the first-bending mode 

i 

Normalized mass of the second-bending mode 

Normalized mass of the third-bending mode 

Normalized mass of the fourth-bending mode 

Total vehicle axial acceleration including aerodynamic 
drag effects 

Modal frequency of the first-bending mode 
Modal frequency of the second-bending mode 
Modal frequency of the third-bending mode 
Modal frequency of the fourth-bending mode 
Damping of the first-bending mode 
Damping of the second-bending mode 
Damping of the third-bending mode 
Damping of the fourth-bending mode 

Displacement of the engine gimbal for the first-bending mode 
Displacement of the engine gimbal for the second-bending mode 
Displacement of the engine gimbal for the third-bending mode 
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Table B-i 


4T 


K 


IT 


2T 


y 3T 

04 T 

n G 

*2G 

0 3G 

^4G 

0 1R 

0 2R 

^3R 

0 4R 

*011 


0fU 

0 O2l 


0 f21 

0 O12 

0 fl2 

0 O22 

0 f22 


Definition of Variables and Coefficients (cont'd) 


Displacement of the engine gimbal for the fourth-bending mode 

Bending slope at the engine gimbal for the first-bending mode 

Bending slope at the engine gimbal for the second-bending mode 

Bending slope at the engine gimbal for the third-bending mode 

Bending slope at the engine gimbal for the fourth-bending mode 

Bending slope at the position gyro for the first-bending mode 

Bending slope at the position gyro for the second -bending mode 

Bending slope at the position gyro for the third-bending mode 

Bending slope at the position gyro for the fourth-bending mode 

Bending slope at the rate gyro for the first-bending mode 

Bending slope at the rate gyro for the second-bending mode 

Bending slope at the rate gyro for the third-bending mode 

Bending slope at the rate gyro for the fourth-bending mode 

First-bending mode displacement at the first- stage oxidizer 
slosh station 

First-bending mode displacement at the first- stage fuel 
slosh station 

,First-bending mode displacement at the second-stage oxidizer 
slosh station 

First-bending mode displacement at the second-stage fuel 
slosh station 

Second-bending mode displacement at the first- stage oxidizer 
slosh station 

Second-bending mode displacement at the first fuel 
slosh station 

Second-bending mode displacement at the second oxidizer 
slosh station 

Second-bending mode displacement at the second fuel 
slosh station 
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Table B-l 

Definition of Variables and Coefficients (cont’d) 


0 Third-bending mode displacement at the first oxidizer 

slosh station 

0 Third-bending mode displacement at the first fuel 

slosh station 

&Q23 Third-bending mode displacement at the second oxidizer 
slosh station 

0 Third-bending mode displacement at the second fuel 

slosh station 

^014 Fourth-bending mode displacement at the first oxidizer 
slosh station 

Fourth-bending mode displacement at the first fuel 
slosh station 

^024 Fourth-bending mode displacement at the second oxidizer 
slosh station 

0f24 Fourth-bending mode displacement at the second fuel 
slosh station 

= Control system autopilot position gain 
= Control system autopilot rate-to-position gain 
= Guidance loop gain 

= Pseudoguidance pole (need not be used) 

c^, b^, b^, a^, a^ = Coefficients for modelling the actuation 
system and control compensation. 


K 


K. 


R 


IG 


' 2 ’ 


Note: The structural bending mode equations require a positive dis- 

placement and slope occurring at the engine gimbal point. 
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VEHICLE 

CENTERLINE 



Figure B-l. 


Definition of Vehicle and Control System Angles 



Figure B-2. Control System Block Diagram 
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Table B~2 

Vehicle Dynamics and Control System Homogeneous Equation Matrix 
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APPENDIX C 


ATLAS/ CENTAUR CONTROL SYSTEM DATA 
1 . INTRODU C TION 

Appendix C contains' data necessary for the analysis of the Atlas/ 
Centaur control system. Figure C-l shows the general vehicle configura- 
tion, including the location of control- system components. Table C-l 
defines the nomenclature used in describing the vehicle dynamics and the 
Atlas/Centaur control system. 

Sections 2, 3, and 4 of this appendix present aerodynamic param- 
eters, fuel-sloshing data, and bending data for this vehicle. Section 5 
presents data and linearized models for the Atlas (booster and sustainer 
stages) and Centaur thrust-vector control (TVC) systems. 


196 



197 
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Figure C-l. Atlas /Centaur General Configuration 
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Table C-l 

Atlas/Centaur Nomenclature 

Vehicle mass excluding propellant slosh masses 

Vehicle moment of inertia 

Vehicle longitudinal c. g. station 

Control thrust 

Total thrust 

Vehicle acceleration 

Control thrust moment arm 

Control moment coefficient 

Vehicle velocity 

Aerodynamic normal force 

Aerodynamic moment arm 

Aerodynamic moment coefficient 

Atlas oxidizer slosh mass moment arm 

Atlas fuel slosh mass moment arm 

Centaur oxidizer slosh mass moment arm 

Centaur fuel slosh mass moment arm 

Atlas oxidizer slosh mass moment arm 

Atlas fuel slosh mass moment arm 

Centaur oxidizer slosh mass moment arm 

Centaur fuel slosh mass- moment arm 

Atlas s oxidizer slosh frequency 

Atlas fuel slosh frequency 

Centaur oxidizer slosh frequency 

Centaur fuel slosh frequency 
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Table C-l 

Atlas/ Centaur Nomenclature (con't) 


“ao 

*af 

^co 

^cf 

"l 

W 2 

W 3 

M 1 

M 2 


fc l 
^ 2 

^ 3 

0 


T1 


T2 


T3 


0 1 

0' 


T1 


T2 


0’ 


T3 


0 ' 

0 ' 

0 > 


lap 

2ap 

3ap 


Atlas oxidizer slosh mode damping 
Atlas fuel slosh mode damping 
Centaur oxidizer slosh mode damping 
Centaur fuel slosh mode damping 
First-bending mode frequency 
Second-bending mode frequency 
Third-bending mode frequency 
First-bending mode mass 
Second-bending mode mass 
Third- bending mode mass 
First-bending mode damping 
Second-bending mode damping 
Third- bending mode damping 

First- bending mode deflection at booster engine gimbal 
station 

Second-bending mode deflection at booster engine gimbal 
station 

Third-bending mode deflection at booster engine gimbal 
station 

First-bending mode slope at booster engine gimbal station 

Second-bending mode slope at booster engine gimbal 
station 

Third- bending mode slope at booster engine gimbal 
station 

First-bending mode slope at Atlas position gyro station 
Second-bending mode slope at Atlas position gyro station - 
Third-bending mode slope at Atlas position gyro station 


199 - 



Table C-l 

Atlas/ Centaur Nomenclature (con't) 


0i 

0t 

0t 

0t 

0i 

0> 

0. 


1 sp 
2sp 
3sp 
lar 
2ar 
3ar 
Irr 


0 

0 


2rr 

3rr 


M 

e 


w TwD 

U DwT 

K 

c 

K ' 
c 


6 

Sc 


First-bending mode slope at spacecraft gyro station 

Second-bending mode slope at spacecraft gyro station 

Third-bending mode slope at spacecraft gyro station 

First-bending mode slope at Atlas rate gyro station 

Second-bending mode slope at Atlas rate gyro station 

Third- bending mode slope at Atlas rate gyro station 

First-bending mode slope at recommended rate gyro 
station 

Second-bending mode slope at recommended rate gyro 
station 

Third- bending mode slope at recommended rate gyro 
station 

Booster engine inertia about gimbal 
Booster engine mass 

Booster engine mass moment arm about gimbal 

Booster engine gimbal station 

Booster engine tail wags dog frequency 

Booster engine actuator frequency 

Booster engine actuator loop velocity gain 

Booster engine actuator loop model gain for low-frequency 
approximation. 

Booster engine deflection angle 
Booster engine command angle 
Booster engine actuation rate 
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2. AERODYNAMIC, MASS PROPERTIES, AND TRAJECTORY- 
RELATED CONTROLS DATA 

Aerodynamic parameters, mass properties data, and trajectory- 
related control parameters for the Atlas /Centaur vehicle are given in 
Table C-2 for the booster portion of flight. Additional trajectory-relhted 
data is given in References C-l and C-2. 

Table C-2 

Aerodynamic, Mass Properties, and Trajectory- 
Related Controls Data for Atlas/Centaur 

Aerodynamic Parameters 


t sec 

0 

40 

80 1J_2 

154 

lb/ rad 

2. 14 

49. 2 x 

10 3 300 X 10 3 134 X 10 3 

11. 8 x 10 

i ft 
P 

42. 3 

42. 2 

37. 9 33. 8 

13. 0 

-2 

u„, sec 
ra. 

0 

0, 070 

3.55 1.56 

0. 075 

Mach 

0 

0.41 

1.65 3.61 

8. 07 

Q lb/ft 2 

0 

193 

801 320 

28 

Mass Properties 





M slugs 

d 

9530 

7310 

5350 3790 

1920 

I slug ft 2 

4. 25 x 10 6 

3. 60 x 

10 6 3. 2 x 10 6 2. 9 x 10 6 

2. 05 x 10 

X in. 
eg 

766 

787 

788 762 

610 

Trajectory and Control Parameters 



t ib 

c 

330 x 10 3 

343 xlO 3 368 xlO 3 378 xlO 3 

382 x 10 3 

T T lb 

389 x 10 3 

407 x 10 3 446 xlO 3 460 xlO 3 

465 x 1 0 3 

A^t ft/ sec 2 

37. 5 

49. 6 

66.9 99.2 

183 

i ft 
c 

37. 2 

35.4 

35.3 37.5 

50. 2 

u sec' 2 
c 

2. 88 

3. 38 

4. 05 4. 88 

9. 33 

V ft/ sec 

3 

456 

1600 . 3570 

8530 

Atlas position gyro station 

_ 

991 in. 


Atlas rate gyro 

station 

= 

584 in. 


Spacecraft gyro 

station 

= 

100 in. 


Recommended rate gyro station = 

708 in. 
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3-. PROPELLANT SLOSHING DATA 


Propellant- sloshing data for the Atlas /Centaur vehicle is given 
in Table C-3 and in Figures C-2 through C-4. 

Table C-3 

Propellant Slosh Parameters for Atlas / Centaur Vehicle 


Propellant Slosh Parameters for Atlas/ Centaur Vehicle 


t sec 

0 

40 

80 

112 

154 

c 

ao ft 

l 

CO 

CO 

-13.4 

-3. 2 

4.4 

17. 5 

1 , ft 
af 

14. 0 

18. 6 

24. 3 

29. 0 

42. 5 

l ft 

CO 

-28. 8 

-31. 0 

-31. 0 

-28. 5 

-15. 9 

l ft 
cf 

-40. 5 

-42.4 

-42.4 

-40. 1 

-27. 5 

M slugs 

100 

374 

374 

374 

330 

slugs 

75 

280 

280 

280 

180 

M slugs 

210 

210 

210 

210 

210 

M ^ slugs 

18 

18 

18 

18 

18 

co rad/ sec 

ao 

4. 68 

4. 27 

4. 96 

6. 04 

7. 32 

co af rad/ sec 

4. 68 

4. 27 

4. 96 

5. 96 

6. 88 

co rad / sec 

co 

5. 59 

6.43 

7. 47 

9. 09 

12.40 

Wc £ rad/ sec 

3. 83 

4. 41 

5. 12 

6. 24 

8. 48 

^ao N * 

0. 0009 

0. 0009 

0. 0010 

0. 0080 

0. 0004 

*af N. D. 

0. 0008 

0. 0008 

0. 0010 

0. 0030 

0. 0008 

t N. D. 
co 

0. 0030 

0. 0030 

0. 0030 

0. 0030 

0. 0030 

£ c f N " D ‘ 

0. 0008 

0. 0008 

0. 0008 

0. 0008 

0. 0008 
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Figure C-2. Slosh Frequencies versus Flight Time 



Figure C-3. Slosh Masses versus Flight Time 



Figure C-4. Slosh Moment Arms versus Flight Time 
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4. VEHICLE BENDING DATA 


Vehicle bending data for the Atlas/ Centaur vehicle is given in 
Table C-4 for three modes at five times of flight during booster flight. 
Additional data will be found in Reference C-2. 

Table C-4 

Bending Mode Parameters 


Bending Mode Parameters 


t sec 

0 

40 

80 

112 

154 

Wj rad/sec 

13. 46 

13. 92 

15. 58 

18. 52 

25. 89 

w., rad/ sec 

32. 37 

37. 29 

38. 45 

38. 71 

45. 97 

Wj rad/ sec 

46. 37 

47. 07 

47.88 

48. 38 

49. 53 

slugs 

3171 

2740 

2085 

1617 

2166 

slugs 

3257 

5044 

10620 

8462 

3087 

Mj slugs 

2122 

1989 

2360 

3249 

12015 

N. D. 

0. 005 

0. 005 

0. 005 

0. 010 

0. 010 

S z N. D. 

0. 005 

0. 010 

0. 010 

0. 010 

0. 010 

£ 3 N. D. 

0. 010 

0. 010 

0. 010 

0. 010 

0. 010 

P T1 ft/ ft 

1. 0 

1. 0 

1. 0 

1. 0 

1. 0 

«> T2 ft/ft 

1. 0 

1. 0 

1. 0 

1.0 

1. 0 

?> T3 ft/ft 

1. 0 

1. 0 

1. 0 

1. 0 

l. 0 

«>' T1 rad/ft 

0. 052 

0. 055 

0. 061 

0. 071 

o. no 

^'t 2 rad / ft 

0. 124 

0. 169 

0. 181 

0. 189 

0. 456 

0- T3 rad/ft 

0. 368 

0.413 

0.470 

0. 511 

0. 403 

V i rad/ft 

0. 045 

0. 046 

0. 049 

0. 054 

0. 070 

®' 2ap rad/ft 

0. 051 

0. 046 

0. 049 

0. 051 

0. 043 

^'Sap rad/£t 

0. 0017 

-0. 0078 

-0. 015 

-0. 021 

-0. 056 

*'ls P rad/ft 

-0. 075 

-0. 063 

-0. 061 

-0. 067 

-0. 141 

P' 2sp rad/ft 

0. 125 

0. 200 

0. 328 

0. 273 

0. 084 

0' 3sp rad/ft 

-0. 067 

-0. 051 

-0. 036 

-0. 018 

-0. 066 

*'lar rad/ft 

-0. 033 

-0. 031 

-0. 026 

-0. 023 

-0. 022 

?' 2ar rad/ft 

-0. 045 

-0. 062 

-0. 107 

-0. 097 

-0. 064 

♦'3ap rad/ft 

0. 040 

0. 026 

0. 0131 

0. 012 

0. 063 

*' lrr raa/ft 

-0. 001 

-0. 001 

-0. 003 

0. 002 

0. 015 

0, 2rr rad/ft 

-0. 072 

-0. 1 

-0. 122 

-0. 096 

-0. 059 

P '3rr rad/ft 

-0. 0139 

-0. OH 

0, 001 

0. 004 

0 04 
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5. THRUST VECTOR CONTROL DATA 

Data and linearized models for the Atlas/ Centaur TVC system are 
given in Tables C-5 through C-7 and in Figures C-5 and C-6. 

Table C-5 

Booster Engine Servoactuator Data 


Booster Engine Parameters For One Engine 

2 


: h 

slug- ft 

= 

377 

M 

e 

slug 

= 

30. 8 

i 

ft 

_ 

2. 52 

e 

x h 

in. 

= 

1212 

W TWD 

rad/sec 

= 

46 

U DWT 

rad/sec 

= 

93 

K 

sec - ^ 

- 

36 

C 

6 

deg 

= 

5 

max 

^a max 

deg/ sec 

= 

28 


Engine actuation model for 6 to = 2° /sec: 


K to 
c DWT 


S~ " (S+27.5) (S + 68. 35 ± j81. 6) 
c 


Low-frequency approximation for 5co - 2° /sec: 

18 


t K 

o c 


c S -5- K 


S + 18 


Engine actuation model for = l°/sec: 


K to 
c DWT 


6 ~ (S + 14. 5) (S + 120 ± j84) 


Low-frequency approximation for 6 to = 1 /sec: 

I 

K 


S-f K 


12. 5 
S + 1 2. 5 
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Table C-6 

Atlas Servoactuator Characteristics 


Symbol 

Definition 

Units 

Booster 

Su stainer 

A - - 

Actuating piston area 

- ft 2 

2 ; 47 x 10~ 2 ■ 

4. 1 x 10" 2 

B 

Bulk modules of hydraulic fluid 

lb /ft 2 

3. 89 x 10 7 

3. 89 x 10 7 

C B 

Coulomb friction coefficient 

ft -lb 

5. 65 x 10 2 

5. 59 x 10 2 

C L 

Discharge coefficient for leakage 
'(bypass) orifice 

ft 2 / sec 
lb/ft2 ' 

4.78 x 10 -6 

2.47 x 10“ 6 

C V 

Viscous friction coefficient 

ft-lb-sec/deg 

' 6. 3 x 10 

3.76 x 10 

I R 

Moment of inertia, engine -gimbaled "wet" 
mass about engine gimbal axis 

slug -ft 2 

3. 77 x 10 2 

4. 20 x 10 2 

K 

a 

Servoamplifier gain 

ma/V 

2.7 

3.42 

K 
. m 

Spring constant of actuator -backup 
structure 

* lb/ft 

1. 138 x 10 6 

1.441 x 10 6 

K t 

Feedback transducer gain 

volt/ deg 

2. 2 

3. 0 

K v 

Servovalve discharge coefficient 

ft 2 / sec 
mA 'J lb/ ft 2 

7. 125 x 10“ 6 

1.68 x 10“ 6 


Distance from center of gravity of 
gimbaled mass to center line of 
engine gimbal axis 

ft 

2. 52 

2.68 
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Table C-6 


Atlas Servoactuator Characteristics (Cont'd) 


Symbol 

Definition 

Units 

Booster 

Sustainer 

■ m r 

Gimbaled mass, one engine "wet" 

slugs 

3. 08 x 10 

3. 19 x 10 

P R 

Hydraulic return pressure from 
servovalve 

lb /ft 2 

7. 2 x 10 3 

7. 2 x 10 3 

p s 

Hydraulic supply pressure to servovalve 

lb/ ft 2 

4. 32 x 10 5 

4.32 x 10 5 

R 

Distance from center line of actuator 
mounting to center line of engine at 
gimbal point 

ft 

1.769 

8. 33 x 10 -1 

V T 

Total volume of hydraulic fluid in one 
actuator cylinder 

ft 3 

7.6x10 3 

4. 14 x 10 -3 

T 

Thrust 

lb 

1. 545 x 10 5 

5. 7 x 10 4 

K 

c 

No-load, open-loop actuator velocity gain 


36. 8 

18. 06 

6 

max 

Engine displacement limit 

deg 

5 

3 

• 

‘-’max 

Engine actuation rate limit 

deg/sec 

28 

- 
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Table C-7 


Centaur Servovalve -Actuator -Engine Parameters- 
RL-10A-3 Propulsion System 


Symbol 

Definition 

Units 

Value 

A 

Actuating piston area 

ft 2 

1.06 x 10~ 2 

B 

Bulk modules of hydraulic fluid 

lb/ft 2 

3. 89 x 10 7 

C b 

Coulomb friction coeffifcient 

ft /lb 

200 

C 

V 

Viscous friction coefficient 

ft lb/ deg/ sec 

3. 62 


Moment of inertia, engine -gimbaled "wet" mass 
about engine gimbal axis 

slug -ft 2 

66. 5 

m r 

Gimbaled mass, one engine "wet" 

slugs 

11. 0 


Distance from center of gravity of gimbaled mass 
to center line of engine gimbal axis 

ft 

2. 04 

R 

Distance from center line of actuator mounting to 
center line of engine at gimbal point 

ft 

1.16 

K 

a 

Servoamplifier gain 

mA/ volt 

3.6 

K t 

Feedback transducer gain 

volt/ deg 

2. 06 

K 

■m 

Spring constant of actuator -backup structux*e 

lb/ ft 

5. 1 x 10 5 

K 

V 

Servovalve discharge coefficient 

. (ft 3 / sec) 
mA y/ lb/ ft 2 

1. 02 x 10“ 6 
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Table C-7 


Centaur Servovalve -Actuator-Engine Parameters- 
RL-10A-3 Propulsion System (Cont'd) 


Symbol 

Definition 

Units 

Value 

P R 

Hydraulic return pressure from servovalve 

lb /ft 2 

7.2 x 10 3 

p s 

Hydraulic supply pressure to servovalve 

lb/ft 2 

1. 512 x 10 5 

T 


lb 

1.5 x 10 4 

V T 

Total volume of hydraulic fluid in one 
actuator cylinder 

ft 3 

1. 33 x 10" 3 

K 

c 

No-load, open-loop actuator velocity gain 

1/sec 

13. 1 

a f 

Area of washout piston DPF 

ft 2 

3.41 x 10" 4 

C F 

Discharge coefficient for DPF orifice 

(ft 3 / sec) 
/ lb /ft 2 

1. 342 x 10“ 7 

K f 

Spring constant-washout piston 

lb/ ft 

9. 36 x 10 3 

K 

s 

Pressure feedback gain 

mA/ lb/ ft 2 

1. 218 x 10~ 4 




Figure C-6. Linearized Block Diagram Centaur Servoactuator 
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APPENDIX D 


SATURN IB/ CENTAUR CONTROL SYSTEM DATA 
1. INTRODUCTION 

Appendix D contains data necessary for the analysis of the 
Saturn IB/ Centaur control system. Figure D-l shows the general 
vehicle configuration, including the location of control-system com- 
ponents. Sections 2, 3, and 4 of this appendix present the aerodynamic 
parameters, fuel- sloshing data, and bending data for this vehicle. 
Section 5 presents data and a linearized model for the Saturn IB thrust 
vector control system. Data for the Centaur TYC system is given in 
Appendix C. 
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VEHICLE STATION (IN.) 


Figure D-l. Saturn IB/ Centaur General Configuration 








2. AERODYNAMIC, MASS PROPERTIES, AND TRAJECTORY-RELATED 
CONTROLS DATA 

Aerodynamic parameters, mass properties data, and trajectory- 
related controls parameters are summarized in Table D-i for the 
Saturn IB/ Centaur vehicle. .Additional aerodynamic data {normal force 
coefficient) is given in Figures D-2 and D-3 for the 260-in. and 154-in. 
diameter fairings. 'Mass properties data as a function of vehicle weight 
and flight time is given in Figures D-4 through D-6. Additional trajec- 
tory data may be found in Reference D-l. 

3. PROPELLANT SLOSHING DATA 

Propellant sloshing data for the Saturn IB/ Centaur vehicle is given 
in Figures D-7 through D-9. 

4. VEHICLE BENDING DATA • 

Vehicle bending data for the" Saturn IB/ Centaur vehicle is given in 
Table D-2 at three times of flight during th'e S-I stage burning and in 
Table D-3 at three times of flight during the S-IVB stage burning. Addi- 
tional data is given in Reference D-2. 


213 



-214 


Table D- 1 


Saturn IB/ Centaur , Summary of Mass Properties, Aerodynamic Parameters, 
and Trajectory-Related Control Parameters 


''pan^ 

S-IB Flight 

S-IVB Flight 

L. O. 

80 

153 

158. 6 

200 W/S* 

200 WO/S* 

625. 6 

Mass 

Properties 








M a 

(slugs) 

39, 831. 

25, 202 

12,423 

8954 , 

8406 

8142 

1880 

I 

(slug-fm) 

44. ’9x10° 

43. 6x1 0 6 

22. 4x10° 

3. 25x10 6 

3. 22x1 0 6 

2. 5 8x10° 

1. 22x10 6 

X cg 

! (in. ) 

718 

724 

1069 

1340 

1349 

1328 

1553 

Trajector and Control 








Parameters 








T c 

• (lb) 

0. 74x1 0^ 

0. 83xl0 8 

0. 6 22x10 6 

2x10 5 

2x10^ 

2x10;? 

2x1 0 S : 

T 

(lb) 

1. 47xl0 6 

1. 66x1 0 6 

1. 24x10 6 

2x10 5 

2x1 0 5 

2x1,0 5 

2x10° 

A T r 

(ft/sec 2 ) 

36. 9 

59. 0 

99. 8 

22. 4 - 

23. 8 

24. 6 

106. 3 


(ft) 

51. 5 

52. 0 

80. 7 

103. 0 

104. 0 

102. 0 

121. 1 

' ^c 

(1 / sec^) 

0. 85 

0. 99 

2. 78 





'V 

(ft/ sec) 

1. 0 

1275. 0 

5924. 0 

5817 

6036 

6036 

24197 

gcosy (ft/sec 2 ) 

32. 2 

29. 63 

20. 08 

19- 54 

15. 71 

15. 71 

0 

Aerodynamic 








Pararrieters 








N a 

; (lb/rad) 

'0 

1. 4xl0 6 

12. 9x1 0 3 






__ 

*p 

• (ft). ■ • 

11. 9 

15. 3 . 

17. 8 

— 

— 


— 

}ia 

(1 / sec 2 ) 

0 

0. 48 

0. 010 

_ — 

_ — 

_ y. 



"Q 

(lb/ft 2 ) 

o 

565. 0 

7. 9 

— 

— 

— 

~ - 

Mach 

0 

1. 317 

5. 837 

- ~ 

— 

- - 

— 


' Shroud 


-1, RAO AND CP/D - CALIBERS FROM STA 100 


STA 



Figure D-2a. Saturn IB/Centaur Three-Stage Vehicle 
260” Diameter Fairing 



Figure D-2b. Normal Force Coefficient! Saturn IB/ Centaur 
260" Diameter Fairing 
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\ Figure D-3b. Normal Force Coefficient Saturn IB/Centaur 
154" Diameter Fairing 
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Horizontal C. G. (int x 10 -i ) (Saturn Coordinate System) 





S-IB FLIGHT TIME (SEC) 



Versus Weight, Liftoff to S-I Burnout 






S-IVB FLIGHT TIME (SEC) 

0 50 inn 150 ?nn ?5 n 500 550 400 450 500 550 



Weight (lb x 10" 3 )' 

■Figure D-5.'. Moment of Inertia and C. G. Versus Weight, 
S-IVB Start to S-IVB Burnout 
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Figure D-6. Moment of Inertia and C.G. Versus Weight, Centaur Separation to Centaur Burnout 
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Figure D-8. Slosh Masses Versus Flight Time 
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Table D-2 

Saturn S-I Stage Bending Mode Parameters 
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Table D-2 

Saturn S -I Stage Bending Mode Parameters (con’t) 
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Table D-3 

Saturn S-IVB Stage Bending Mode Parameters 


ro 

DO 


" Flight Time 

158. 6 

200 sec 

200 sec 

625 

Bending """ 

Mode Parameters 

sec 

w/ Shroud 

wo/ Shroud 

sec 

Combined Mode Bending Freq. 





0 )^ (rad/ sec) 

92. 27 

99. 60 

35. 59 

41. 31 

0)2 (rad/sec) 

197. 03 

238. 97 

109. 89 

172. 71 

0)3 (rad/sec) 

■ 260 . 92 

272. 08 

181. 22 

228. 97 

Mass (Slugs) 





Mi 

M- 

8954 

8406 

8142 

1880 

8954 

8406 

8142 

1880 

M 3 

8954 

8406 

8142 

1880 

Damping (N. D. ) 





I 1 

s 2 

0 . 01 

0 . 01 

0 . 01 

0 . 01 

. 0 . 01 

0 . 01 

0 . 01 

0 . 01 

h 

0 . 01 

0 . 01 

0 . 01 

0 . 01 

Modal Deflection (ft/ft) 





?T1 At Station 

JT2 x = 1186 

HT3 

0 . 90 
0 . 82 
0. 36 

0 . 90 
0. 38 
0 . 22 

0. 45 
0. 535 
0. 53 

1 . 2 
2 . 0 
2 . 0 











w 

ro 

Ui 


Table D~3 

Saturn S-IVB Stage Bending Mode Parameters ( cont ' d ) 



200 sec 
w / Shroud 

200 sec 
wo/ Shroud 

625 

sec 

0. 163 

0. 0711 

0. 0533 

0. 125 

0. 156 

0. Ill 

0. 238 

0. 2 

0. 122 

-0. 0945 

0. 0522 

0. 05 

-0. 0511 

-0. 178 

0. 0222 

0. 200 

-0. 104 

-0. 089 

-0. 105 

-0. 32 

-0. 145 

0. 23 

0. 32 

0. 34 

0. 36 

-0.70 

0. 33 








5. THRUST VECTOR CONTROL DATA 


The linearized TVC Model for the Saturn IB control system is 
given in Figure D-10 and the associated data for the S-I and S-IVB stages in 
Table D-4. Table D-5 gives the corresponding TVC transfer functions. 



Figure D-iO. Linearized Block Diagram of Saturn iB 
Thrust Vector Control System 
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Table D-4 

Parameters for S-IB and S-IVB Stages 


Parameter 

S-I * 

(H- 1 Engines) 

S-rVB 

( J-2 Engine) 

I 

n 

(slug-ft 2 ) 

1015. 0 

2660 

d 

(in. ) 

27. 5 

11. 625 

m l 

lb-sec 

in 

16. 1 

117. 0 

1 

n 

(ft) 

2. 657 

2.425 

K 1 

( 1/ sec) 

17. 142 

18. 76 

K 2 

( in/ lb) 

0. 704437 x 10~ 5 

0. 6573 x 10 -5 

K 

o 

(lb/ in) 

4. 35 x 10 5 

12. 889 x 10 5 

k e 

(lb/ in) 

0. 668 x 10 5 

3. 9101 x 10 5 

B 

(lb- sec/ in) 

116. 0 

538 

w 2 

(rad / sec) 

20.0 

12. 0 


H-l Engine data is only for one engine. 


Table D-5 

Thrust Vector Control System Transfer Function 


TVC 

System 

f(s> 

c 

Approximated (s) 

c 

S-I 

S ’S 

_ + 1 
Uo + L J 

(^> + 0 

( S , lV 3 V S , lA 

( S I lV s ! ill 

1,20.674 J^16.770 ^5. 8073 ±j 127.45 j 

( 20. 674 + l6. 770 + 1 j) 

S-IVB 

s 

-~+ 1 
V 12 / 

(u + 1 ) 


f s . . V s . ,"V s . i\ 




Note: The effect of engine reaction is neglected. 


S-IB actuator limit ±8° S-IVB actuator limit ±7° 

actuation rate limit ±24°/ sec actuation rate limit ±8 /sec 
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APPENDIX E 


SATURN V (VOYAGER PAYLOAD) CONTROL SYSTEM DATA 
1. INTRODUCTION 

Appendix E contains data necessary for the analysis of the Saturn V 
control system and the Voyager spacecraft. Figure E-l shows the gen- 
eral Saturn V vehicle configuration with the Voyager payload. Section 3 
presents the fuel sloshing data for this vehicle configuration. Vehicle 
bending data is given in Reference E-l. Section 3 presents data on the 
Saturn V thrust vector control (TVC) system for the S-IC, S-II and S-IVB 
stages. Controls data fox - the Voyager spacecx-aft is presented in Sec- 
tion 4^ 
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Figure E-l, 


i i. 

Saturn V (Voyager Pay- \ 
load)* Vehicle Configuration ' 
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Figure E-3. Slosh Masses Versus Flight Time 



Figure E-4. Slosh Mass Stations Versus Flight Time 
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NORMALIZED AMPLITUDE RATIO (DB) 


3 - SATURN V THRUST VECTOR CONTROL SYSTEM -DATA 

3. 1 SATURN V CONTROL COMPUTER FREQUENCY RESPONSE 

The Saturn V control computer mixing amplifier has a frequency 
response-as shown in Figure E-5. The linear equation that describes 
this response is as follows: 

J o 1 666 , 

I. 

m. 


L n = Input current in mA 
I = Servovalve or output current in mA 


(396)' 


1 . 70 
(396) 


S + 1 


20 


10 


-10 


-20 


-30 


-40 


1 1 1 i — i 1 1 1 

i i iiii n 

i 1 i nn . 



AMPLITUDE RATIO 






PHASE LAGV “ 
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Figure E-5. Saturn V Control Computer Mixing 
Amplifier Frequency Response 
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3- 2 s-IC STAGE TVG DATA 

The linear block diagram for the S-IC stage TVC system is shown in 
Figure E-6. The nonlinear block diagram and system parameters are 
shown in Figure E-7 and Table E-l, respectively. The frequency response 
and time response specifications are shown in Figures E-8, E-9, and 
E- 10. 

The linear equation describing the S-IC Stage thrust vector control 
system is as follows: 



where: 

/3 is in degrees 
j3 is in milliamperes 
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0.04 in- lb/ ma • 

2940 cis/ in- lb 
3.32 x 10“^ sec ' • 

57 in 2 '.,.- : . 

205.000 lbs/in . 

239.000 lbs/ in, . ;; 


94.5 lb-sec 2 /in 


-- 505 Ib-sec/in 

in -lb 

6.03688 x 10"^. psi 


0.53 sec 


— 0.25 sec 


f . . 


t— 0.3479. in-lb/in 


6e max = 5 * 15 de § 
^ e max = 7 de S^ Sec 
R =5.31 ft 

n 

I n = 56,760 slug ft 

w DWT = “n = 38 rad / sec 
U TWD = ^ rad / sec 


Figure E-6. S-IC Stage Thrust Vector Control System 
Block Diagram and System Parameters 





































Table E-l 
S-IC Parameters 


k tm 

= 

0. 04 -P' ~ Ib 
ma 

K 

V 

= 

2940 

m. -lb 

T 

V 

= 

3.32x10 ^ sec 

T d 

= 

0. 25 sec 

H 1 

= 

0. 3479 ~ lb 

m. 

H 2 

= 

3. 69 x 10 -5 ~~ 

psi 

H 3 

= 

2. 51 x 10 -5 ~ 

psi 

H 4 

- 

5. 116 in - - lb 
m. 

A a 

- 

57 in. 2 

M 

e 

= 

94. 5 lb -. sec2 
m. 

D 

- 

505 lb : sec 
in. 

p s 

= 

1800 psi (nom) 

k l 

- 

239, 000 ~ 
m. 

k t 

= 

205, 000 ~ 

m. 

(X ) max, 
a 

= ± 0. 05 in. 

< q s2 » 

max. = ±740 cis 

.(^) max. 

= ±5. 74 in. 

F 

g 

= 

12, 000 lb 


= 

66, 000 + 30, 000 


Q, 


Q 


s2 


Q c 


Q 


s 2 


-(AP - P) 
s 

I AP - P 


(AP + P) 
s_ 

I AP - P I 


AP 


- 1 


AP 


1/2 


1/2 


+ 1 
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FREQUENCY (Hz) 


Figure E-9. S-IC Frequency Response Specification (phase lag) 
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POSITION (DEG) 



TIME (SEC) 


Figure E-10. S-IC Time Response Specification 
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3.3 S-II STAGE TVC DATA 

The linear block diagram for the S-II Stage TVC system is shown in 
Figure E-ll. The nonlinear block diagram and system parameters are 
shown in Figure E-12 and Table E-2, respectively. The frequency response 
specification is shown in Figure E-13. 

The linear equation describing the S-II Stage thrust vector control 
system is as follows: 


^e 

0. 159 [0. 08 S + l] 

i 

r s i 

11. 02 + 1 

. 

f 


r s i 

— + j 

421 


S 2 2(0. 399)S 

2 33 22 + 1 , 

(33.73) ' 

\ 


1 


S 2(0, 987 )S 

j (94. 22) 2 94.22 

wher e : 

/? is in degrees 

/? c is in milliamperes 
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MeS + DS 


in-lb I ?d K dpf S 


T d S +-1 


in- lb 


6.035 in - lb/ma 

K, 312,000 lb/ in 


220 cis/in - lb 

M 123 lb “ sec 2 / in ' 

* e 

<5e =7 dee 

max K 

3 <,18 x 3 . 0" 3 sec 

■q 625 lb - sec/ in 

^max = 9l 6 de 8 /sec 


' • ' 

R = 0. 991 ft 

13 in 2 '.' 

. . - . v ‘ — = 0.08 sec 

• • a. ■ ' 

n 

260,000 Tb/in 

’ “O 

- - 7 . , .L. 10 -„< 10 . in - lb'/psl 

„ dpf - . . . ' 

I n = 2690 slug ft 

> 

•. K-s; 1.06 in lit /in 

“dwt = “n = 34 ra -d/ sec 


31 , , ' • » ’ 

«Twd “ rad/ sec 

Figure E- 1 1. 

S-II Linear Thrust Vector Control System 
Block Diagram and System Parameters 





































Table E-2 

S-II TVC System Block Diagram Parameters 


k tm 

- 

no r " 1^5 

0. 035 a 

mA 

K ' ’ 

V 

= 

220 . C1S ^ 
in. -lb 

T 

V 

= 

3. 18x10 ^ sec 

T d 

= 

8 x 10”^ sec 

H 

= 

1.06“- - lb 
m. 

A 

a 

= 

. 2 
13 m. 

M 

e 

= 

,, 2 
. „ ^ ib-sec 

125 : 

in. 

D 

= 

625 lb .- sec 
m. 

P 

s 

= 

3500 psi (nom) 

k l 

= 

312, 000 lb /in. 

K t 

- 

260, 000 lb/in. 


(X^) max. = ±0. 26 in. 
(Q 9 ) max. = ±71 cis 
(/3^) max. = ±1. 51 in. 


g 


= 9000 lb 

= 20, 000 + 40, 000 sin (500 t) 

1/2 


°5 

(AP + P s ) 

AP 

Q s2 

|AP + P s | 

P + 1 

s 

Q 4 

-(AP - 

AP , 

Q s2 

AP - P s | - 

— - 1 
s 


1/2 
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Figure E-13. 1 S-II Frequency Response Specification 
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3.4 S-IVB STAGE TVC DATA 


The linear block diagram for the S-IVB Stage TVC system is shown 
in Figure E-14 The nonlinear block diagram and system parameters are 
shown in Figure E-15 and Table E-3, respectively. The frequency 
response specification is shown in Figure E-16. 

The linear equation describing the S-IVB Stage thrust' vector 
control system is as follows: 


o. 1496 [ 0.08 S + l] 


32. 257 


+ 1 


10. 49 


+ 1 


405 


+ 1 


J L 


X 


174 


+ 1 


[(34. 68) 


_ + 2(0. 327) s 
2 + (34. 68) + 


where: 


is in degrees 
is in milliamperes 
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125 . lb'^sec'V in 
538. lb-sec/ln 

0.08 sec. 
i r\ ' 1 ’*1 r\"*5' .. „ it 



R = 0.966 ft 
n 


; K , 'I0 x "io: 5 '-in-ib/psi V I n = 2660 slu § ft 


1.129 in-lb/in . . “ W T = “n" 31 rad/ sec 
' f • * “> TWD = 24 rad/ sec 


Figure E-14. S-IVB Linear Thrust Vector Control System 
Block Diagram and System Parameters 














in. -lb. 


H 



Figure E-15. S-IYB Thrust Vector Control Syst 
Nonlinear Block Diagram 
















Table E-3 

S-IVB TVC System Block Diagram Parameters 


k tm 

= 

„ m. -lb 

0. 035 . 

mA 

K 

V 

- 

149 5 c * s 
1 in. -lb 

T 

V 

= 

3.18x10 2 sec 

T 

d 

= 

_2 

8x10 sec 

H 

- 

1. 129 in '-— b 
m. 

A 

a 

= 

11. 78 in. 2 

M 

e 

= 

„ 2 
, lb -sec 

l/LD. 1 

m. 

D 

= 

538 lb ; sec • 
m. 

P 

s 

= 

3650 psi (nom) 

k l 

' - 

284, 000 lb/in. 

k t 

= 

233, 000 lb/in. 

(X^) max. 

= 

±0. 26 in. 

(Q _) max. 
s ^ 

= 

±40. 4 cis 

(/?.) max. 

= 

±1 . 45 in. 

F 

g 

= 

9000 lb 


q 5 (AP + 


15, 000 + 30, 000 sin (500 t> 

jV At , , 1/2 

P P ^ X 
s s 


Q 4 -(AP 


- P ) 
s 
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amplitude ra - 


20 


40 



FREQUENCY (Hr) 


Figure F-I6. S-IVB Frequency Response Specification 
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4. VOYAGER SPACECRAFT DATA 


Control system terminology used in analyzing the Voyager control 
system is given in Table E-4. Values of the control system parameters 
are given in Tables E-5 and E-6 for two spacecraft configurations (capsule 
on and off) for four different times of flight during the interplanetary 
pha s e . 

The Voyager powered flight linear actuation control system model 
is given in Figure E-17 and the associated data in Table E-7, 
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Table E-4 

Voyager Analysis Data Terminology 


Mass Parameters 

2 

I = vehicle moment of inertia, slug ft 
M = vehicle mass, slugs 
Engine Parameters 

T = engine thrust, lb 

= engine mass, slugs 

= distance between engine c. g. , and engine gimbal, ft 

= distance between engine gimbal and vehicle c. g. , ft 

2 

Ir = inertia of engine about the gimbal point, slug ft 
Slosh Parameters 


L BO’ 

L BF’ 

L CO’ 

l 'cf 

M BO’ 

m bf 

M 

’ CO’ 

M 

CF 

K BO’ 

K BF’ 

K CO’ 

rr 

CF 

W BO’ 

CO 

BF’ 

W CO’ 

co CF 

^BO’ 

^BF’ 

^CO’ s 

CF 


Bending Parameters 

M r M 2 , M 3 , M 4 

• i0 v u> 2 , co 3 , w 4 


1’ S 

2’ h 

’ ^4 


IT’ 

<P ZT , 

0 3T ’ 

*4T 

IT’ 

0 2T’ 

0 3T’ 

^4T 


bus and capsule, fuel and oxidizer, 
slosh mass distances from vehicle 
c. g. , (positive aft of c. g. , towards 
gimbal), ft 


slosh masses, slugs 
slosh mode spring constants, lb -ft ^ 
slosh mode frequencies, rad/sec 
slosh mode damping, no dimensions 


bending mode masses (four modes), 
slugs 

bending mode frequencies, rad/sec 

bending mode damping, no dimensions 

bending displacements at the gimbal 
point (positive values), ft 

bending slopes at the gimbal point 
(positive values), ft/ft 
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Table E-4 

Voyager Analysis Data Terminology (con't) 

0Boi > 0 BO2’ 0 BO3’ 0 BO4 ~ bendin S displacements at bus 

oxidizer slosh mass station, ft 


0 BFl’ 0 BF2’ 0 BF3’ 0 BF4 


bending displacements at bus 
fuel slosh 


0 rQr 0 ro2 > ^C03’ 0 CO4 = bend ing displacements at capsule 

oxidizer slosh mass station, ft 


0 1 G’ 0 2G’ 0 3G’ 0 4G 


= bending slopes at the position 
gyro station, ft /ft 


Variables 


rigid body vehicle attitude angle, rad 
bus oxidizer slosh mass displacement, ft 
bus fuel slosh mass displacement, ft 
capsule oxidizer slosh mass displacement, ft 
capsule fuel slosh mass displacement, ft 
first bending mode amplitude, rad 
second bending mode amplitude, rad 
third bending mode amplitude, rad 
fourth bending mode amplitude, rad 
engine displacement angle, rad 
engine command angle, rad 
attitude error angle, rad 

integrator feedback attitude command, rad 


Station Numbers 


Engine gimbal point 

Bottom of bus propellant tanks 

Gyro location 


111 in. 


118. 25 in. 


136 in . 


Capsule slosh mass attach points 
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Table E-5 


Voyager Spacecraft Data 
Parameter Values, Vehicle With Capsule On. . 



MCC 

Start 

MOI 

End 

MOI 

MOT 

I, slug ft 2 

26, 157 

26, 623 

19, 094 

•19, 094 

M, slugs 

636 

595 

286 

286 

T,' lb 

1, 050 

7, 750 

7, 750 

1, 050 

M g , slugs 

12.4 

12. 4 

12. 4 

12.4 

L, , ft 
e 

1. 88 

1. 88 

1. 88 

1. 88 

l t , ft 

4. 91 

4. 91 

4. 9L. 

4. 91 

I H , slug ft 2 

50 

50 

‘ 50 

50 

L BO = ^BF’ ft 

1. 68 

1. 68 

5. 8 

5. 8 

L CO = L CF’ ft 

-9. 15 

-9. 15 

-7. 35 

-7. 35 

m bq > slugs • 

75. 6 

75. 6 

6. 23 

6. 23 

M BF’ s1u s s 

47. 2 

47. 0 

3. 87 — 

3. 87 

M co , slugs 

2. 7 

2. 7 

2. 7 

2. 7 

M cf , slugs^ 

0. 9 

0. 9 

0. 9 

0. 9 

K bo > lb -ft 

105 

840 

76 

10. 4 

l^BF’ lb “ft 

66 

555 

47 

6. 4 

K CO' lb “ft 

9.7 

78 

160 

22 

k cf , lb -ft 

3. 2 

26 

53 

7. 4 

"’BO = W BF J rad/sec 

1. 18 

3. 35 

3. 48 

1. 29 

a> C Q, rad/sec 

1. 89 

5. 37 

7. 69 

2. 87 

co CF , rad/sec 

2. 30 

6. 55 

9. 40 

3. 49 

£ BO = £ BF 

0. 001 

0. 001 

0. 001 

0. 001 

^CO = ^CF 

0. 01 

0. 01 

0. 01 

0. 01 
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Table E-5 


Voyager Spacecraft Data 

Parameter Values, Vehicle With Capsule On (con't) 



MCC 

Start 

MOI 

End 

MOI 

MOT 

M 1 = M 2 = M 3 
= M 4 , slug ft 2 

1 

1 

1 

1 

o)j, rad/sec 

24. 70 

24. 70 

24. 95 

24. 95 

10 ^, rad/sec 

24. 97 

24. 97 

25. 16 

25. 16 

rad/sec 

37. 50 

37.50 

39- 21 

39- 21 

rad/sec 

47. 13 

47. 13 

47. 89 

47. 89 

^1 = ^2 = £ 3 = ^4 

0. 01 

0. 01 

0. 01 

0. 01 

0 rjn > 

0. 00506 

0. 00506 

0. 0117 

0. 0117 

0 2T’ ft 

0. 00216 

0. 00216 

0. 00431 

0. 00431 

0 3T’ ft 

0. 00347 

0. 00347 

0. 00958 

0. 00958 

04T ’ ^ 

0. 00262 

0. 00262 

0. 00706 

0. 00706 

0 1T 

0. 000592 

0. 000592 

0. 001089 

0. 001989 

0 2T 

0. 000480 

0. 000480 

0. 00070 

0. 00070 

0 3T 

0. 000254 

0. 000254 

0. 00068 

0. 00068 

^4T 

0. 000543 

0. 000543 

0. 00094 

0. 00094 

0 BO1 = - 0 BFl’ ft 

0. 00353 

0. 00353 

0. 00885 

0. 00885 

0 BO2 = 0 BF2’ ft 

0. 000092 

0. 000092 

0. 0025 

0. 0025 

0 BO3 = 0 BF3’ ft 

0. 00282 

0. 00282 

0. 00781 

0. 00781 

0 BO4 = 0 BF4’ ft 

0. 00122 

0. 00122 

0. 0046 

0. 0046 

0 CO1’ 0 CF1 ’ ft 

-0. 00332 

-0. 00332 

-0. 00366 

-0. 00366 

0 CO2’ 0 CF2’ ft 

-0. 00466 

-0. 00466 

-0. 0056 

-0. 0056 

0 CO3’ 0 CF3’ ft 

-0. 000025 

-0. 000025 

0. 00013 

0. 00013 

0 CO4’ 0 CF 4’ ft 

-0. 00525 

-0. 00525 

-0. 00648 

-0. 00648 
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Table E-5 


Voyager Spacecraft Data 

Parameter Values, .Vehicle With Capsule On (con't) 




MCC 

Start 

MOI 

End 

MOI 

MOT 

*1G 


0. 000595 

0. 000595 

0. 00109 

0. 00109 

^2G 


0. 00048 

0. 00048 

0. 000070 

0. 000070 

?> 3 G 


0. 00025 

0. 00025 

0. 000068 

0. 000068 

-0 

■s. 


0. 00054 

0. 00054 

0. 000094 

0. 000094- 

03'p 


0. 00129 

0. 00129 

0. 00172 

0.00172 

^BOl 

= 0 BF1 , ft 

0. 0058 

0. 0058 

0. 01112 

0. 01112 

0 BO2 

= ^BF2’ £t 

0. 00348 

0. 00348 

0. 00107 

0. 00107 

0 BO3 

= 0 BF3’ ft 

0. 00382 

0. 00382 

0. 00118 

0. 00118 

0 1G 


0. 000626 

0. 000626 

0. 00112 

0. 00112 

0 2G 


0. 00159 

0. 00159 

0. 00137 

0. 00137 

0 3G 


0. 00129 

0. 00129 

0. 00172 

0; 00172 

MCC 

duration 

« 380 sec for Mars arrival 

separation 


MOI 

duration 

~ 380 sec 




MOT 

duration 

=s undetermined- 
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Table E-6 


Voyager Spacecraft Parameter Values, Vehicle with Capsule Off 



MCC 

Start 

MOI 

End 

MOI 

MOT 

I, slug -ft 2 

14, 500 

13, 523 

4, 547 

4, 547 

M, slugs 

550 

500 

208 

208 

T, lb 

1, 040 

7, 750 

7, 750 

1, 040 

M^, slugs 

12. 4 

12. 4 

12. 4 

12. 4 

L , ft 
e 

3. 19 

3. 05 

2. 5 

2. 5 

l t , ft 

3. 19 

3. 05 

2. 5 

2 5 

I H , slug ft 2 

50 

50 

50 

50 

L BO “ ^BF’ 

0. 995 

1. 142 

1. 71 

1. 71 

M Bq> slugs 

75. 6 

75. 6 

6. 23 

6. 23 

M BF’ - s1u § s 

47. 2 

47. 0 

3. 87 

3. 87 

K bq , slugs 

121 

997 

104 

14 

K bf , slugs 

76 

624 

63 

8. 7 

“bo = “bf* rad / sec 

1. 26 

3. 65 

4. 07 

1. 50 

S BO = £ BF 

0. 001 

0. 001 

0. 001 

0. 001 

M =M 2 = M 3 = M 4 

1 

1 

1 

1 

u)^, rad/sec 

25. 32 

25. 32 

25. 32 

25. 32 

< o 2 > rad/sec 

25. 76 

25. 76 

25. 76 

25. 76 

u> 3 , rad/sec 

49. 27 

49. 27 

49. 27 

49. 27 

^1 = ^2 = ^3 

0. 01 

0. 01 

0. 01 

0. 01 

»1T- £t 

0. 00741 

0. 00741 

0. 0140 

0. 0140 

0 2T > ft 

0. 00678 

0. 00678 

0. 00462 

0. 00462 

d 3T 5 ^ 

0. 00118 

0. 00118 

0. 00563 

0. 00563 

*1T 

0. 000626 

0. 000626 

0. 00112 

0. 00112 

0 2T 

0. 00159 

0. 00159 

0. 00137 

0. 00137 
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Figure E-17, Block Diagram of Voyager Powered Flight 
Control System Model 
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